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ABSTRACT 


Deep dielectric charging and subsequent electrostatic discharge in composite materials 
used on spacecraft have become greater concerns since composite materials are being used more 
extensively as main structural components. Deep dielectric charging occurs when high energy 
particles penetrate and deposit themselves in the insulating material of spacecraft components. 
These deposited particles induce an electric field in the material, which causes the particles to 
move and thus changes the electric field. The electric field continues to change until a steady 
state is reached between the incoming particles from the space environment and the particles 
moving away due to the electric field. An electrostatic discharge occurs when the electric field is 
greater than the dielectric strength of the composite material. The goal of the current 
investigation is to investigate deep dielectric charging in composite materials and ascertain what 
modifications have to be made to the composite properties to alleviate any breakdown issues. A 
1-D model was created. The model is given the space environment calculated using the 
Environmental Workbench software, the composite material properties, and the electric field and 
voltage boundary conditions. The output from the model is the charge density, electric field, and 
voltage distributions as functions of the depth into the material and time. Analysis using the 
model show that there should be no deep dielectric charging problem with conductive 
composites such as carbon fiber / epoxy. With insulating materials such as glass fiber / epoxy, 
Kevlar, and polymers, there is also no concern of deep dielectric charging problems with average 
day-to-day particle fluxes. However, problems can arise during geomagnetic substorms and 
solar particle events where particle flux levels increase by several orders of magnitude, and thus 
increase the electric field in the material by several orders of magnitude. Therefore, the second 
part of this investigation was an experimental attempt to measure the continuum electrical 
properties of a carbon fiber / epoxy composite, and to create a composite with tailorable 
conductivity without affecting its mechanical properties. The measurement of the conductivity 
and dielectric strength of carbon fiber / epoxy composites showed that these properties are 
surface layer dominated and difficult to measure. In the second experimental task, the 
conductivity of a glass fiber / epoxy composite was increased by 3 orders of magnitude, 
dielectric constant was increased approximately by a factor of 16, with minimal change to the 
mechanical properties, by adding conductive carbon black to the epoxy. 
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CHAPTER 1 


INTRODUCTION 


The charging of spacecraft, and its possible role in spacecraft anomalies due to 
electrostatic discharges, has been a well known problem since the beginning of space 
flight some thirty years ago. Charging is caused by energetic particles in the space 
environment: electrons, protons and positively charged heavy ions. There are three types 
of charging: entire vehicle charging, surface charging, and internal charging, also known 
as deep dielectric charging. Entire vehicle charging is when the entire potential of the 
spacecraft is raised. Surface charging is when only the potential of the spacecraft surface 
is raised; however this may also occur locally where only part of the surface has its 
potential raised due to geometric and material considerations. Deep dielectric charging is 
like surface charging, except that the potential increase is not on the surface of the 
spacecraft component but inside the material of the component. The last two types are a 
concern for composite material structures, and the last type, deep dielectric charging, is 
the focus of this research. 

Most of the work done up until now on charging has been on surface charging. It 
is caused by low energy electrons (less than 30-50 keV) which do not penetrate the 
surface of the external structural material. These particles accumulate on the surface of 
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the structure, and can lead to large differential charging between various parts of the 
spacecraft. This charging can grow large enough to cause an electrostatic discharge, 
which can lead to surface damage or spacecraft anomalies. These spacecraft anomalies 
are caused by an induced internal current created by the discharge which can cause 
electronics damage or cause erroneous commands to be processed by the spacecraft. 

Deep dielectric charging occurs when high energy electrons or ions penetrate the 
surface of, and deposit charge within, a insulating material. If the deposition of incoming 
charged particles is greater than the charge leakage through the material, a large potential 
difference can build up in the material and lead to an electrostatic discharge. The 
discharge can occur within the material or from the interior of the material to one of its 
surfaces. In either case a current path is created and the material is locally damaged 
around the discharge site. There seems to be a correlation between these discharges and 
periods of increased flux, as spacecraft anomalies have been observed when there are 
more high energy particles in their local space environment. 

The goals of this research are to develop a model for the analysis of deep 
dielectric charging of insulating materials and composite materials, determine relevant 
composite electrical properties, and to develop technology to alleviate charging problems 
by changing the properties of the composite. A schematic of the model is shown in 
Figure 1.1. The model was used to create a Fortran Code, The Composite Deep 
Dielectric Charging Analyzer (CoDDCA). The code will have as inputs the space 
environment and material properties. The space environment is given in terms of electron 
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Figure 1.1 : Schematic of the deep dielectric charging model 
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and proton fluxes versus energy; also geomagnetic substorms and solar particle events 
may be specified. The material properties are divided up into mechanical, electrical, and 
geometrical properties. The outputs from the code are given as functions of depth into 
the material from the surface exposed to the space environment and time. They include 
the charge density profile, the electric field profile, and the voltage profile. 

The model was used to examine the effects of the various parameters input into 
the model. These sensitivity studies revealed which parameters are important to deep 
dielectric charging and which play minimal roles. The model was also used to examine 
specific case study orbits, including orbits where there have been suspected incidents of 
deep dielectric charging. 

Composites are treated as a homogenous material in the above model. Therefore 
composites on a micro-scale were modeled separately to determine the effects of their 
inhomogeneity on the electric field. The different arrangements investigated include the 
effects on the electric field due to a conductive fiber surrounded by an insulative matrix, 
unlikely fiber arrangements which can create points of electric field intensification, and 
the effects of an insulative surface layer of pure epoxy often found on top of conductive 
carbon fiber / epoxy composites. 

Experimental work was done to measure the continuum electrical properties of 
composites. Carbon fiber / epoxy composites were made of various thicknesses, areas, 
and layups and were used to investigate the through-thickness conductivity and dielectric 
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strength. Values typical of previous work were measured, but the conductivity and 
dielectric strength did not behave as continuum properties; they displayed non-classical 
thickness and surface layer dependencies. 

Another part of the experimental work was to create a conductivity tailorable 
composite. This work was undertaken because the sensitivity studies showed that the 
deep dielectric charging problem is dependent on conductivity, and that increasing 
conductivity can have a significant effect in reducing the likelihood of an electrostatic 
discharge. Therefore, a composite system was developed that had conductivity that could 
be tailored without greatly affecting its mechanical properties. A glass fiber / epoxy 
composite was used, the conductivity was increased by adding carbon black powder to 
the epoxy. Laminates made from this composite were investigated to determine the 
effects of the carbon black on the electrical and mechanical properties of the composite. 
Control of electrical properties over large ranges (3 orders of magnitude) without 
significant effects on mechanical properties, were achieved. 

Previous work relevant to the current research is described in Chapter 2. This 
includes a description of the space environment, analytical work on surface charging and 
deep dielectric charging, and some general information on composites. The problem 
statement and approach for the current research is presented in Chapter 3. The analytical 
methodology which is used in the Composite Deep Dielectric Charging Analyzer 
(CoDDCA) Code is developed in Chapter 4, as are micro-mechanical modeling of 
composite electrical properties, and the data reduction equations used in the experimental 
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analysis. The results of the code, including the parametric studies and the case studies are 
presented and discussed in Chapter 5. Chapter 6 describes the test matrices used, the 
procedures used to manufacture the composite laminates, and the experimental 
procedures used to measure the mechanical and electrical properties of the composite 
samples. The results of the experimental analysis are presented and discussed in Chapter 
7. Finally, conclusions and recommendations for future work are presented in Chapter 8. 
The appendices include a description of the CoDDCA code and how to use it, as well as 
the code output profiles and the experimental data. 
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CHAPTER 2 


BACKGROUND 


This chapter presents background information on the space environment, charging 
in the space environment, and composite materials. The space environment is broken up 
into different sources of radiation: the Van Allen radiation belts, galactic cosmic rays, 
solar particle events and geomagnetic substorms. The discussion of spacecraft charging 
in the space environment includes both surface charging and deep dielectric charging. 
Finally, a general description of composites and composite mechanics relevant to this 
work are given. 

2.1 RADIATION SOURCES 

The main sources of energetic particles in the space environment are trapped 
radiation of the Van Allen radiation belts, galactic cosmic rays consisting of 
interplanetary protons and ionized heavy nuclei, and particles associated with solar 
particle events and substorms. 1 2 The first two sources are fairly constant, while the third 
is highly time dependent. The Van Allen belt radiation commonly leads to deep 
dielectric charging, while radiation from the other sources is more likely to lead to single 
event phenomena (SEP) in electronic components. This occurs because the belts contain 
certain types of particle; the particle types determine the phenomena. The are two types 
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of single event phenomena, hard and soft, depending on if the damage which occurs is 
permanent or temporary. Soft single event phenomena are also known as single event 
upsets (SEU). They result in such problems as a change in stored data or wrong 
commands. Hard single event phenomena can result in permanent damage such as 
burnout of power MOSFET’s, gate rupture, latchup, or the freezing of bits. The simplest 
solution to single event phenomena is to use shielding. A major source of shielding is the 
spacecraft structure. However, not all wires and components can be easily shielded, and 
shielding can add significantly to the spacecraft mass. 

The Sun’s activity occurs in an 11 year cycle know as the solar cycle. The 
activity is quantified by the sunspot number (typically called the R value), the higher the 
number the more active the Sun is. Solar minimum defines the beginning of the cycle 
when the sunspot number is at its lowest. Solar maximum is when the Sun is most active 
and the sunspot number is at its highest, as shown in Figure 2.1. Each cycle can have 
maximums and minimums which vary by a factor of 4 from one cycle to the next. 

2.1.1 Van Allen Radiation Belts 

The Van Allen radiation belts consists primarily of energetic electrons and 
protons, with a smaller percentage of heavy ions like oxygen (0 + ). They are trapped by 
the Earth’s magnetic field lines, and thus gyrate about the poles along the magnetic field 
lines. They form toroidal belts around the Earth as shown in Figure 2.2. There are two 
belts of high concentrations of particles, the inner belt and the outer belt. The inner belt 
extends from approximately 1000 km to 6000 km in altitude, and is populated by very 
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Figure 2.2: 


Omnidirectional Electron Flux (#/cm 2 -s) 
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Electron and proton Van Allen radiation belts 4 
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high energy (10’s of MeV) protons and high energy (1-10 MeV) electrons. The outer belt 
extends from approximately 20,000 km to 32,000 km in altitude, and is populated 
primarily by high energy electrons. For this reason, most spacecraft have orbits with 
altitudes in the hundreds of kilometers, below the inner belt, and at geosynchronous orbit 
(GEO) which is above the upper belt, at 35,782 km in altitude. However, the Van Allen 
radiation belts continue beyond the high concentration belts out to an altitude of 60,000 
km, with smaller amounts of particles. It should be noted that when the particles travel 
near the poles, higher energy particles travel to lower altitudes producing the aurora 
borealis and aurora austrailialis. Therefore, spacecraft in low altitude polar orbits are 
subjected to the same density of high energy particles found in much higher equatorial 
orbits. 


2.1.2 Galactic Cosmic Rays 

Galactic cosmic rays are energetic particles which originate outside our solar 
system. They are believed to be ejected at high energies from nova and supernova 
explosions, solar flares from other solar systems, or quasars. The rays consists mainly of 
protons and Helium nuclei but may also consist of other heavier nuclei. Even though 
they amount to only a small portion of the total dose a spacecraft receives, galactic 
cosmic rays can lead to single event upset phenomena due to their high energies. The 
Van Allen radiation belts shield the Earth from galactic cosmic rays, and thus spacecraft 
in geosynchronous orbit are more susceptible to single event phenomena than those in 
lower altitude orbits. 
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2. 1 .3 Solar Particle Events 


Solar particle events are similar to galactic cosmic rays but are ejected from our 
Sun instead of from other solar systems’ suns. The Van Allen radiation belts also protect 
Earth from solar particle events. These events are fairly infrequent, on average a few per 
year. They consist mainly of protons which are ejected from the Sun during solar flares. 
A solar flare is created by a highly concentrated explosive release of high energy 
radiation. Events may last from a couple of hours to over a week. Typically, the effects 
last 2 or 3 days. Solar particle events can be a significant contributor to the total dose 
received by a spacecraft, in addition to creating single event phenomena. During a solar 
particle event, the flux levels of protons in the Van Allen radiation belts increase by a 
couple of orders of magnitude. The belts will retain high levels of particles in them until 
the particles dissipate; typically these high levels last from several hours to several days. 

2.1.4 Geomagnetic Substorms 

A geomagnetic substorm occurs when an interplanetary disturbance, a giant 
magnetic cloud containing hot gas from the Sun’s corona which is ejected from the Sun, 
strikes the Earth’s magnetic field. The magnetic cloud collides with, and envelopes, the 
Earth’s magnetic field, it compresses the magnetic field on the day side and stretches it 
out on the night side, thus causing the geomagnetic substorm. This phenomena is 
different from the solar particle event mentioned earlier. Substorms can last from several 
hours to a couple of days, depending on the severity of the storm. There exists a 
correlation between the geomagnetic activity index and the 1 1 year sunspot cycle. The 
geomagnetic activity index has a peak which occurs during the declining phase of the 
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sunspot cycle, and a secondary peak near solar maximum. Therefore, spacecraft can be 
affected by substorms during and for several years after an extreme solar maximum. 
Substorms consist of three distinct phases: growth, expansion, and recovery. The net 
result of a substorm is to increase the electron flux seen by spacecraft by 2 to 3 orders of 
magnitude. 5 

2.2 SURFACE CHARGING 

During the 1970’s, 1980’s, and continuing into the 1990’s, protection techniques 
have been developed which have basically taken care of the surface charging problem. 6 7 
Several of these have been published as design guidelines, “Design Guidelines for 
Assessing and Controlling Spacecraft Charging Effects” 8 , “Space Environmental Effects 
on Spacecraft: LEO Materials Selection Guide” 9 , “On Orbit Charging: Current TWR 
Design Requirements” 10 , and M1L-STD 1541 A". One of the powerful engineering tools 
used to analyze surface charging is NASA’s Charging Analyzer Program (NASCAP). 8 12 
It analyzes the surface charging of a three dimensional spacecraft surface as a function of 
time, for a given space environment, surface potentials, and material properties. By 
calculating and locating the regions of high surface voltage gradients, the areas where 
discharges are likely of occur are identified. Therefore, the materials or geometry of the 
spacecraft in these regions can be varied to minimize the likelihood of an electrostatic 
discharge. 

The general approach used by all of the above mentioned guidelines, is to select 
where possible the materials and conductive coating to be used on the surface of the 
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spacecraft, as to minimize differential charging of spacecraft parts. Differential charging 
is minimized when the voltage across the surface is as uniform as possible. If the above 
mentioned selection is not possible, the approach is then to use alternatives such as 
special filtering, cabling, or grounding. Once these have been selected, the charging 
codes are run to simulate the charging effects and to analyze whether the chosen 
configuration will work. However, since anomalies continue to occur on spacecraft, there 
must be another source of electrostatic discharges, such as deep dielectric charging. 
Therefore, an investigation of deep dielectric charging is required. 

2.3 DEEP DIELECTRIC CHARGING 

Vampola 11 showed that anomalies on spacecraft due to deep dielectric charging 
are linked to geomagnetic substorms, where the number of high energy electrons 
increases greatly. He was able to link the anomalies to deep dielectric charging by 
providing evidence that not all anomalies could be linked to surface charging. At 
geosynchronous orbit, deep dielectric charging is not affected by the spacecraft’s local 
time, while surface charging is known to occur between pre-midnight (22:00) and dawn 
(8:00), when the surface potentials become high enough for a discharge to occur. 
Therefore, between pre-noon (10:00) and dusk (20:00) no discharges should be observed. 
Since discharges do occur as shown in Figure 2.3, they must be due to deep dielectric 
charging instead of surface charging. 

Several spacecraft failures have been associated with electrostatic discharges 
resulting from deep dielectric charging, including the $300 million Telesat Canada 
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communications satellites Anik El and E2, Intelsat K, the European Space Agency (ESA) 
telecommunications spacecraft Olympus, and recently AT&T’s Telstar 401. Both Aniks 
failed a day apart in January 1 994, when they suffered a temporary loss of altitude control 
due to problems with their momentum wheels. It is believed that a discharge from deep 
dielectric charging disabled key circuitry in the momentum wheel control systems. Anik 
El was out of service for a few hours, but E2 was out for over 6 months and is now 
operating with a reduced lifetime on orbit, due to failure of the redundant momentum 
wheel control system. On the same day as the Anik El failure, Intelsat K’s momentum 
wheel control circuitry experienced an operation anomaly which caused it to lose attitude 
control. Control was reestablished when the backup circuitry was engaged. Analysis of 
the radiation environment at the time of the failures indicates that the amount of high 
energy electrons was greatly elevated. In August 1993, Olympus experienced a critical 
anomaly which lead to its eventual failure and out-of-service status. Telstar 401 
experienced a massive power failure in January 1997, rendering the spacecraft completely 
inoperable. This occurred after an intense geomagnetic substorm. An interesting point to 
note is that as Intelsat K, both Anik E’s, and Telstar 401 were all made by Lockheed 
Martin, these anomalies may be due to a similar design feature. 

Other spacecraft have experienced switchings or anomalies due to electrostatic 
discharges resulting from deep dielectric charging, including Voyager 1 as it passed by 
Jupiter in September 1977. There, it experienced a series of power on resets (POR). 
Each time a reset was executed, the onboard clock was offset. This offset caused a 
sequence of camera commands to be executed late and almost caused the loss of one of 
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the objectives of the photographic mission. The European communications spacecraft 
ECS-2 and ECS-4 experienced several automatic reconfiguration mode switchings in the 
early 1990’s. Anomalies where also reported on the US Air Force Defense Support 
Program (DSP) spacecraft, which affected the star-sensor shutter, and also on Meteosat-1 . 
These anomalies occurred on both spacecraft following periods of increased solar 
activity. Most anomalies are suspected to be due to discharges occurring in insulation 
around cables exposed to the outside of the spacecraft. 

There have also been electrostatic discharges on the Combined Release and 
Radiation Effects Spacecraft (CRRES). CRRES was sent up to measure outer-zone 
electrons and internal discharges during 1990 and 1991. 1416 The internal discharge 
experiment (IDM) exposed 16 samples of standard insulating materials with electrodes to 
the Earth’s radiation environment. The samples were composed of two different 
geometries, cables and printed circuit boards (PCB), and with various connections; 
grounding, floating, or semi-conducting elements. The samples were shielded with a thin 
aluminum foil, thus stopping all electrons below 150 keV. The orbit used was a high 
elliptical transfer orbit so that the spacecraft was exposed to all the magnetosphere 
environment. Over the 14 months of the flight, approximately 4300 spontaneous 
discharges were recorded from the samples, thus proving that deep dielectric charging did 
exist and that it is of concern. The results are of the form of the flux above a certain 
energy level on the days that discharges were recorded. It should also be noted that while 
some spacecraft have been affected by deep dielectric charging, other spacecraft of 
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similar design and in operation at the same time have not been affected. This 
demonstrates that more research into the issue of deep dielectric charging is required. 

A general overview of deep dielectric charging was done by Garrett and 
Whittlesey 17 , and Soubeyran 18 , who also gives some basic design guidelines to follow to 
minimize deep dielectric charging. Some early modeling of deep dielectric charging was 
done by Berkley 19 and Frederickson 20 . Berkley analyzed the charging of thin polymer 
films, using multiple electron energy ranges. The model can incorporate both open and 
short circuit boundary conditions, and can include non-linear transport terms, as well as 
time-varying radiation induced conductivity. The electron-irradiation parameters, the 
radiation induced conductivity and the electron deposition profile, are computed using a 
single-scattering Monte Carlo model of the polymer. The time-dependent charge- 
transport equations are solved by direct discretization of the partial differential equations. 
With his code, Berkley was able to study the effects of a time-varying electron beam, 
conductivity, and the influence of various boundary conditions, using both linear and 
non-linear transport. The results showed that charging reaches a steady-state quickly, and 
that the model has excellent qualitative and in many cases quantitative agreement with 
the experiments he performed using an electron beam to irradiate samples. 

Frederickson used basically the same model as Berkley, but improved the 
calculations to include the dose effects of secondary X-rays and bremsstrahlung radiation, 
and to include the effect of the electric field on the trajectories of the incident penetrating 
electrons. A Monte Carlo simulation is used for the electron trajectories of the mono- 
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energetic electron beam, and a computerized iteration technique is used for the motion of 
the space charge and the induced electric fields. Two different boundary conditions are 
used, both the front and rear surfaces are held at zero potential (grounded), and the front 
surface allowed to float free. The first boundary condition represents the spacecraft 
surface in full sun which allows solar photoemission to keep the surface at zero potential. 
The second represents the spacecraft surface in full shade where without solar 
photoemission the potential can vary. It was found that the bremsstrahlung and x-rays 
produced very little effect, and that the space charge electric fields do not significantly 
alter the trajectories of the incident electrons. Both Berkley’s and Frederickson’s models 
work well, but do not simulate the space environment, as they work only for an electron 
beam of mono-energetic electrons. 

Recently, Soubeyran 21 of Matra Marconi Space developed ESA-DDC" 2 ’ 
(European Space Agency Deep Dielectric Charging), a numerical tool to analyze deep 
dielectric charging in the space environment, for the European Space Agency. The 
electron and proton transport is calculated using Monte Carlo transport codes to track the 
path followed by the electrons and protons. This provides the deposited dose and charge 
within the material for a given energy distribution or mono-energetic particles. Ampere’s 
equation with Ohm’s law are transformed into equivalent electrical circuits and sent to a 
circuit analyzer. The circuit analyzer calculates the induced conductivity, current, 
electrical potential, and electric field. The code is restricted to 1-D geometries with the 
upper surface exposed to the space environment and the lower surface grounded to the 
structure potential. The material between the surfaces can be made up of multiple layers 
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of materials, both conducting and insulating. Subsequently, the code has been modified 
for the analysis of 2-D geometries. 

Another method, which was proposed by Whittlesey 24 , is to manufacture the 
spacecraft and before it is launched, to test it to see if electrostatic discharges occur. 
Testing is performed using the MIL-STD 1541 sparking device, which simulates the 
effects of space electrostatic discharges. Thus if discharges occur modification of the 
spacecraft is required after the spacecraft is designed and assembled. This technique can 
be effectively used to double-check spacecraft which have been designed to have minimal 
electrostatic discharges, but is grossly inefficient for finding and fixing problems. 

Most previous work on deep dielectric charging has dealt with the charging of 
spacecraft components such as wire insulation and printed circuit boards, where the 
electrostatic discharges can lead to anomalies in spacecraft operation. However, with the 
increasing use of composite materials, which can be dielectrics, deep dielectric charging 
is becoming a problem with the structure of spacecraft. 

2.4 COMPOSITES 

Composites are starting to replace metals, such as aluminum, as the structure of 
spacecraft, due primarily to their higher stiffness-to-weight ratios. Composites are made 
up of multiple layers or plies, which are stacked at various angles to get the desired 
material properties. These plies are in turn made up of fibers and matrix material which 
surrounds the fiber. The dominant fibers used are carbon (sometimes referred to as 


45 



graphite), glass, and Kevlar. The dominant matrix materials for space applications are 
epoxies, cyanates, and poly-ether-ether-keytones (PEEK). 

By varying the angles of the plies and the stacking sequence of the plies (refer to 
Figure 2.4), one can tailor the properties of the laminate or composite structure. The 
reason changing the ply angle has an effect on the material properties is that the plies are 
anisotropic. The material properties in the longitudinal direction can be very different 
from the properties in the transverse and through-thickness directions (refer to Figure 
2.5). The ply angle is defined as the angle between the geometric coordinate system and 
the ply coordinate system. The range of angles is from +90° to -90°. The geometric 
coordinate system is arbitrarily assigned to a structural direction, for example the length 
of a solar panel array, and the ply coordinate system is aligned with the fiber direction, as 
shown in Figures 2.4 and 2.5. Therefore a 0° ply will have its longitudinal properties 
aligned with the principal direction of the structure. The laminate properties are based on 
the ply properties and the ply angles, and can be calculated using Classical Laminated 
Plate Theory (CLPT). For more information refer to Jones 25 and Tsai 26 . 

In modeling composites in general, and specifically for the charging problem, the 
homogenized equivalent properties for the composite ply are used instead of the 
individual properties of the fibers and matrix. This is an acceptable simplification for a 
first order solution of the problem, as the thickness of the material is much greater than 
the diameter of the fibers for most material systems of interest. To determine the ply 
properties, one has to combine the fiber properties with the matrix properties. The ply 
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Figure 2.4: Typical composite in the geometric coordinate system 
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properties are dependent on the volume fraction of fiber used in the composite material, 
the matrix properties, and the fiber properties. There are many models which are used to 
determine the ply properties, all of which have modeling limitations. Typically the 
simplest models are those which are used to determine ply properties in the longitudinal 
direction. The most basic models are the “rule of mixtures” (ROM) and the “inverse rule 
of mixtures” (IROM), which result in relations similar to the equations for resistors or 
springs in series and in parallel. These models are independent of the details of the fiber 
and matrix geometry. More complicated models are required for the more difficult to 
model properties, such as the transverse and through-thickness directions properties. 
These models are can be based on experimental or analytical work, and therefore different 
researches can develop greatly different models for the same property. More specific 
details are given in section 4.5. 

However, the details of the fiber and matrix can matter, especially in cases where 
the properties of the fiber and matrix differ radically. This is the case with electrical 
properties. Geometry can be important, as in the case of a low resistance percolation 
path, as shown in Figure 2.6. Another case is the insulative epoxy rich surface layer on 
composite laminates, as shown in Figure 2.7. These details will be examined on a case 
by case basis in section 4.4. 

The materials used in this research are carbon fiber/epoxy, glass fiber/epoxy, and 
Kapton. The epoxy, glass fibers, and Kapton are insulators, while the carbon fibers are 
conductors. The electrical properties of these materials are tabulated in Table 2. 1. 27 ' 36 
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Figure 2.6: Typical composite with a percolation path 
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Table 2.1: 


Material properties of materials used in this research 


Material 

Density 

(kg/m 3 ) 

Conductivity 

(1/Q-m) 

Dielectric 

Constant 

Dielectric 
Strength (MV/m) 

Carbon Fibers 

1384- 2200 

2.0 x 10 4 - 10 6 


0.0032 - 0.0044 

Epoxy 

1052 - 2187 

10‘ 3 ~ 10' 8 

2.78 - 5.2 


Glass/Epoxy 

1550 - 2076 

< io '° 

4.2 - 5.68 

17.7 -21.7 

Carbon/Epoxy Long. 

1577 - 1700 

374 - 47,600 



Carbon/Epoxy Trans. 

1577 - 1700 

1.5 - 2000 



Carbon/Epoxy Thick. 

1577 - 1700 

0.1 - 106 



Kapton 

1420 - 1670 

7 x 10" 16 ~ 10' 15 

2.7 - 3.5 

154 - 303 


52 































CHAPTER 3 


PROBLEM STATEMENT AND APPROACH 

3.1 PROBLEM STATEMENT 

In this work, we will develop a model of the deep dielectric charging of 
composites such that, given the orbit data (apogee height, perigee height, inclination, and 
solar cycle), the composite electrical properties (conductivity, dielectric constant, and 
dielectric strength), the density, and the thickness, we can calculate the distribution of 
charge density, electric field, and voltage with respect to position through the thickness of 
the composite. We will make suggestions to reduce the probability of electrostatic 
discharges from occurring due to deep dielectric charging. We will also experimentally 
measure the continuum electrical properties of carbon fiber / epoxy composites. Finally, 
a composite material system will be developed with conductivity properties that can be 
tailored without affecting its mechanical properties. 

3.2 APPROACH 

Previous literature on spacecraft charging is used to develop an understanding of 
the space environment and the charging phenomena. This understanding is used to 
develop the analytical models. The analysis has three fundamental goals. The analysis is 
used to identify key parameters associated with deep dielectric charging and the resulting 
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electrostatic discharges, as well as sensitivity of the analysis to these parameters. 
Secondly, the analysis provides insight into the interaction of the fiber and matrix on a 
micro-scale. This is required since in the deep dielectric charging model the composite is 
treated as a homogeneous material, when in fact it is made up of fibers surrounded by 
matrix. The analysis will also be used to examine specific orbits as case studies, where 
spacecraft are suspected to have suffered problems due to deep dielectric charging. 

The purpose of the experimental program is gain an insight into the parameters 
which control the through-thickness electrical properties of composites. The critical 
properties to be investigated are determined from the analytical program; they include the 
conductivity and the dielectric strength. These properties are initially assumed to be 
continuum properties. The second purpose on the experimental program is to investigate 
a conductivity tailorable composite. This is done to aid at reducing the chance of 
electrostatic discharges occurring, since increasing the conductivity decreases the electric 
field. Therefore, by increasing the conductivity, the likelihood of an electrostatic 
discharge is reduced. This material system is to have tailorable electrical properties 
without changing the mechanical properties significantly. 

3.3 ANALYTICAL TASKS 

The analysis is made up of two parts, the deep dielectric charging analytical 
model and the modeling of the fiber/matrix micro-scale effects. The deep dielectric 
charging model comprises several parts. First comes the acquisition of the environmental 
data from the Environmental Workbench software. In the actual charging model, the 
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composite is treated as a homogenous material and the fiber/matrix micro-scale effects 
are studied separately. The model includes the effects of the space environment on the 
conductivity of the composite, including radiation induced conductivity and high field 
conductivity. The model is used to perform sensitivity studies on various parameters to 
determine which ones affect the likelihood of discharges occurring. These results gave 
direction to the experimental program, as to which properties should be examined. The 
model is also used to evaluate the likelihood of discharges occurring on various case 
study orbits. 

The composite is modeled as a homogenous material in the charging model, when 
in actuality it can be highly anisotropic and locally inhomogeneous. It is made up of 
fibers surrounded by matrix; therefore various effects of the fiber/matrix micro-scale are 
modeled. These include the effects of a conductive fiber surrounded by an insulating 
matrix, unlikely fiber arrangements, and the pure epoxy surface layer inherent on 
composites. 

3.4 EXPERIMENTAL TASKS 

There are two major experimental tasks. The first task is to get good 
measurements of through-thickness conductivity and dielectric strength of carbon fiber / 
epoxy composites. The effects of laminate thickness, surface area, and ply stacking 
sequence are investigated. The second task is create a conductivity tailorable composite. 
This task will involve adding conductive carbon black powder to a glass fiber / epoxy 
composite to control the conductivity. The material properties, both electrical and 
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mechanical, are measured and plotted against percentage of carbon black to create 
preliminary design charts. 
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CHAPTER 4 


THEORY 


This chapter presents the theory used in the development of the composite deep 
dielectric charging analyzer. The geometry being analyzed is presented first, followed by 
the basic electrostatic equations used. The development of the deep dielectric charging 
model is then presented. Additional models used to model micromechanical details and 
composite electrical micromechanics are described. Finally, the equations used to reduce 
the experimental data are presented. 

4.1 MODEL GEOMETRY 

The geometry of the composite is a 3-D plate with longitudinal, transverse, and 
through-thickness material properties, as shown in Figures 2.4 and 2.5. However, for the 
deep dielectric charging model the composite is modeled as a 1 -D plate in the through- 
thickness direction with homogenized material properties, as shown in Figure 4.1. The 
incoming radiation from the space environment is assumed to be perpendicular to the 
surface, and uniform across the surface. The front surface of the composite is the surface 
that is exposed to the space environment, and the back surface is the surface that faces the 
inside of the spacecraft. The depth (x) is defined as starting at the front surface and 
ending at the back surface. 
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Figure 4. 1 : Geometry used in the CoDDC A model 
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4.2 BASIC ELECTROSTATIC EQUATIONS 

The equations used in modeling deep dielectric charging are derived in this 
section. The continuity equation 37 relates the current density to the charge density. The 
charge density includes both the incoming charge density rate from the space 
environment and the charge density in the material. 

dp 

v ' J= ~aT +p - (l) 

where J is the current density vector (A/m 2 ), t is the time (s), p is the charge density 
(#/m 3 ), and p in is the incoming charge density rate (#/m 3 -s) from the space environment. 
The current density can be represented by: 

J = ct E (2) 

where a is the conductivity (1/Q-m) of the material. It should be noted that the 
conductivity can be a tensor for composites, and E is the electric field vector (V/m). 
Rewriting equation 1 for the 1-D case, and substituting equation 2 for J: 

dp . 3<yE 

dt’ = Pin__ ^T (3) 

where x is the depth (m) into the material. 

Gauss’ Law in differential form 37 relates the electric field to the charge density: 

VE = ip (4) 

where e is the permittivity (C 2 /N-m 2 ) of the material, defined as: 

e = Ke 0 (5) 
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where k is the dielectric constant of the material, and e„ is the permittivity of free space 
(C 2 /N-m 2 ). Rewriting equation 4 for the 1-D case: 


3E 1 
dx ~ e P 


( 6 ) 


Poisson’s equation 37 relates the voltage to the charge density: 

1 

V 2 V = — p 
8 

where V is the voltage (V). Rewriting equation 7 for the 1-D case: 

d 2 V 1 


(7) 


( 8 ) 


It can be seen that the voltage and electric field are related by 

V V = -E 

Rewriting equation 9 for the 1 -D case: 


In the space environment, the conductivity of a material is not constant, it is 
affected by the incoming radiation and the electric fields produced within the material. 
Therefore the conductivity can be expressed as: 

<7 = a 0 k E +a R +cr D (11) 

where a 0 is the dark conductivity (1/Q-m) of the material, k E is the conductivity scaling 
factor for the effect of high electric fields on conductivity, a R is the radiation induced 


(9) 


( 10 ) 
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conductivity (1/Q-m), and ct d is the delayed conductivity (1/Q-m). The conductivity in 
materials occurs because atoms have electrons which are free to move through the 
material. With no electric field applied, the electrons move in random directions, while 
when a field is applied to the material, the electrons drift in the opposite direction of the 
field. 


4.2.1 Conductivity In High Electric Fields 

Additional conduction is caused by the generation of charge carriers (electrons) by 
the ionization of neutral centers (atoms) in the bulk of the material. This ionization is due 
to a thermal activation process which may be modified by an applied electric field: at low 
fields there is minimal effects. The conductivity of insulating materials exposed to high 
electric fields is modeled based on a ‘classic dielectric’ approach, 38 as follows: 

°E= CT o k E (12) 

where k E is defined as: 


k E = 


2 + cos h(p F VE /2kT) 


f 2 kT (eE5^ 
sinhl 


eE5 


2kT; 


(13) 


where k is Boltzmann’s constant (1.380658 x 10' 23 J/K), T is temperature (K), e is the 
charge of an electron (C), 8 is the jump distance (m), and P F is the Frenkel parameter 39 

( -\/C - Vn -m) defined as: 


Pe = 


( 





(14) 


In k E the first term in brackets is the field dependence of carrier concentration and the 
second term is the field dependence of carrier mobility. It has been shown 38 that because 
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of small jump distance, only a few Angstrom units, the rise in the carrier mobility may be 
neglected up to fields of approximately 100 MV/m. Since this value is above the 
dielectric strength of most materials, the second term will be set to one in the charging 
model. It should be noted that there is extremely good agreement between this model and 
measure values of high field conductivity. 

4.2.2 Radiation Induced Conductivity 

When insulators are exposed to radiation, it is known that the conductivity of the 
material increases. 40 ' 42 This occurs because the number of additional conducting ion- 
electron pairs produced in the material is proportional to the energy deposited by incident 
radiation. This radiation induced conductivity is defined as: 

a R =k R D A (15) 

where k R is the coefficient of radiation induced conductivity (s/Q-m-rad), D is the 
radiation dose rate (rad/s), and A is the radiation induced conductivity exponent which is 
a material dependent parameter. Both k R and A have been empirically determined for 
many materials by many investigators. k R is relatively small for most organic insulators, 
and can be 2 or more orders of magnitude greater for inorganic insulators. It should also 
be noted that k R for a given material can vary by as much as 2 orders of magnitude 
between different sources. This could be due to older experimental techniques used in 
some of the older sources and the use of samples manufactured at different times with 
different manufacturing procedures. The radiation induced conductivity exponent (A) 
usually lies between 0.5 and 1 .0; with most organic substances it lies around 1 .0. 


62 



4.2.3 Delayed Conductivity 

The delayed conductivity is defined as: 

a D =k D Dt (16) 

where k D is the coefficient of delayed conductivity (1/rad-Q-m). Frederickson 20 found 
that delayed conductivity is negligibly small, and will therefore be assumed to equal zero 
and will not be included. 


4.2.4 Conductivity Summary 

The expression for the conductivity that will be used in the charging model is as 
follows: 


a = a. 


2 + cosh(p F VE/2kT)^ 


+ k R D A 


(17) 


4.3 DEEP DIELECTRIC CHARGING MODEL 

The basic problem is that high energy particles from the space environment 
penetrate the surface of the composite material, and start losing energy until they stop 
somewhere within the material. These stopped particles induce an electric field within 
the material. This electric field causes the particles to move, thus causing a current in the 
material. This current in turn influences the electric field. The electric field continues to 
grow until equilibrium is reached between the incoming particles from the space 
environment ( p jn ) and the particles leaving due to the current (J). This is also known as 
the steady state. Another possibility is that the electric field exceeds the dielectric 
strength of the material before equilibrium is reached, at which point a breakdown and 
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subsequent electrostatic discharge occurs. A model of this process will be described. A 
Fortran code, CoDDCA (Composite Deep Dielectric Charging Analyzer), based on this 
model, was written. It calculates the electric field, voltage, and charge density as 
functions of depth into material and time, until one of the two above mentioned 
possibilities occur. Appendix A contains the user manual for CoDDCA. Magnetic field 
effects are assumed to be negligible because of the low current density (< 10 nA/cm 2 ) and 
low velocity of charge carriers in insulators (< 105 m/s « c). 18 

4.3.1 Environment 

The charged particle environment is acquired from the Environment Workbench 
(EWB) software, written by the S-Cubed division of Maxwell Laboratories, Inc. 41 The 
program runs the electron and proton models AE8 and AP8, and gives empirical 
omnidirectional flux data at different energy levels for different orbits and solar cycle 
conditions (maximum or minimum). The output from the Environmental Workbench is 
unidirectional penetrating flux, EWB converted it internally. The E and P in the model 
names AE8 and AP8 refer to electron and proton, and 8 is the version number of the 
models. These models of the near Earth radiation regime were created by the National 
Space Science Data Center, from data compiled from many different spacecraft. AE8 
gives the fluxes of electrons in the energy range 50 keV to 7 MeV, and AP8 gives the 
fluxes of protons in the energy range 50 keV to 500 MeV. The time dependent variations 
of radiation fluxes such as those due to geomagnetic substorms or short term solar 
modulations such as solar particle events are not included. 
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The depth the particles penetrate into the material is calculated using empirical 
information on the stopping power of silicon. The stopping power is given as the linear 
energy transfer (LET) which is a function of energy, as shown in Figure 4.2. The 
stopping power in silicon can be used, even thought the materials used here are 
composites, because the linear energy transfer scales with density. This allows it to be 
used with any material as long as its density is known. Using the particle flux versus 
energy from EWB and the LET curves, the particle flux versus depth of penetration into 
the material can be calculated. EWB produces a distribution of particles as a function of 
energy, as shown in Figures 4.3. The LET curve provides the depth a particle penetrates 
into a material as a function of energy. These two sets of data are combined graphically 
to produce a distribution of particle flux as a function of depth into the material, as shown 
in Figures 4.4. This is accomplished by determining equations for the electron and 
proton LET, as shown in Figures 4.5, which can be used to convert the EWB energy data 
into penetrating depth data, using the following equation: 


x = 


LET p mech 


(18) 


where Eg is the particle energy (MeV), and p mech is the density (kg/m 3 ) of the material. 
Figures 4.6 show the relationship between the electron and proton energy and the depth 
the particles penetrate. 


Finally, to determine the incoming particle charge density rate, as shown in Figure 
4.7, the incoming particles per area (particle flux) must be converted to incoming 
particles per volume (charge density rate). This is easily accomplished since the flux data 
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Figure 4.2: Stopping powers in silicon for various particles 44 
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Proton Flux (10 9 #/m 2 -s) Electron Flux (10 9 #/m 2 -s) 




Figure 4.3: Electron and proton fluxes as a function of energy 
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Figure 4.4: Electron and proton fluxes as a function of depth 
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Figure 4.5: Electron and proton LET as a function of energy 
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Figure 4.7: Electron and proton charge density rates as a function of depth 
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is grouped in discrete bins, each bin is at a different depth in the material. Therefore, by 
assuming that the flux of particles that penetrate to a certain depth are evenly distributed 
from midway back to the previous depth to midway forward to the next depth (Ax), as 
shown in Figure 4.8. Thus mathematically the charge density rate of a certain depth bin 
is equal to the flux associated with that depth bin divided by Ax. The difference between 
the electron charge density rate and the proton charge density rate gives the net incoming 
charge density rate, p in , as a function of penetrating depth, as shown in Figure 4.9. 

This is a simpler approach than what Soubeyran 21 did with a Monte Carlo 
deposition model. However, a simpler model is acceptable for our work since the model 
is primarily being used to find order-of-magnitude estimates of the problem and to give 
direction to the experimental portion of this research. The added accuracy of the Monte 
Carlo deposition would not change the general conclusions of our model or the choice of 
parameters to be examined experimentally. The most likely materials to be susceptible to 
deep dielectric charging will still be identified as well as which parameters affect the 
likelihood of an electrostatic discharge. The main difference between the Monte Carlo 
model and our model is that in our model all particles of a certain energy penetrate to a 
certain depth, while with the Monte Carlo model there is a distribution as to where the 
particles of a certain energy penetrate to. This would have the effect of spreading out the 
‘bins’ used in our model and smoothing out the distributions of particles with depth 
somewhat. 


72 



Incoming charge 
density (#/m 3 -s) 


Flux (#/m 2 -s) 


Figure 4.8: 



Representation of how the particle per area data is converted into particle 
per volume data 




Figure 4.9: Incoming charge density rate as a function of depth 
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For certain calculations, such as the radiation induced conductivity, the dose rate 
is required. The dose rate is defined as: 45 

6 _ f <*E„)dE.(dE,/dx)QF(E.) ()9) 

E « P mech 

where <D(E g ) is the charged particle flux (#/m 2 -s) with kinetic energy between E g and E g + 
dE g , dE g /dx is the stopping power (MeV/m) of the material, and QF(E g ) is the quality 
factor for particles of energy E g . The quality factor is defined as 1 for electrons and 10 for 
protons with energy less than 14 MeV. Equation 19 can also be written for numerical 
data as: 




O, (dE g / dx) QFj 

P mech 


( 20 ) 


where i represents the parameter at energy i. This equation can be simplified by 
observing that: 


dE„ /dx 
— k = LET 

P mech 


( 21 ) 


therefore, equation 20 becomes: 

D = Z,°i LET , Q F i (22) 


4.3.2 Closed Form Solution 

A closed form solution to the deep dielectric charging problem can be attempted if 
it is assumed the conductivity of the material is constant in space (the x direction). This 
is the case if high field conductivity is neglected. Therefore equation 3 becomes: 




( 23 ) 
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Substituting equation 6 into equation 23: 


dp ct 

^T = P.n -TP 
dt 8 


(24) 


this equation can be solved for the charge density. Since p is a function of x and t, and 
p jn is only a function of x. Re-writing, equation 24 becomes: 


dp(x,t) ct 

— — — + -p(x,t) = p in (x) 
dt 8 


(25) 


Assuming p(x,t) is separable: 

p(x,t) = X(x) T(t) 

Substituting equation 26 into equation 25, and re-arranging: 


d T a p in 

— + — T - — 
dt 8 X 


(26) 


(27) 


This differential equation is general first order linear, and can be solved as: 


p(x,t) = -p,„(x) + X(x)e E ' 
a 


X(x) can be solved using the initial condition that at t = 0, p(x,0) - 0: 


X(x) = --p in 


Therefore, equation 28 becomes: 


(28) 


(29) 


P( X >0 = “ Pin(x) _e £ ' 


( 30 ) 


Thus the steady state value of the charge density, p(x,oo) is: 
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( 31 ) 


p(x,oo) = -p in (x) 

CT 

The time to steady state, defined as when p reaches 99% of its actual steady-state value, 
is defined as follows: 


t = - — ln(0.01) 
a 


(32) 


The electric field is calculated using equation 6: 


E = ijpd, 


Therefore, by integrating equation 30, the electrical field is: 


E = 


1 


a 

J - e E ' Jpin(x)dx + C, 


(33) 


(34) 


This equation is integrated numerically using the trapezoidal rule. The voltage is 
calculated using equation 10: 


V = - J Edx 


(35) 


Therefore, by integrating equation 34, the voltage is: 


V = -- 

CT 


1_e ^ , ll PinOO^ +C,X + C 2 


(36) 


This equation is also integrated numerically using the trapezoidal rule. 


It should be noted that only the incoming charge density rate needs to be 
integrated, as the rest of the equation is not a function of x. This allows the integrals to 
be evaluated only once for the electric field and voltage at any desired time. The 
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boundary conditions C, and C 2 can be calculated based on either one electric field and 
one voltage boundary condition, or two voltage boundary conditions. 


4.3.3 Finite Difference Solution 

If the high field conductivity is not ignored, the material conductivity is not 
constant in space, and cannot be removed from the differential in equation 3. Therefore, 
by using the product rule for differentials, equation 3 becomes: 



SE do 
dx dx 


(37) 


Substituting equation 6 into equation 37: 


dp o do 

¥ = |i "'-7 p - E ^ 


This equation is solved numerically with a central difference in space, forward difference 
in time routine. 


The electric field can be calculated using equation 33: 

E = ^Jpdx + C, (39) 

This equation is integrated numerically using the trapezoidal rule. The voltage can be 
calculated using equation 35: 

V = -^ jjpdx + C, x + C 2 (40) 

This equation is also integrated numerically using the trapezoidal rule. It should be noted 
that since the charge density changes, these integrals need to be evaluated for every time 
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step. The boundary conditions C, and C 2 can be calculated based on either one electric 
field and one voltage boundary condition, or two voltage boundary conditions. 

Equation 39, equation 40, and equation 38 are solved sequentially, in a time 
stepping routine. This continues until either steady state is reached or the electric field 
exceeds the dielectric strength of the material and an electrostatic discharge occurs. 

4.3.4 The Analytical Model 

To summarize, a tool now exists which can determine the charge density, electric 
field, and voltage as functions of depth into the composite material and time, in a 
specified orbital environment, and thus predict if an electrostatic discharge will occur. 
The program can run either model, and can be run for any orbit by specifying the orbit 
parameters (the apogee height, perigee height, inclination, and solar cycle condition) in 
EWB. Substorms and solar particle events can also be user specified, as well as the 
material parameters. The material is characterized electrically by its conductivity, 
dielectric constant, dielectric strength, density, and boundary conditions. The structure is 
characterized by its thickness, temperature, and boundary conditions. There are also 
options to include radiation induced conductivity and high field conductivity. 

4.4 MICROMECHANICAL DETAILS 

This section deals with the modeling of the interactions between the fibers and 
matrix of the composite material and their effects on the electric field. The first model 
examines the electric field in the matrix as a function of the radial distance from a 
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conductive fiber, in which the electric field is zero. The second model examines the 
insulative pure epoxy surface layer of carbon fiber / epoxy composites. 

4.4,1 Micro-scale Model 

At the micro-scale, a concern is that higher electric fields could be possible due to 
geometric details of the fibers imbedded in the matrix. This model examines the steady- 
state electric field in an insulative matrix which passes around a conductive fiber, as 
shown in Figure 4.10. Here r is the radial distance (m) from the center of the fiber, 0 is 
the angular distance counter-clockwise from the x axis, R f is the radius of the fiber (m), 
and E„ is the steady-state electric field (V/m) calculated using the composite deep 
dielectric charging model. The electric field is assumed to be constant far away from the 
fiber, since the scale of the problem is much smaller than the scale used in the composite 
deep dielectric charging model. 


To solve the problem, Poisson’s equation (equation 7) needs to be expressed in 


radial coordinates: 


J__3_ 
r 2 3r 


1 


3V 


dx ) r 2 sin(9) 30 


' / v avi 

sin(0) -rrr I + 


3V 1 


V 


30 J r 2 sin 2 (0) 3<|> e 


(41) 


assuming p = 0 for simplicity, equation 41 becomes Laplace’s equation: 


\_d_ 

X 2 dx 


f 


V 


3V 
3r , 


+ 


r 2 sin(0) 30 


3 ( • ^ 5V 1 

sm(e) — + 


3V 


30 J r 2 sin 2 (0) 3<f> 


= 0 


(42) 


assuming V is independent of <j), and separable in r and 0, the solution to equation 42 is: 
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X 



Figure 4. 10: Schematic of the fiber-matrix micro-scale model 
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( 43 ) 


V(r,0) = |;fA i r i +4-1^(008(0)) 

i=0 V r y 

where A ; and B, are arbitrary constants associated with the solution of 2 nd order 
differential equations, and P,(x) is defined by the Rodrigues formula as: 

Pi(*) = ^f;r) f* 2 - 1 ) 1 (44) 

2 i! ydxy 

The first boundary conditions required to solve for the constants is the voltage is zero at 
the fiber-matrix interface since the fiber is a conductor. The second boundary condition 
is that far from the fiber the voltage approaches the far field voltage created by the far 
field electric field (E 0 ). The far field voltage is equal to the electric field multiplied by the 
distance from the center of the fiber. The boundary conditions in circular coordinates are: 

V = 0 @ r = R f 

, . (47) 

V => -E 0 r cos(0) @ r » R f 

Applying the first boundary condition, and observing that the cosine term cannot 
be always equal to zero, the second constant is determined in terms of the first constant: 

Bj = -Aj R f 2i+I (48) 

Therefore, equation 43 becomes: 

00 ( 2i + l \ 

V(r,0) = lA, r' -- P^co^)) (49) 

i=o V r / 
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Applying the second boundary condition, and observing that at r » R f the second 
term in the brackets of equation 49 is negligible, the first constant can be solved for by 
equating like terms. The first constant is solved as: 

Aj = -E 0 (50) 

Therefore the voltage solution is: 


V(r,9) = -E„ 


r - 


R f 




r 2 / 


cos! 


(e) 


(51) 


To calculate the electric field equation 9 is used in radial coordinates: 

E(r,0) = - 


3v A i av - 

"a7 r + 7 ae 6 


Therefore the electric field is: 
E(r.0)=E„ 

The magnitude of the electric field is: 


cos(9) 

f r 3 A 

1 + 2 3 r - sin(0) 

i- R : 

0 


\ X J 

l x ) 



/ 

f R, 3 R f b ) 


f R f 3 R, 6 ^ 

cos 2 (0) 

1 + 4 +4 

+ sin 2 (©) 

1-2 + 

\ 

l r r J 


l r r )) 


(52) 


(53) 


(54) 


4.4.2 Epoxy Surface Laver 

Composite laminates are made up of a fairly even distribution of fibers and 
matrix, however at the top and bottom of the laminate there exists a surface layer 
comprised of mainly epoxy, as shown in Figure 4.11. This surface layer is highly 
insulative, and can lead to problems in conductive composites such as carbon fiber / 
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Figure 4.11: Schematic of the pure epoxy surface layer 
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epoxy. Problems occur because while the bulk composite is conductive, right at the 
surface it is insulative, resulting in a large (many orders of magnitude, as shown in Table 
2.1) jump in the conductivity. Therefore, the time to steady state also jumps by many 
orders of magnitude. 


To calculate the electric field in the surface layer, the layer is modeled as a 
resistor and capacitor in parallel, as shown in Figure 4.12. Here C SL is the capacitance (F) 
of the surface layer, R SL is the resistance (Q) of the surface layer, h SL is the thickness (m) 
of the surface layer, V is the voltage (V) across the surface layer, i jn is the incoming 
current (A) from the space environment, i c is the current (A) flowing through the 
capacitor, and i R is the current (A) flowing through the resistor. The capacitance of the 
surface layer can be related to its dielectric constant by: 


r = 

SL h 


s A 


(55) 


-SL 


where A is the cross sectional area (m 2 ). The resistance can be related to the conductivity 


of the surface layer by: 


R - 

Rsl ‘ 0 a 


(56) 


The electric field across the surface layer is defined as: 


E = 


(57) 


i SL 
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Using Kirchhoff s conservation of charge rule, the current in the model can be 
expressed as: 

*in = *c + *r 

the incoming current can be related to the incoming charge density from the space 
environment at the front surface by: 


iin=Pi„ h SLA (59) 

Using the equations for the current flowing through a capacitor, the current flowing 
through a resistor, and equation 59, equation 58 becomes: 

dV 1 


Pin^st-A - ^"SL + 


(60) 


SL 


Substituting equations 55 and 56 into equation 60: 

sA dV ctA 

Pi " hsLA = ^"dT + hT 

Re-arranging and simplifying, equation 61 becomes: 


(61) 


dV a 

dt + 8 V_Pi " 8 


SL 


(62) 


This differential equation is general first order linear, and can be solved as: 


V = P in 


‘SL 


a 

• + C, e E ' 


(63) 


where C 3 can be solved for with the initial condition that the voltage is zero at time equal 
zero. Therefore equation 63 becomes: 
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( 64 ) 



Substituting equation 64 into equation 57, the electric field across the surface layer is: 

E = p, n ^( 1 - e ''j (65) 

Therefore, the steady-state electric field is: 

E = P„^ (66) 

The time to steady state, defined once again as when E reaches 99% of its actual steady- 
state value, is defined as: 

t = -— ln(0.01) (67) 

a 

It should be noted that this is the same equation as for the steady-state time for the 
composite deep dielectric charging model closed form solution (equation 32). 

4.5 COMPOSITE MICROMECHANICS 

4.5.1 Continuous Fiber Composites 

This section will give an overview of the micromechanics of composite electrical 
properties. As referred to earlier, micromechanics is the determination of ply properties 
based on the constituent properties of the fiber and matrix. Micromechanics is used to 
homogenize material properties; this may or may not be a good idea depending on the 
usage of the material property. A general overview of the topic is given by Kaddour et 
al. 34 An important definition in micromechanics is the fiber volume fraction, defined as: 
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( 68 ) 


V Ulf 

v, = Vol ( +Vol„ 

where V f is the fiber volume fraction, Vol f is the volume of fibers (m 3 ) in the composite, 
and Vol m is the volume (m 3 ) of matrix in the composite. 

The composite conductivity in the longitudinal direction (ct l ) is based on a simple 
rule of mixtures as: 46 

c L =«,V, + CT„(l-V r ) (69) 

where cr f is the conductivity (1/Q-m) of the fiber, and CT m is the conductivity (1/Q-m) of 
the matrix. 


The composite conductivity in the transverse direction (a x ) is more complicated. 
One model expresses it as: 46 




a f -a 
= 2 V f — 


a T +a f 


(70) 


The equations given above for the longitudinal and transverse conductivity were derived 
analytically using a continuum approach. The self consistent model assumes that the path 


of each fiber can be divided into a series of short straight segments. Each segment is 
treated as a long thin rod for the purpose of calculating the average electric field within a 
fiber. The presence of the other fibers is taken into account by assuming the segment is 
embedded directly in a matrix having the bulk properties of the sample. It should be 
noted that equations 69 and 70 can also be used to calculated the dielectric constant of the 
lamina because of the identical form of the conduction laws controlling both properties. 
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Therefore the conductivity variables can be replaced with dielectric constant variables in 
the given equations. 


To calculate the in-plane electrical properties (dielectric constant and 
conductivity) of a plain weave fabric composite, refer to Chou and Ning 47 . For an angled 
ply, the conductivity in the principal directions are expressed as: 48 


a, = a L cos 2 (0) + a T sin 2 (0) 

a 2 = a L sin 2 (0) + a T cos 2 (0) 
where 0 is the ply angle (refer to Figure 2.4). 


(71) 

(72) 


Volpe 49 suggested an empirical expression for estimating the conductivity of 
laminates made of the family of [0/±45/90] plies, in the fiber directions as: 


= 


K + {n„) 


n t 


(73) 


where N e is the number of plies in the 0 direction, N 45 is the number of plies oriented at 
±45 degrees to the 0 direction, N T is the total number of plies, and a L is the longitudinal 
conductivity of a unidirectional laminate. This expression assumes that the laminate is an 
electrically homogenous material, the conductivity is not a function of thickness or 
stacking sequence, the plies perpendicular to the 0 direction are negligible, and that each 
±45° oriented ply with respect to the 0 direction is equivalent to a single ply in the 0 
direction. This expression was verified with actual tests with good agreement. 
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The conductivity in the through-thickness direction is theoretically the same as the 
transverse conductivity. An additional model, appropriate to high fiber conductivity to 
matrix conductivity ratios, is based on the longitudinal conductivity, and approximated 
by the following expression: 50 

= V,' a L (74) 

where (3 is defined as: 

P = ^ (75) 

p D f 

where hp ly is the thickness (m) of the ply, and D f is the diameter (m) of the fibers. 
Equation 74 assumes that the conductivity is not a continuum property, and that the 
matrix is an insulator and its conductivity is zero. This equation is significantly different 
from the equations presented above for the longitudinal and transverse conductivity. The 
main reason for this difference is that the equation for the through-thickness conductivity 
was developed using an empirical approach with measured experimental data which 
showed that the conductivity depends mainly on the fiber conductivity and the amount of 
fiber contact (i.e. percolation). It should be noted that equation 74 underestimates the 
experimental results of through-thickness conductivity, especially when used on thick 
plies (more than 10 fiber diameters thick). 

Percolation occurs when sufficient fibers exist in the composite as to create a path 
of touching fibers from the top surface of the composite to the bottom surface (refer to 
Figure 2.6). Percolation theory predicts that below a certain fiber volume fraction, the 
transverse conductivity is negligible compared to the longitudinal conductivity, as there 
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are not sufficient fibers to link up together to create a percolation path. Thus the 
conductivity changes by several orders of magnitude after crossing the percolation 
threshold. This critical volume fraction is defined theoretically as 0.45 for unidirectional 
laminates. 51 Joy et al 51 ' 53 found that transverse ply conductivity decreases with increasing 
thickness up to a thickness of 1 50 fiber diameters, above which it was considered to be of 
infinite thickness. The higher conductivity at small thicknesses is due to additional 
conduction paths. Therefore, there is theoretically a critical ply thickness, above which 
the transverse conductivity becomes negligible as compared to the longitudinal 
conductivity for a certain fiber volume fraction. 


It is interesting to note the widely varying results from different micromechanics 
models when the conductivity of the fibers is much greater than the conductivity of the 
matrix, as in carbon fiber / epoxy composites. If the assumption that the conductivity of 
the matrix is effectively equal to zero is extended to the previous models, equation 69 
becomes: 

a L =a f V f (76) 


and equation 70 becomes: 

ct t = (2 V f - l) CT f (77) 

This is a “weak” percolation based model which shows that the transverse conductivity 
becomes zero below a critical volume fraction (V f ) of 0.5, and rises linearly after that. 


Finally combining equation 74 and 76 gives: 


° thick = Vf P+ ' cr f 


(78) 
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which is a “stronger” percolation based model. These three models are plotted in Figure 
4.13, the fiber conductivity used is 1,000,000 1/Q-m, the fiber diameter used is 5 pm, and 
the ply thickness used is 0. 1 25 mm. 

4.5.2 Particulate Filled Composites 

In particulate filled composites, particles are used as the reinforcing medium 
instead of continuous length fibers. The particulates can be short chopped fibers, flakes 
and microballoons, or powders. One of the main uses of these composites in to create 
conductive composites out of insulative matrices. 


Percolation phenomena in particulate filled composites occurs when sufficient 
conductive filler is added to an insulating polymer matrix, such that the composite 
transforms from an insulator to a conductor. This occurs because there is sufficient filler 
to create a continuous linkage of filler particles from one edge to another, as shown in 
Figure 4.14. The critical volume fraction, V c , is when the change occurs. 


The critical volume fraction can be approximated by the packing fraction of the 
filler powder: 54 

V c *4> p -0.05 (79) 

where the packing fraction (<J> P ) can be calculated using: 


*)> p — 


m 


nn 


Volf, || Pmi 


(80) 
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Figure 4.13: Various conductivity micromechanical models 
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Figure 4.14: Critical fiber volume fraction in conductive composites 54 
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where m nn is the measured mass (kg) of the filler powder, Vol nn is the volume (nr 3 ) of 
tapped filler powder, and p fiM is the theoretical density (kg/m 3 ) of the filler powder. 

Therefore the rules of particle packing can be applied to the critical volume 
fraction, allowing particle filled composites to be thought of as conductive skeletons. 
Fine particles pack less densely than larger particles, and thus result in lower critical 
volume fractions. This occurs due to an increase in surface area, lower particle mass, and 
a greater significance of the weak short range attractive/repulsive forces. Therefore, the 
critical volume fraction increases with particle size. Also, the more spherical the powder, 
the less interparticle friction and the higher the critical volume fraction. Fiber and flake 
fillers do not pack densely due to extensive bridging, and thus their packing fraction and 
critical volume fraction decrease as the aspect ratio of the filler increases. 

In particulate composites, the electrical conductivity is not only a function of the 
constituent conductivities and the filler volume fraction, but also of the shape or aspect 
ratio of the particulates, and the relative particle arrangement. 

4.6 DATA REDUCTION EQUATIONS 

This section will present the equations which will be used to reduce the data from 
the experimental analysis. The geometry of the electrical test sample is as shown in 
Figure 4.15, where A is the cross-sectional area (m 2 ), and h is the thickness (m). The 
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Figure 4.15: 


Electrode 



Geometry of the electrical test sample 
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geometry of the mechanical test sample is as shown in Figure 4.16, where 1 is the length 
(m), and A is the cross-sectional area (m 2 ). 


4.6.1 Conductivity 

The conductivity (1/Q-m) of the sample is calculated as: 

h 

a ~ RA 

where R is the resistance (Q) across the sample. 


(81) 


The high resistance samples act as RC circuit, as shown in Figure 4. 1 7, where C is 
the capacitance (F) of the sample, and R is the resistance (Q) of the sample. The voltage 
in the circuit as the capacitor discharges is given by: 

t 

V(t) = V 0 e * (82) 

where V(t) is the voltage (V) at time t, V 0 is the initial voltage (V), and x is the RC circuit 
time constant (s), defined as: 

x = RC (83) 


To calculate the conductivity of the surface layer, equation 81 is used replacing 
the sample thickness (h) with the surface layer thickness (h SL ): 


a = 


SL 

RA 


(84) 
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Figure 4.17: 



RC circuit model of high resistance conductivity samples 
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4.6.2 Dielectric Constant 


The dielectric constant of the sample is defined as: 

Ch 

K = 

e 0 A 


( 85 ) 


4.6.3 Dielectric Strength 

The dielectric strength (V/m) of the sample is defined as: 

E = — (86) 

where V BD is the voltage applied to breakdown the material (V). To calculate the 
dielectric strength of the surface layer, equation 86 is used replacing the sample thickness 
(h) with the surface layer thickness (h SL ): 

E =^D (87) 

max v v ' 

hsL 


4.6.4 Density 

The density (kg/m 3 ) is defined as: 


P mech 


m 

Vol 


( 88 ) 


where m is the mass (kg) of the sample, and Vol is the volume (m 3 ) of the sample. 


4.6.5 Young’s Modulus 

The Young’s modulus (MPa) is defined as: 

p ^ mech ^ 

^ mech “ " Ac 

mech & mech 


(89) 
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where a mech is the stress (MPa) applied to the material, e mech is the strain due to the applied 
stress, and P is the load (N) applied to the material. 

4.6.6 Poisson Ratio 

The Poisson ratio is defined as: 

v = - -*=*- (90) 

£ x 

mcch 

where e y ^ is the strain perpendicular to the direction of the applied load, and e x is 
the strain in the direction of the applied load. 


4.6.7 Failure Stress 

The failure stress (MPa) is the maximum stress in the material before failure, it is 
defined as: 


a 


maxmcd, 



where P fajl is the load (N) applied to the sample at failure. 


(91) 


4.6.8 Failure Strain 

The failure strain is the strain when the material fails, it is defined as: 

_ Al 

Em “ mC d, — j 

where Al is the change in length, stroke, (m) of the sample at failure. 


(92) 
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CHAPTER 5 


ANALYTICAL RESULTS AND DISCUSSION 


This chapter presents the result from the deep dielectric charging model and the 
fiber-matrix interaction models. A base case was used with the deep dielectric charging 
model, from which all other runs are just modifications. A summary of the results are 
presented here. The plots from individual code runs can be found in Appendix B. 

5.1 DEEP DIELECTRIC CHARGING MODEL 

5.1.1 Base Case Results 

In this section, results are presented for a typical case, the base case from which 
all other code runs are variations. The orbit used has an altitude of 7,000 km 
(approximately one earth radii), an inclination of 0°, and was set during the solar cycle 
condition of solar maximum. This orbit is not a particularly practical orbit. It was 
chosen for two reasons. It is in the middle of the Van Allen belts (refer to Figure 2.2), 
where the highest concentrations of electrons are located, and geostationary spacecraft 
must fly through this region to get to geosynchronous orbit. The orbit environment is 
shown in Figure 4.3, it has a total electron flux of 2.8 x 10 s #/cm 2 -s with a total dose rate 
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of 7.15 rad/s, and a total proton flux of 2.9 x 10 7 #/cm 2 -s with a total dose rate of 771 


rad/s. 


The material properties used in the base case are those for a typical organic 
polymer, and are shown in Table 5.1 . The base case was run using both the closed form 
solution and the finite difference solution, without including radiation induced 
conductivity and the high field conductivity. The two boundary conditions used were an 
electric field and voltage of zero on the back surface. The electric field of zero on the 
back surface represents the back surface being a conductively grounded backplane. The 
results shown in Figures 5.1 through 5.3 for the closed form solution and Figures 5.4 
through 5.6 for the finite difference solution, show that both solutions give similar results 
with the finite difference solution taking slightly longer to reach steady state than the 
closed form solution. This is due to the definition of when the finite difference solution 
has converged— in this case when the residual is less than 0.1. The residual is defined as 
the summation of the difference between the incoming charge density rate and the current 
density for every node. If a larger residual was chosen for when steady state is reached, 
the time to steady state would be lower. This validates both models since they give the 
same result distributions, with a slight difference in the time to steady state. The results 
show that it took an average of 1 .3 seconds to reach steady state, with a maximum electric 
field of approximately 4400 V/m and a maximum voltage of approximately 3 volts. The 
electric field result is several orders of magnitude below the breakdown value of 15 
MV/m, therefore an electrostatic discharge will not occur. 
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Table 5.1: Base case material properties 


Property 


Value 


Density, p mech 

Thickness, h 
Dielectric Constant, k 
Electrical Conductivity, a 
Dielectric Strength, E max 
Radiation induced conductivity 
exponent, A 

Coefficient of radiation induced 
conductivity, k R 

Temperature, T 


1 600 kg/m 3 
2.5 mm 
3.0 

1 x 10' 10 1/Q-m 
15 x 10 6 V/m 


1.0 


1.1 x 10' 14 s/Q-m-rad 


120 °C 



Charge Density (10 15 #/m 3 ) 


Time (sec) 



0.122 

0.245 

0.367 

0.489 

0.734 

0.979 

1.223 


*1 I 


Depth (mm) 


Figure 5.1 : Closed form base case charge density 
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Figure 5.2: Closed form solution base case electric field 






Voltage (V) 



Time (sec) 

0.122 

0.245 

0.367 

0.489 

0.734 

0.979 

1.223 


I- ' 

2 


Depth (mm) 


Figure 5.3: Closed form solution base case voltage 
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Figure 5.4: Finite difference solution base case charge density 
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Figure 5.6: Finite difference solution base case voltage 


5.1.2 Base Case With High Field Conductivity 


The base case was run including high field conductivity. The temperature used 
was 120 °C, as indicated in Table 5.1. The temperature was chosen arbitrarily within the 
possible range of spacecraft temperatures. The analysis showed that the high field 
conductivity scaling factor has the same form as the electric field profile, as shown in 
Figure 5.7. However, the variations are quite small and the scaling factor can be 
approximated as 1. Thus the results showed that there is no appreciable change to the 
shape of the profiles, and the steady-state charge density, electric field, and voltages were 
unchanged. It should also be noted that the time to steady state did not change. 

Since the high field conductivity had no effect for the given steady-state electric 
fields of the base case, sensitivity studies were performed on the high field conductivity 
equation to determine the effects of dielectric constant and temperature. The range of 
values used are shown in Table 5.2. The dielectric constant was varied through the range 
of possible values for all materials, from 1 to 9, as shown in Figure 5.8. The temperature 
was varied through the possible range of temperatures seen by spacecraft in Earth orbit, 
from -200 to 150 °C, as shown in Figure 5.9. The results show that there is no effect due 
to high field conductivity until an electric field of at least 100,000 V/m, after which the 
high field conductivity scaling factor increases rapidly. Also the high field conductivity 
scaling factor decreases significantly with increasing dielectric constant and temperature. 
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High Field Conductivity Scaling 
Factor 



Depth (mm) 


Figure 5.7: High field conductivity scaling factor 



Table 5.2: Sensitivity studies parameter values 


Property 

Range 

Conductivity, a 

1 0‘ 8 to 10' 12 1/Q-m 

Dielectric Constant, k 

1.0 to 9.0 

Thickness, h 

0.01 to 5.0 mm 

Density, p mech 

600 to 2400 kg/m 3 

Temperature, T 

-200 to 150 °C 

Radiation Induced Conductivity 
Exponent 

0.5 to 1.0 

Coefficient of Radiation Induced 
Conductivity 

10’ 15 to 1.1 x 10' 14 s/Q-m-rad 

Solar Cycle Condition 

Minimum & Maximum 

Substorm Intensity 

1 to 1000 

Solar Paticle Event Intensity 

1 to 1000 
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Figure 5.8: High field conductivity scaling factor dielectric strength sensitivity 
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Figure 5.9: High field conductivity scaling factor temperature sensitivity 
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5.1.3 Base Case with Radiation Induced Conductivity 


The base case was run including radiation induced conductivity, the coefficient of 
radiation induced conductivity (k) and the radiation induced conductivity exponent (A) 
where selected (refer to Table 5.1) to produce the worst case possible (the highest 
possible radiation induced conductivity within the range of values for polymers and 
composites). The results showed no change to the shape of the profiles, and a slight 
decrease in the steady-state charge density, electric field, and voltage as compared to the 
base case. The reason for the slight decrease is that the radiation induced conductivity is 
9.2 x 10' 12 1/Q-m which is an order of magnitude less than the dark conductivity of the 
material. Thus when it is added to the dark conductivity, the conductivity increases 
slightly to 1.09 x 10' 10 1/Q-m. Therefore since the conductivity increased, the time to 
steady state also decreases slightly from 1.22 seconds to 1.12 seconds. It should be noted 
that even though the dark conductivity of the material could decrease below the value of 
the radiation induced conductivity, the value obtained for the radiation induced 
conductivity is a worst case value and in reality it may even be lower. Also, the effect of 
radiation induced conductivity is to decrease the value of the electric field and thus helps 
in preventing electrostatic discharges from occurring, so it is conservative to ignore it. 

To observe the effects of the coefficient of radiation induced conductivity and the 
radiation induced conductivity exponent, sensitivity studies were performed on the 
radiation induced conductivity equation. The range of values used are shown in Table 
5.2. The radiation induced conductivity exponent was varied through all possible values, 
0.5 to 1.0, and the coefficient of radiation induced conductivity was varied through the 
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range of possible values for composites and polymers, from 10' 15 to 1.1 x 10 14 s/Q-m-rad. 
The results, as shown in Figures 5.10 and 5.11, show that the radiation induced 
conductivity increases with dose rate, the radiation induced conductivity exponent. A, and 
the coefficient of radiation induced conductivity, k. 

5.2 SENSITIVITY STUDIES 

To understand the effects of the various parameters input into the deep dielectric 
charging model, sensitivity studies were performed on the material properties: the 
conductivity, dielectric constant, thickness, and material density. Sensitivity studies were 
also performed on the space environment: solar cycle condition, geomagnetic substorms, 
and solar particle events. 1 he range of values used in each sensitivity study is shown in 
Table 5.2. No sensitivities were performed on the dielectric strength even though it is 
one of the input material parameters. This is because the only role the dielectric strength 
plays is to determine when an electrostatic discharge will occur. Therefore changing the 
dielectric strength will only have an effect if it is decreased below the maximum steady- 
state electric field, thus allowing an electrostatic discharge to occur. 

5.2.1 Conductivity 

The conductivity was varied within the range of insulating composites and 
polymers, from 10' 2 to 10‘ 8 1/Q-m, as shown in Figure 5.12. The result show that both 
the maximum steady-state electric field and the time to steady state scale inversely with 
conductivity. Also there was no change in the shape of the charge density, electric field, 
and voltage distributions, they were just decreased by the inverse of the conductivity. 
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These results are expected, as can be seen in equations 31, 32, 34, and 36 where the 
conductivity is in the denominator. 

5.2.2 Dielectric Constant 

The dielectric constant was varied within the range of possible values for 
materials, from 1 to 9, as shown in Figure 5.13. The results show that the maximum 
steady-state electric field is unchanged with increasing dielectric constant, and the time to 
steady state increases linearly with dielectric constant. Also there was a slight change in 
the shape of the charge density distribution, and no change at all to the electric field and 
voltage distribution. These results are expected, as can be seen in equations 31 and 32 
where the dielectric constant is in the numerator, and in equations 34 and 36 which show 
that the electric field and voltage are not a function of the dielectric constant. 

5.2.3 Thickness 

The thickness was varied within the range of possible structural thicknesses, from 
0.01 to 5.0 mm, as shown in Figure 5.14. The results show that the maximum steady- 
state electric field increases with thickness until a constant value is obtained after which 
the electric field does not change with increasing thickness. The time to steady state is 
unchanged with increase in thickness. There were also changes to the charge density, 
electric field, and voltage distributions. These results occur because at small thicknesses 
there are less particles in the material, as the higher energy particles have enough energy 
to pass through the material. While as the thickness increases, more particles are stopped 
in the material increasing the electric field. This continues until the material is thick 
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Figure 5.13: Dielectric constant sensitivity 
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enough to stop the majority of the particles, and thus the maximum electric field does not 
increase with any further increase in thickness. 

5.2.4 Density 

The density was varied within the range of possible densities for composites and 
polymers, from 600 to 2400 kg/m 3 , as shown in Figure 5.15. The results show that the 
maximum steady-state electric field increases with density until it reaches a constant 
value. The time to steady state is unchanged with increase in density. Also, the charge 
density, electric field, and voltage distributions change with density. The explanation for 
the electric field trend is similar to that for the thickness sensitivities, since the density of 
the material controls how far a particle penetrates for a given energy, as shown in Figure 
5.16. Therefore, at low densities the higher energy particles pass through the material, 
thus resulting in a lower electric field. As the density increases, more particles are 
stopped by the material and the maximum electric field increases, until the density is such 
that the majority of the particles are stopped and the maximum electric field is unaffected 
by further increases in density, as shown in Figure 5.17. 

5.2.5 Solar Cycle Condition 

The base case was run with the environment of the same orbit but at the solar 
cycle condition of solar minimum. As can be seen in Figures 5.18, the solar minimum 
condition does not affect the shape of the flux profiles, it reduces the electron flux by 
about 3.5 times and does not affect the proton flux. This results in the decrease of the 
negative charge density, the maximum electric field, and the maximum voltage by 
approximately the same amount, as shown in Figures 5.19 through 5.21. However, since 
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Figure 5.15: Density sensitivity 
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Figure 5.17: Steady-state electric field for various densities 
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Figure 5.18: Base case orbit environment 
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Figure 5.19: Solar condition charge density sensitivity 
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Figure 5.20: Solar condition electric field sensitivity 
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Figure 5.21 : Solar co 




the proton flux was not affected, most of the positive charge density was not affected. 
Thus there is no change to shape of the electric field near the front surface (Depth (x) = 0 
mm). It should also be noted that the time to steady state did not change between solar 
maximum and solar minimum; this is because the time to steady state is only dependent 
on the conductivity and the dielectric constant of the material. 

5.2.6 Geomagnetic Substorm 

As mentioned earlier, geomagnetic substorms increase the amount of electron flux 
in orbit. Therefore the substorm intensity factor is used to scale the incoming electron 
flux. The intensity factor was incremented from 1 to 1000, representing no substorms up 
to a worst case substorm, as shown in Figure 5.22. The results show that the maximum 
steady-state electric field increased linearly with the substorm intensity factor, while the 
time to steady state remained unchanged. Also, the charge density, electric field, and 
voltage distribution changed with the substorm intensity factor, as while the number of 
electrons was increasing the number of protons remained constant. Thus as the intensity 
factor increased, the effect of the protons becomes more and more negligible. This is 
shown in Figure 5.23, where the increase in electric field due to protons near the front 
surface (refer to Figure 5.2) is unobservable. 

5.2.7 Solar Particle Event 

As mentioned earlier, solar particle events increase the proton flux in orbit. 
Therefore the solar particle event intensity factor is used to scale the incoming proton 
flux. The intensity factor was incremented from 1 to 1000, representing no solar particle 
event up to a worst case solar particle event, as shown in Figure 5.24. The results show 
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Figure 5.23 : Electric field of base case with a substorm intensity of 1 000 
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Figure 5.24: Solar particle event sensitivity 
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that the maximum steady-state electric field remains constant until the intensity factor is 
high enough to start increasing the magnitude of the (negative) electric field. The time to 
steady state remains unchanged with increasing intensity factor. Also, the charge density, 
electric field, and voltage distributions changed with solar particle event intensity factor, 
as while the number of protons increased, the number of electrons remained constant. 
The reason for the initial constant electric field is the boundary condition of zero electric 
field on the back surface. With no solar particle event, the protons cause an increase in 
electric field of approximately 400 V/m, while the electrons cause a decrease in electric 
field of approximately 4400 V/m. Therefore since the electrons penetrate further than the 
protons, to obtain an electric field of zero on the rear surface, there must be an electric 
field of approximately 4000 V/m on the front surface. As the proton flux increases, it 
creates a greater increase in electric field, however until the increase is greater than 
approximately 4000 V/m, the only change is to decrease the electric field on the front 
surface, the maximum electric field remains the same. As the intensity factor increases, 
the influence of the electrons becomes more and more negligible, as shown in Figure 
5.25, where the decrease in electric field due to the electrons is unperceivable. 

5.2.8 Summary 

The time to steady state depends inversely on the conductivity and linearly on the 
dielectric constant. The steady-state electric field is inversely dependent on the 
conductivity and dependent on the amount of particles present in the material. The 
amount of particles present in the material is controlled by the sample thickness, material 
density, and the substorm and solar particle event intensity factors. The shape of the 
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Figure 5.25: Electric field of base case with a solar particle event intensity of 1 000 



charge density, electric field, and voltage distributions are dependent on the location of 
particles and on the relative amounts of protons and electrons. The location of particles is 
controlled by the sample thickness and the material density, and the relative amounts of 
particles is controlled by the substorm intensity factor and the solar particle event 
intensity factor. 

5.3 MICROMECHANICAL DETAILS 

This section presents the results of the modeling of the interactions between fibers 
and matrix on the electric field. The two models examined are the micro-scale model of a 
fiber surrounded by matrix, and the insulative pure epoxy surface layer of carbon fiber / 
epoxy composites. 

5.3.1 Micro-scale Model 

At the micro-scale, a concern is that higher fields can be possible due to geometric 
details of the fibers imbedded in matrix. The micromechanics of the fiber surrounded by 
epoxy has been looked at, to see if a substantial increase in the electric field can occur as 
the electric field in the epoxy matrix approaches the highly conductive carbon fibers. 

The data used in this model comes from the base case used with the deep 
dielectric charging model. The far field electric field used is the maximum steady-state 
electric field from the base case of 4326 V/m. The far field electric field is approximately 
6 to 8 times the fiber radius, this represents a worst case matrix rich zone. For ease of 
display the results, as shown in Figure 5.26, are plotted as electric field magnification 
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Figure 5.26: Electric field magnification around a fiber 
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versus normalized radius. The electric field magnification is the calculated electric field 
divided by the far field electric field, and the normalized radius is the radius divided by 
the fiber radius. 

Since the magnitude of the electric field is used the results are quadrant 
symmetric, the results are identical above and below the x axis and to the left and right of 
the y axis. The results show that the electric field increases in the matrix as the field 
approaches the fiber parallel to the far field, with the maximum occurring at the 
fiber/matrix interface. When the field approaches the fiber perpendicular to the far field 
direction the electric field decreases, and reaches zero at the fiber/matrix interface. Thus 
the greatest possible increase is approximately a factor of 3. It should be noted that this 
solution is an approximation to the problem as the charge density is assumed to be zero. 
Also, the typical distance between fibers for a 0.6 volume fraction carbon fiber / epoxy 
composite is approximately 1 .3 times the fiber radius, thus invalidating the far field 
electric field assumption. 

Another related problem is that unlikely arrangements of multiple fibers may 
create severe distortions in the electric fields, as shown in Figure 5.27, possibly resulting 
in increases as high as 50 times. 55 It should be noted that the electric field is small 
enough that even an increase of 50 times would not be sufficient to create an electrostatic 
discharge in the base case. However, in other cases where the electric field is close to the 
dielectric strength of the material an increase of 50 times could cause an electrostatic 
discharge to occur. 
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Figure 5.27: Schematic of an unlikely fiber arrangement which can lead to high field 

magnification 
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5.3.2 Epoxy Surface Laver 

A related problem is that of the electric field in carbon fiber / epoxy with an epoxy 
rich surface layer. This layer, approximately 20 pm thick, has a much lower (several 
orders of magnitude) electrical conductivity than the bulk composite. The analysis was 
accomplished by modifying the deep dielectric charging analyzer to include two different 
conductivities, the dark conductivity of the bulk composite and the surface layer 
conductivity. 

The results showed that the electric field in the bulk composite reached steady 
state as it would have in the case of only carbon fiber / epoxy with no epoxy surface 
layer, a maximum steady-state electric field of 4.33 x 10' 6 V/m and time to steady state of 
1.22 x 10‘ 9 sec. However, the electric field in the surface layer kept increasing. To solve 
for the electric field and steady-state time in the surface layer, the R-C circuit model is 
used (refer to Section 4.4.2), with the data shown in Table 5.3. A simple R-C circuit 
analysis can be used since there is such a difference in the time scales of the bulk 
composite and the epoxy surface layer. The results show that the time to steady state is 
1.22 sec and the steady-state electric field across the surface layer is 473 V/m. These 
values are very close to the values that would have been obtained if the entire sample was 
made of epoxy, i.e. the base case described earlier. 
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Table 5.3: Data used in epoxy surface layer model 


Property 

Value 

Surface layer conductivity, a 

10' 10 1/Q-m 

Bulk composite conductivity, a 

0.1 1/Q-m 

Dielectric constant, k 

3 

Surface layer thickness, h SL 

25 pm 

Incoming charge density, p ln 

1.18 x 10 16 #/m 3 -s 
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5.4 CASE STUDIES 


Case studies were performed for the base case material configuration on actual 
orbits used by spacecraft where deep dielectric charging may be of concern. These orbits 
include geosynchronous (GEO) orbit, the Global Positioning System (GPS) orbit, a low 
altitude retrograde or polar orbit, the Molniya orbit, and the space shuttle parking orbit. 
The environment data was once again obtained from the Environmental Workbench 
software for both solar maximum and solar minimum solar cycle conditions. The orbit 
parameters used are shown in Table 5. 4. 56 

5.4. 1 The Geosynchronous Orbit 

According to the Environmental Workbench, at geosynchronous orbit there is no 
distinction between solar maximum and solar minimum. The electron and proton fluxes 
can be seen in Figures 5.28. While the electron flux distribution looks similar to the base 
case, the proton flux does not have the high energy (greater than 1 MeV) particles. The 
results, as shown in Figures 5.29 through 5.31, are very similar to the base case except 
that they are slightly less than an order of magnitude smaller. Also, the electric field and 
voltage distributions penetrate deeper into the material before becoming negligible. 

5.4.2 Global Positioning System Orbit 

The Global Positioning System orbit environment, as shown in Figures 5.32, has 
approximately half the electron flux at solar minimum as compared to solar maximum, 
and the electron flux has higher energy particles, up to 3 MeV. Also, the proton flux is 
unchanged between solar maximum and solar minimum, and once again it does not have 
the higher energy particles like the base case does. The results, as shown in Figures 5.33 
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Table 5.4: Case study orbit data 


Orbit 

Apogee Altitude 

Perigee Altitude 

Inclination 

Geosynchronous 

35,782 km 

35,782 km 

0° 

GPS 

20,222 km 

20,222 km 

60° 

Polar 

322 km 

322 km 

98° 

Molniya 

40,172 km 

272 km 

63.4° 

Shuttle Parking 

322 km 

322 km 

28.4° 
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Figure 5.28: Geosynchronous orbit environment 
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Figure 5.29: Geosynchronous orbit charge density distribution 
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Figure 5.30: 
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Figure 5.3 1 : Geosynchronous orbit voltage distribution 
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Figure 5.33: Global Positioning System orbit charge density distribution 
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Figure 5.34: Global Positioning System orbit electric field distribution 
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Figure 5.35: Global Positioning System orbit voltage distribution 
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through 5.35, show that the charge density distributions for the two solar conditions are 
very similar to each other. Also, electric field and voltage distributions once again 
penetrate deeper into the material as described in the geosynchronous orbit case. It is 
interesting to note is that the ramp up portion of the electric field distribution for the solar 
minimum case is the same as that for the solar maximum case except that it is shifted 
down by approximately 200 V/m. 

5.4.3 The Low Altitude Retrograde Orbit 

The low altitude retrograde orbit environment, as shown in Figures 5.36, has a 
typical electron flux profile with a few high energy particles, and the solar maximum 
values are approximately 2.5 times greater than the solar minimum values. The proton 
flux is more interesting; the solar minimum values are greater than the solar maximum 
values, as expected for low altitude orbits. Also, the distributions are dissimilar; there are 
proportionally more high energy particles at solar minimum than at solar maximum. 
Another interesting fact is there is approximately 3 orders of magnitude more electrons 
than protons, compared to the previous orbits where there is only at best a 1 order of 
magnitude difference. 

The results, as shown in Figures 5.37 through 5.39, show that even though there 
are more protons at solar minimum, the steady-state charge density, electric field, and 
voltage distributions are still larger at solar maximum. This occurs because although 
there are more protons at solar minimum, they are still 3 orders of magnitude less than the 
number of electrons. This can be noticed in the charge density distribution where the 
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Figure 5.36: Low altitude retrograde orbit environment 
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Figure 5.37: Low altitude retrograde orbit charge density distribution 
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Figure 5.38: Low altitude retrograde orbit electric field distribution 
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Steady-state Voltage (mV) 





positive charge density near the front surface due to the protons is essentially zero. It is 
also the reason for the constant electric field near the front surface. As in previous case 
studies where there are high energy electrons (greater than 1 MeV) the distributions are 
none zero through almost all of the thickness. Also, because of the lower fluxes present 
in low altitude orbits all the maximum values are significantly lower than previous cases. 

5.4.4 The Molniya Orbit 

The Molniya orbit environment, as shown in Figures 5.40, has no exceptional 
features. The proton flux is within an order of magnitude of the electron flux. The 
electron flux has 3 times the particles at solar maximum as compared to solar minimum. 
The proton flux is the same for solar maximum and solar minimum, and has a relatively 
small amount of high energy particles. The results, as shown in Figures 5.41 through 
5.43, exhibit the typical behavior seen earlier, the steady-state charge density, electric 
field, and voltage distributions are greater at solar maximum than at solar minimum by 
approximately a factor of 3 representing the difference in electron flux. Once again, as in 
the case studies with few high energy electrons the electric field and voltage distributions 
penetrate deeper into the material before becoming negligible. 

5.4.5 The Space Shuttle Parking Orbit 

The space shuttle parking orbit environment, as shown in Figures 5.44, has a 
typical electron flux distribution, the solar maximum values are approximately 2 times 
greater than the solar minimum values, and no high energy particles. The only difference 
to the base case environment is that the flux levels are much smaller, as would be 
expected in a low altitude orbit. The proton flux distributions are very different from 
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Figure 5.41 : Molniya orbit charge density distribution 
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Figure 5.42: Molniya orbit electric field distribution 
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Figure 5.43: Molniya orbit voltage distribution 
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Figure 5 .44: Space shuttle parking orbit environment 
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anything seen earlier. The distribution is not contained within the usual energy 
boundaries of 0.01 to 100 MeV. However, this is not a concern since as shown in Figure 
4.6, particles of these energies will penetrate over 100 mm of the base case material, and 
our sample material is only 2.5 mm thick. Therefore, these missing extremely high 
energy particles will just pass through the material. Also, as expected with low altitude 
orbits, the proton flux at solar minimum is much greater than at solar maximum. As with 
the low altitude retrograde orbit, the electron flux is several orders of magnitude greater 
than the proton flux. 

The results, as shown in Figure 5.45 through 5.47, show than once again since 
there are many more electrons than protons, the fact that there are more protons at solar 
minimum than maximum has no consequence. The steady-state charge density, electric 
field, and voltage distributions are greater for solar maximum than solar minimum. The 
fact that there are far greater amounts of electrons than protons can also be seen in charge 
density distribution where the positive portion near the front surface, due to the protons, 
is almost zero. It can also be seen in the electric field distribution where near the front 
surface it is constant. It should also be noted that since there are no high energy 
electrons, the electric field and voltage distribution reach essentially zero before the back 
surface of the material. Also, since the flux levels are very low in low altitude orbits, the 
maximum values of the results are also very low, as compared to the base case. 
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Figure 5.45: Space shuttle parking orbit charge density distribution 
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Figure 5.46: 


Space shuttle parking orbit electric field distribution 
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Figure 5.47: Space shuttle parking orbit voltage distribution 


169 




5.4.6 Summary 


Even though all the case studies showed that the electric field was not large 
enough to create an electrostatic discharge in the base case material, this does not mean 
that deep dielectric charging is not a concern. An extremely high intensity geomagnetic 
substorm or solar particle event could increase the flux levels sufficiently to create an 
electrostatic discharge. Also, other materials with either a lower conductivity or lower 
dielectric strength could also be susceptible to electrostatic discharges due to deep 
dielectric charging. The worst orbits are the geosynchronous orbit, the Global 
Positioning System orbit, and the Molniya orbit. These all produced electric fields in the 
100’s of volts per meter for the base case material. 

5.5 INTERNATIONAL SPACE STATION 

This section deals with an actual problem. There was a concern that the Kapton 
layer covering the main power bus of the International Space Station would be affected 
by deep dielectric charging. The power bus is used to connect electrically the solar cell 
wing to the space station. 

The basic problem is that during operation there is a differential of 160 volts 
across the Kapton Layer. The space environment is assumed to be at plasma ground (i.e. 
0 volts) and the power bus is carrying 1 60 volts from the solar cells to the space station. 
The space station orbit used had a constant semi-major axis of 6771 km, this is equivalent 
to an apogee and perigee altitude of 393 km, with an inclination of 51.6 0 . 57 Since the 
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exact type of Kapton being used was unknown, an analysis was performed on three types 
of Kapton, type HN, VN, and FN films. 58 The properties used are summarized in Tables 

5.5 and 5.6. 


The composite deep dielectric charging analyzer was run on each of the different 
types of Kapton listed in the above mentioned tables to determine which type would 
produce the worst results, the type that was most likely to undergo an electrostatic 
discharge. Once these were done, the worst case material type was analyzed to determine 
the effect of changing various parameters such as the solar cycle condition, including 
radiation induced conductivity, including high field conductivity, a negative voltage 
applied through the power bus, and no voltage applied through the power bus. Table 5.7 
shows the settings used during each test case. 

The space station orbit environment data, as shown in Figures 5.48, has typical 
flux distributions for low altitude orbits. However, it should be noted that unlike the low 
altitude cases investigated earlier, the higher proton flux during solar minimum had a 
noticeable effect on the results. This is because the charging in entirely dominated by the 
protons as all the electrons pass straight through the Kapton film, due to its small 
thickness. 


The results, as shown in Figures 5.49 through 5.51 for Test Case No. 1, show that 
in each case where there is a voltage applied through the power bus, the problem is 
dominated by the voltage difference applied across the Kapton film. The voltage varies 
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Table 5.5: Kapton type HN and VN film properties 


Property 

100VN, HN 

200 VN, HN 

300VN, HN 

500 VN, HN 

Density (kg/m 3 ) 

1420 

1420 

1420 

1420 

Thickness (pm) 

25 

50 

75 

125 

Dielectric Strength (V/m) 

303 x 10 6 

240 x 10 6 

205 x 10 6 

154 x 10 6 

Dielectric Constant 

3.4 

3.4 

3.5 

3.5 

Resistivity (Q-m) 

1.5 x 10 15 

1.5 x I0 15 




Table 5.6: Kapton type FN film properties 


Property 

120 FN 616 

150 FN 019 

250 FN 029 

Density (kg/m 3 ) 

1530 

1670 

1570 

Thickness (pm) 

30 

38 

63 

Dielectric Strength (V/m) 

272 x 10 6 

197 x 10 6 

197 x 10 6 

Dielectric Constant 

3.1 

2.7 

3.0 

Resistivity (Q-m) @ 23 °C 

1.4 x 10 15 

2.3 x 10 15 

1 .9 x 10 15 
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Table 5.7: Space station test case analysis parameters 


Test Case 

Material 

B.C. 

Solar Cycle 

R.I.C. 

H.F.C 

1 

100 VN, 100 HN 

0, 160 V 

max 

no 

no 

2 

200 VN, 200 HN 

0, 160 V 

max 

no 

no 

3 

300 VN, 300 HN 

0, 160 V 

max 

no 

no 

4 

500 VN, 500 HN 

0, 160 V 

max 

no 

no 

5 

120 FN 616 

0, 160 V 

max 

no 

no 

6 

150 FN 019 

0, 160 V 

max 

no 

no 

7 

250 FN 029 

0, 160 V 

max 

no 

no 

8 

100 VN, 100 HN 

0, 160 V 

max 

yes 

no 

9 

101 VN, 100 HN 

0, 160 V 

max 

no 

yes 

10 

100 VN, 100 HN 

0,-160 V 

max 

no 

no 

11 

100 VN, 100 HN 

0,0 V 

max 

no 

no 

12 

100 VN, 100 HN 

0, 160 V 

min 

no 

no 


B.C.: Voltage boundary conditions 
R.I.C.: Radiation induced conductivity 
H.F.C.: High field conductivity 
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Figure 5.48: Space station orbit environment 
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Figure 5.49: Space station Test Case No. 1 charge density distribution 
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Electric Field (MV/m) 






Figure 5.5 1 : Space station Test Case No. 1 voltage distribution 
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linearly from 0 V at the front surface to 160 V at the rear surface. This results in the 
electric field being simply the change in voltage divided by the film thickness. The only 
effect of the incoming radiation is to perturb the electric field about the voltage induced 
field value. However, since the amount of incoming particles is so small compared to the 
base case, the perturbation is only approximately 100 V/m, which is insignificant against 
the voltage induced electric field in the millions of volts per meter. The constant portion 
of the electric field distribution near the front surface occurs because there are no particles 
deposited in the film there, as even the smallest energy protons penetrate at least 0.5 pm. 

When the solar cycle was changed to solar minimum, as shown in Figures 5.52 
and 5.53, the only observable change to the results was the charge density and the 
radiation induced electric field perturbation increased by approximately an order of 
magnitude. This occurred because, as mentioned earlier, during solar minimum there is 
approximately an order of magnitude more protons. The effect of including the radiation 
induced conductivity, using the base case radiation induced conductivity parameter 
values, was to slightly decrease the electric field perturbation and the charge density, this 
occurred since the conductivity only increased from 6.67 x 10' 16 1/Q-m to 7.20 x 10 16 
1/Q-m. The same result was seen for the high field conductivity using a temperature 
120°C, for the same reason, the conductivity increased by a factor 1.48 from 6.67 x 10 16 
1/Q-m to 9.85 x 10 16 1/Q-m. By applying a negative voltage of -160 V through the 
power bus, the only effect was to make the voltage induced electric field become positive 
of the same magnitude. The no voltage through the power bus case showed the effect of 
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Figure 5.52: 
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Figure 5.53: Space station electric field distribution 
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only the radiation induced electric field, as shown in Figure 5.54. In this case the voltage 
created a parabolic shape, starting and finishing at 0 V, as shown in Figure 5.55. 
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Figure 5.55: Space station power off voltage distribution 
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CHAPTER 6 


EXPERIMENTAL PROCEDURE 


This chapter presents the experimental procedures used in this research. The test 
matrices are explained, followed by information on the manufacturing process used to 
produce the composite samples. Finally the measurement techniques used are described. 

6.1 TEST MATRICES 

There are two separate experimental efforts in this research. The first part is the 
investigation of the conductivity and the dielectric strength of a carbon fiber / epoxy 
composite. The second part is an investigation into creating a conductivity tailorable 
composite material, using a fiberglass / epoxy composite. The two material systems 
used in this investigation are: Hercules AS4/3501-6 tape prepreg, and a woven fabric 
fiberglass / epoxy wet layup. The first material system, AS4/3501-6, is a unidirectional 
prepreg and has a first generation brittle (low strain-to-failure) 3501-6 matrix, which has 
been widely used in the aerospace industry. The second system is comprised of a 
unidirectional woven fiberglass cloth type 1543-38 prepared in a wet layup with Rutapox 
L20/SL resin. This system is being used in the aircraft industry by Grob Aerospace of 
Germany in their general aviation aircraft. 
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6.1.1 Composite Electrical Properties 

The through-thickness electrical properties of composites made of carbon fiber / 
epoxy prepreg were investigated. The composite parameters that were looked at were the 
laminate thickness, the laminate stacking sequence, and the sample cross sectional area, 
as shown in Table 6.1. Three different thicknesses were used, 4 plies, 8 plies, and 32 
plies, corresponding to approximately 0.5 mm, 1.0 mm, and 3.8 mm thick, respectively. 
Two different stacking sequences were used, unidirectional and quasi-isotropic. With 4 
plies a quasi-isotropic laminate cannot be produced, therefore [±45] s was used as an 
approximation. Therefore six different laminates need to be manufactured. The cross 
sectional areas used for the conductivity samples are 25.4 mm x 25.4 mm (1" x 1"), 25.4 
mm x 50.8 mm (1" x 2"), and 50.8 mm x 50.8 mm (2" x 2"), corresponding to 
approximately 645 mm 2 , 1290 mm 2 , and 2580 mm 2 . All the dielectric strength samples 
were 25.4 mm x 25.4 mm. 

One 152.4 mm x 177.8 mm (6" x 7") laminate of each layup was manufactured, 
and 3 samples of each geometry were used. Six different laminates, and 3 different 
sample areas were required for the conductivity investigation, resulting in 18 different 
sample types and a total of 54 samples, as shown in Table 6.2. For the dielectric strength 
investigation, six different laminates and samples with only one area were required, 
resulting in 6 different sample types and a total of 18 samples, as shown in Table 6.3. 
Each sample was measured 3 times for electrical properties. The thickness was measured 
in 5 different locations, and the area was calculated by measuring the length and width in 
3 different places, as shown in Figure 6.1 . 
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Table 6.1 : Layups used for the composite electrical properties measurements 


Thickness 

Stacking 

Unidirectional 

Sequence 

Quasi-Isotropic 

4 Plies 

[0] 4 

[±45] s 

8 Plies 

[0]« 

[0/±45/90] s 

32 Plies 

[0] 32 

[0/±45/90] 4S 


Table 6.2: Electrical conductivity measurement test matrix 


Area 

Ui 

[0] 4 

lidirectioi 

[0] 8 

nal 

(01,6 


Quasi-Isotrc 

[0/±45/901 s 

•pic 

[0/±45/90) 4S 

Total 

BO 

3 

3 

3 

3 

3 

3 

18 

OB 

3 

3 

3 

3 

3 

3 

18 


3 

3 

3 

3 

3 

3 

18 

Total 

9 

9 

9 

9 

9 

9 

54 


Table 6.3: Dielectric strength measurement test matrix 


Area 

Ur 

[0] 4 

lidirectioi 

[0] 8 

lal 

(01,6 


Quasi-Isotro 

[0/±45/90] s 

pic 

[0/±45/90] 4S 

Total 

l"x 1" 

3 

3 

3 

3 

3 

3 

18 
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6. 1 .2 Conductivity Tailorable Composite 


To investigate the development of a conductivity tailorable composite, a glass 
fiber / epoxy wet layup procedure was used. Carbon black was added to the epoxy in 4 
different percentage by mass: 5, 10, 15, and 20 percent. The carbon black used was 
Columbian Conductex Ultra High Conductivity carbon black. The carbon black filled 
epoxy resins are identified by Penn Color as RD 38663-50 for the 5%, RD 38662-50 for 
the 10%, RD 38661-50 for the 15%, and RD 38660-50 for the 20%. The electrical 
properties measured included the conductivity, dielectric constant, and dielectric strength, 
and the mechanical properties measured included the density. Young’s modulus, Poisson 
ratio, failure stress, and failure strain. Three 304.8 mm x 355.6 mm (12" x 14") 
laminates of each percentage were manufactured, and 2 samples were taken from each 
laminate. Therefore 6 samples of each percentage of carbon black were manufactured, 
for a total of 30 samples for each test. As shown in Table 6.4, the test matrix is fully 
populated. The size of the samples used for the electrical properties are 25.4 mm x 25.4 
mm (1" x 1 "). The density was measured using the conductivity samples. The size of the 
tensile test specimen used were 50.8 mm x 355.6 mm (2" x 14") with 50.8 mm x 76.2 
mm (2” x 3") tapered glass loading tabs at each end on both sides. Each sample was 
tested once. The geometry of the electrical property samples were measured as shown in 
Figure 6.1. For the tensile specimen the thickness was measured 9 times, the width 3 
times, and the length between tabs twice, as shown in Figure 6.2. 
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Table 6.4: Conductivity tailorable composite test matrix 


Property 


Percent Carbon Black By Mass 

Total 


0% 

5% 

10% 

15% 

20% 


Conductivity 

6 

6 

6 

6 

6 

30 

Dielectric Constant 

6 

6 

6 

6 

6 

30 

Dielectric Strength 

6 

6 

6 

6 

6 

30 

Density 

6 

6 

6 

6 

6 

30 

Young's Modulus 

6 

6 

6 

6 

6 

30 

Poisson Ratio 

6 

6 

6 

6 

6 

30 

Failure Stress 

6 

6 

6 

6 

6 

30 

Failure Strain 

6 

6 

6 

6 

6 

30 

Total 

48 

48 

48 

48 

48 

240 
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6.2 MANUFACTURING PROCEDURES 

The manufacturing procedure used for the AS4/3501-6 laminates is the TELAC 
net resin cure, 59 and is summarized here. Procedures for manufacturing the glass fiber / 
epoxy laminates was documented by Ed Wolf 60 and the manufacturer, and are also 
presented here. Two different layup techniques were used in this investigation: a prepreg 
layup for the carbon fiber / epoxy material, and a wet layup for the glass fiber / epoxy 
system. Latex gloves were worn whenever the raw materials were handled directly, and 
care was taken to avoid unnecessary contact with the uncured materials at all times. 

6.2. 1 Carbon Fiber / Epoxy Prepreg Layup 

The carbon fiber / epoxy prepreg tape came in the form of rolls 305 mm wide. 
The rolls of prepreg were stored in sealed bags in a freezer at a temperature lower than 
-18 °C. In preparation for layup, the prepreg was taken out of the freezer and allowed to 
sit at room temperature in the sealed bag for at least 1 hour, or until the material was no 
longer cold to the touch. If the bag had been opened while the material was still cold, 
unwanted moisture could have condensed on it. The prepreg was then unrolled and the 
appropriate plies were cut using Teflon coated aluminum patterns and a sharp utility 
knife. The plies were cut in such a way that any seams within a ply were parallel to the 
fiber direction. Once cutting was completed, the prepreg roll was resealed and placed 
back in the freezer. 

To make the smaller-than-standard laminates, 152.4 mm x 177.8 mm (6" x 7"), 
the prepreg was cut using the standard 304.8 mm x 355.6 mm (12" x 14") patterns. 
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Following this, each ply was cut into 4 pieces using small patterns. The individual plies 
were stacked in the proper sequence using an L-shaped aluminum jig to aid in alignment. 
The corner of the laminate which was situated directly in the comer of the aluminum jig 
was assumed to have the most accurate ply stacking and was marked for future reference. 
Between each application of a ply, the paper backing was removed from the laminate. 
When the laminate was completed, it was trimmed with a sharp utility knife and an 
aluminum pattern to its nominal 1 52 mm wide by 1 78 mm long size, and florocarbon 
release film was placed on both the top and bottom surfaces. If the laminates were to be 
cured within 24 hours, they were sealed in a vacuum bag and left out at room 
temperature. Otherwise, they were sealed in a bag and put into the freezer. 

In preparation for the cure, an aluminum cure plate was sprayed with Mold Wiz® 
release agent and covered with a sheet of nonporous Teflon. On top on this, three layers 
of cork dam were used to make six 152 mm by 178 mm (6" x 7") rectangular vessels for 
the laminates, followed by another cork dam made of one layer of cork, this second dam 
was used to create a step down from the 3 layer dam height. A sheet of nonporous Teflon 
was placed in the bottom of the vessel, to line the bottom and sides of the vessel. This 
was followed by the prepreg laminate with florocarbon release film on both sides, and a 
sheet of porous Teflon. The location of the “good” comer was marked. The aluminum 
top plate, sprayed with Frekote 700, and wrapped in nonporous Teflon, was then placed 
on top of the laminate. The assembly was then covered with a single sheet of porous 
Teflon, several layers of glass breather, and finally, the vacuum bag, which was sealed to 
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the caul plate with vacuum tape. The cure plate was now ready for the autoclave. Figure 
6.3 shows a schematic of the cure layup. 

6.2.2 Glass Fiber / Epoxv Wet Layup 

The glass fiber / epoxy system requires a wet layup procedure. The raw materials 
used in the preparation of this system were a 1 meter wide roll of woven unidirectional 
(the transverse weave is just used to keep the cloth together, it has no strength) glass cloth 
type 1543-38, a drum of Rutapox L20 epoxy resin and another smaller container of liquid 
SL hardener for the resin. The Columbian Conductex Ultra High Conductivity carbon 
black was added directly into the L20 epoxy resin by Penn Color Inc. Therefore the 
percent mass of carbon black is based on the mass of the resin only. The resin and 
hardener were mixed in the correct proportion. This mixture was then added to the dry 
fabric plies, producing a layup ready to be cured in the autoclave under vacuum. 

The wet layup was performed in the following manner. First, the necessary plies 
of glass fabric were cut from the roll with sharp scissors. This was accomplished by 
tracing out the plies to be cut using a felt pen and a specially designed Teflon coated 
aluminum pattern. While cutting the fabric, some shearing of the fabric occurred; this 
was straightened out when the plies were place on the curing plate, by carefully flattening 
out all the wrinkles in the fabric. The glass plies to be used in each laminate were 
carefully weighed. The appropriate mass of liquid matrix was prepared by mixing the 
resin and hardener in a 100 to 34 mass ratio to give a total mass of matrix equal to the 
total mass of the fabric plies. The layup calls for the combining of equal masses of 
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Figure 6.3: Schematic of AS4/3501-6 cure layup 
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matrix and fibers. For the carbon black filled resin, the total mass of filled resin used was 
increased to keep the mass of actual resin without carbon black the same, as shown in 
Table 6.5. Once mixed together the epoxy was well stirred to ensure complete mixing of 
the resin and hardener. 

In preparation for the wet layup and cure, the cure plate was sprayed with Mold 
Wiz® release agent and covered with a sheet of nonporous Teflon. On top on this, three 
layers of cork dam were used to make three 305 mm by 356 mm (12" x 14") rectangular 
vessels for the laminates, followed by another cork dam made of one layer of cork, this 
second dam was used to create a step down from the 3 layer dam height. A sheet of 
nonporous Teflon was placed in the bottom of the vessel, to line the bottom and sides of 
the vessel, followed by a sheet of peel-ply and then the layers of dry fabric. After the last 
ply of fabric was placed in the vessel, the matrix was poured on top and allowed to soak 
into the plies. The matrix was spread across the top ply using a half inch acid brush, 
starting at the center of the laminate and working outward radially toward the edges. By 
working from the center outward, the amount of trapped air was minimized. Care was 
taken not to disturb the alignment of the fibers while spreading the matrix with the brush. 

When wetted, the layup changed appearance from white in color to transparent. 
However, this only occurred with the pure epoxy as the epoxy with carbon black is black 
in color. With the higher percentages of carbon black, the epoxy did not soak through the 
plies very well. In these cases it was necessary to proceed to the next step without 
waiting for the epoxy to soak through the plies completely. When no areas of white were 
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Table 6.5: Rutapox L20/SL resin and hardener mixing ratios used 


Carbon Black 

Resin 

Resin Mass 

Hardener Mass 

0% 

100% 

92.16 g 

31.34 g 

5% 

95% 

97.01 g 

31.34 g 

10% 

90% 

102.40 g 

31.34 g 

15 % 

85 % 

108.42 g 

31.34 g 

20% 

80% 

11 5.67 g 

31.34 g 
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visible, any remaining matrix was spread evenly over the laminate, and a sheet of peel- 
ply was smoothed on top, followed by a sheet of porous Teflon. The aluminum top plate, 
sprayed with Frekote 700 and wrapped in nonporous Teflon, came next. Sheets of 
bleeder paper were rolled into dams and added around the cork vessels and the vacuum 
ports to absorb any extra liquid epoxy and to prevent it from flowing into the vacuum 
system. The assembly was then covered with a single sheet of porous Teflon, several 
layers of glass breather, and finally, the vacuum bag, which was sealed to the caul plate 
with vacuum tape. The cure plate was now ready for the autoclave. Figure 6.4 shows a 
schematic of the cure layup, and Figure 6.5 shows a typical cure plate with the bleeder 
paper rolls around the laminate vessels to absorb the excess epoxy. 

6.2 3 AS4/3501-6 Cure 

The AS4/3501-6 material was cured according to the TELAC net cure procedure, 
in the TELAC autoclave. The TELAC autoclave is a 0.91 m (3 ft) diameter Baron 
Blakeshee Inc. model BAC 35 autoclave, as shown in Figure 6.6, it has a 1.68 m (5.5 ft) 
length, a maximum temperature capability of 427 °C (800 °F), a maximum pressure 
capability of 1 .72 MPa (250 psi), and a maximum vacuum capability of 762 mm Hg (30 
in Hg). Six laminates were cured at one time in the autoclave. The position of the 
laminates on the cure plate is shown in Figure 6.7. 

A vacuum check was performed on the cure plate before it was rolled into the 
autoclave. This was done by applying a vacuum of 762 mm Hg (30 in Hg) and then 
shutting off the vacuum source. The vacuum seal was considered satisfactory if less than 
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Figure 6.4: Schematic of glass fiber / epoxy cure layup 
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Figure 6.5: Typical cure plate with bleeder paper rolls around laminate vessels 
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Figure 6.6: TELAC’s Baron Blackshee autoclave 
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1 27 mm Hg (5 in Hg) was lost in 5 minutes. Otherwise, attempts were made at trying to 
find and repair leaks in the vacuum bag, or the entire bag was replaced. Once a 
successful vacuum check was performed, the cure plate was rolled into the autoclave and 
the vacuum was rechecked. 

In the autoclave, a vacuum of 762 mm Hg (30 in Hg) was applied to the cure 
plate. The autoclave pressure was raised to a gage pressure of 0.59 MPa (85 psi) and 
held. Once the pressure was achieved, the autoclave temperature was raised at a rate of 1 
to 3 °C per minute until 1 16 °C (240 °F) was reached. This temperature was held for 1 
hour, following which the temperature was raised at the same rate to 177 °C (350 °F) and 
held for 2 hours. Finally the temperature was decreased at a rate 3 to 5 °C per minute to 
80 °C (176 °F) at which point the pressure was released. The temperature, pressure, and 
vacuum cycles are shown in Figure 6.8. The laminates were postcured in an oven at 177 
°C (350 °F) for 8 hours with no applied pressure or vacuum. 

6.2,4 Glass Fiber / Epoxy Cure 

The glass fiber / epoxy was cured according to the manufacture’s specification in 
the TELAC autoclave, refer to section 6.2.3 for description of the autoclave. Three 
laminates were cured at one time in the autoclave. The position of the laminates on the 
cure plate is shown in Figure 6.9. 

After a successful vacuum check, refer to section 6.2.3 for explanation, a vacuum 
of 762 mm Hg (30 in Hg) was applied to the cure plate in the autoclave and the 
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temperature was raised to 60 °C (140 °F) at a rate of 5 to 6 °C per minute. No external 
pressure was applied to the laminates. This temperature was held for 8 hours, after which 
the vacuum was released and the laminates were removed from the autoclave. The 
temperature, pressure, and vacuum cycles are shown in Figure 6.10. The laminates were 
postcured in an oven at 80 °C (176 °F) for 15 hours with no applied pressure or vacuum. 

6.2.5 Final Preparation 

All laminates were milled into coupons with a Van Norman model 22-L milling 
machine mounted with a water-cooled 220 grit diamond cutting wheel. The diameter of 
the cutting wheel is 254 mm and a spindle rate of 1 100 rpm was used. The milling table 
feed rate was 279 mm per minute. First, the edges of the laminate were cleaned by 
trimming approximately 3 mm (1/8 in) from each side. For the AS4/3501-6 samples, two 
25.4 mm (1 in) wide strips were cut off the shorter side, making sure to start with the 
“good” corner side, after which a 50.8 mm (2 in) wide strip was cut off. The strips were 
then cut into the required sized samples, three 50.8 mm by 50.8 mm, three 25.4 mm by 
50.8 mm, and six 25.4 mm by 25.4 mm. As shown in Figure 6.1 1, 9 conductivity and 3 
dielectric strength samples were cut. For the glass fiber / epoxy samples, three 50.8 mm 
(2 in) wide strips were cut off the shorter side, followed by a 25.4 mm (1 in) wide strip. 
The first 50.8 mm strip was not used, in an effort to eliminate curing defects located near 
the edge of the laminates. The next two were used for tensile tests. The 25.4 mm wide 
strip was cut into six 25.4 mm by 25.4 mm samples. The first 2 from the edge were 
discarded for edge defects, the next 2 samples were used for the conductivity 
measurement, and the remaining 2 were used for the dielectric strength measurement, as 
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Figure 6.1 1 : Schematic of the cutting patterns for the AS4/3501-6 laminates 




shown in Figure 6.12. The same cutting process was used to manufacture the 50.8 mm x 
76.2 mm (2" x 3") glass loading tabs, from sheets of purchased fiberglass. Once the glass 
tabs were cut, one edge was tapered to approximately 45° on a belt sander. 

6.3 CONDUCTIVITY MEASUREMENT 

ASTM D257-93 Standards summarizes the issues related to measurement of DC 
resistivity in insulating materials. 61 To measure the conductivity of the samples, the 
current had to be uniformly distributed across the surface of the samples. Therefore the 
application of a conductive coating was required. After experimenting with various 
coatings, including gold sputtering, vapor deposited silver, and conductive epoxy with 
various curing processes, it was decided to use conductive epoxy with an aluminum foil 
surface. To prevent moisture from penetrating all the samples, they were stored in an air- 
tight jar with desiccant. 

6.3.1 Conductive Epoxy Cure 

The conductive epoxy used was Epo-tek 410LV two part silver epoxy, 
manufactured by Epoxy Technology. The resin and the hardener were mixed at a specific 
ratio by mass of 15.0 resin to 2.1 hardener. The resin was thoroughly mixed before being 
used to ensure an even distribution of the silver. Once the two parts were thoroughly 
mixed together, it was applied to the surface of the conductivity samples with a half inch 
acid brush. The sample was then placed on a piece of aluminum foil, and a second piece 
of aluminum foil was place on top of the sample. A weight was placed on top of the 
samples to ensure an even distribution of the epoxy. The samples were then cured in an 
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Figure 6. 12: Schematic of the cutting patterns for the glass fiber / epoxy laminates 
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oven. The AS4/3501-6 samples were cured for 1 hour at 100 °C (212 °F), and because of 
the lower cure temperature of the glass fiber / epoxy, it was cured for 4 hours at 60 °C 
(140 °F). The edges of the samples were then sanded to remove any epoxy the may have 
ended up on the sides of the sample creating a conductive path between the two surfaces. 

6.3.2 Surface Laver Removal 

After the first set on conductivity measurements, it was postulated that the epoxy 
surface layer was affecting the measurement. Therefore the conductivity samples were 
sanded to remove the surface layer. The sanding process was done with a hand-held 
orbital sander, the sample were held to the table with double-sided tape. Once the 
sanding was complete, the sample thicknesses were re-measured and a new conductive 
epoxy layer was applied. 

6.3.3 Carbon Fiber / Epoxy Measurement 

The conductivity of the carbon fiber / epoxy samples was determined by 
measuring the resistance across the sample and calculating the conductivity base on 
equation 81. The resistance was measured using the “direct method” of measurement. 
This refers to the application of a known voltage, and measuring the resulting current. 
The instrument used was a Keithley 6517 High Resistance Electrometer with an applied 
voltage of 0.05 Volts. It was borrowed from the Adaptive Materials and Structures 
Laboratory (AMSL). The sample was placed between to electrodes in a sample holder 
made of Lexan, a highly resistive material, as shown in Figure 6.13. The sample holder 
was placed in a shielded enclosure, and a two probe technique specified by the 
manufacturer was used to reduce signal noise. Figure 6.14 shows the test setup used. 


210 




Figure 6.13: Schematic of the conductivity sample holder 
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Figure 6. 14: Experimental setup used for the conductivity measurements 
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Note the batteries and switch-box assembly are used for the glass fiber / epoxy 
measurements only. 

6.3.4 Glass Fiber / Epoxv Measurement 

The glass fiber / epoxy samples required a different procedure because of their 
high resistance. When direct measurements were tried with these samples, the conductive 
coatings on the surfaces made the sample act like a capacitor - the resistance reading kept 
increasing as the current able to flow decreased. To get the resistance of the sample, the 
sample was assumed to be modelable as a resistor and capacitor in parallel, and the RC 
circuit discharge constant was measured (refer to section 4.6.1). By using two 6 volt 
batteries in series, 12 volts was applied across the sample. The voltage source was 
removed, and 500 voltage and time data points were recorded by the Keithley 6517 
electrometer. The experimental setup can be seen in Figure 6.14. Once the data was 
collected a curve was fit to the data, thus giving the time constant from which the 
conductivity was calculated using equations 83 and 81. 

6.4 DIELECTRIC CONSTANT MEASUREMENT 

The dielectric constant was calculated by measuring the capacitance of the 
sample. The conductivity samples were used since they were already coated with 
conductive epoxy and they were not damaged by the conductivity measurements. The 
samples were placed in the above mentioned Lexan sample holder, and the capacitance 
was measured using a Continental Specialties Corporation model 3001 Capacitance 
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Meter borrowed from AMSL. Figure 6.15 shows the experimental setup used. The 
dielectric constant was calculated using equation 85. 

6.5 DIELECTRIC STRENGTH MEASUREMENT 

ASTM D 149-94 Standards summarizes the issues related to measurement of 
dielectric strength of insulating materials. 62 The testing procedure described here was 
developed by Aaron Bent of AMSL. 63 The testing was performed using a Philips model 
PM5138 Function Generator, which can output a DC source voltage from 0 to 1.0 in 0.01 
increments. The two voltage amplifiers used were a TREK model 663A H.V. power 
supply and model 662 ±10,000 volt amplifier and a Kepco Bipolar Operational Power 
Supply/ Amplifier model BOP 1000M, which supplies 1000 volts. This test setup was 
also borrowed from AMSL. When using the Kepco amplifier the voltage increment was 
10 volts, while when using the TREK amplifier the voltage increment was 100 volts. 
Figure 6.16 shows the experimental setup used. The samples were placed in a silicone oil 
test fixture, as shown in Figure 6.17. Silicone oil surrounded the specimen to prevent 
flashover and partial discharges, as per ASTM D 149-94. The oil was kept at room 
temperature. The sample was held between two 0.25" diameter hemispherical electrodes 
(electrode type #5 in ASTM D 149-94), which were connected to the high voltage 
amplifier. Hemispherical electrodes make contact with a discrete point of the sample, in 
contrast to electroded samples which would result in measurement of the “weakest link” 
dielectric strength of the entire sample. To prevent moisture absorption, the samples 
were stored in an air-tight jar with desiccant before testing. 


214 




215 


Figure 6.16: Experimental setup used to measure the dielectric strength 
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Figure 6. 1 7: Schematic of the dielectric strength sample holder 
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The Method B, Step-by-Step testing method of ASTM D 149-94 was used in the 
application of voltage. A DC voltage was applied to the specimen, in single increments 
of the voltage generator. At each time step, the voltage was held for 5 seconds soak time, 
and then immediately dialed to the next voltage. The starting voltage for each sample 
was approximately 50% of the breakdown voltage. For a test to be valid at least 5 
increments need to be made before breakdown. The recorded breakdown voltage is the 
highest level reached where the sample survived for the entire 5 second duration. 
Breakdown occurs when the voltage is sufficiently high to allow a current path through 
the sample material. The test was first done using the Kepco amplifier. If breakdown did 
not occur by the 1000 volt maximum of the amplifier, the TREK amplifier was then used. 
The dielectric strength was calculated using equation 86. 

6.6 DENSITY MEASUREMENT 

The density measurement was performed on the conductivity samples before the 
conductive epoxy surface layer was applied, since the area and thickness of the samples 
was already measured. The sample mass was then measured using an OHaus model 
TS4K0 digital balance, and the density was calculated using equation 88. 

6.7 TENSILE TESTS 

The tensile test was performed to determine the Young’s modulus, Poisson ratio, 
failure stress, and failure stain. The testing apparatus used was a 1 10,000 lbs MTS model 
31 1.21 testing machine with an Instron model 8500 Plus Controller, and a Macintosh IIx 
for data acquisition. Figure 6.18 shows the experimental setup used. The tests were 
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performed using stroke control with a rate of 0.02 in/sec, and a gripping pressure of 500 
psi. The load range used was ±10,000 lbs, and the stroke range used was ±1". The data 
acquisition was done using the “TELAC Data Acquisition” Labview program using a 
sampling frequency of 2 Hz. The procedure followed is documented in the “TELAC 
Manufacturing Course Class Notes”. 59 All the tensile test data can be found in Appendix 
C. 


Prior to testing, the samples were prepared by bonding tapered glass loading tabs 
to each end of the samples on both sides. The tabs were bonded using Epoxi-Patch 0151 
Clear, a 2 part epoxy which cures at room temperature in 48 hours with steel weights 
placed on the tabs to hold them in place. Following the tab application, 2 Measurements 
Group EA-06-125AD-120 strain gages were applied to the tensile specimen using a M- 
Bond 200 strain gage adhesive kit. More detailed directions are given in the “TELAC 
Manufacturing Course Class Notes”. 59 The gages were placed near the center of the 
sample, one in the longitudinal direction and the other in the transverse direction, as 
shown in Figure 6.19. 

6.7.1 Young’s Modulus Measurement 

The Young’s modulus was calculated by plotting the applied stress versus the 
longitudinal stain and graphically measuring the slope of the linear portion of the curve, 
as shown in Figure 6.20. The Young’s modulus was calculated using equation 89. 
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Figure 6.19: Illustration of the location of the strain gages on the tensile test sample 
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6.7.2 Poisson Ratio Measurement 


The Poisson ratio was calculated by plotting the transverse strain versus the 
longitudinal strain and graphically measuring the slope of the linear portion of the curve, 
as shown in Figure 6.21 . The Poisson ratio was calculated using equation 90. 

6.7.3 Failure Stress Measurement 

The failure stress was calculated from the maximum load applied to the sample. 
The load can be converted into stress using equation 91. 

6.7.4 Failure Strain Measurement 

The failure strain is the strain in the material at the maximum load. 
Unfortunately, by the time the sample failed the strain gages had long since stopped 
functioning. Therefore, the strain was calculated from the stroke using equation 92. This 
value can typically be inaccurate, as using stroke as a measurement of strain is a gross 
approximation. However, in some cases, as shown in Figure 6.20, there can be quite 
good agreement between the strain gage strain and the stroke strain. 


6.8 MICROSCOPY 

Microscopy was performed on the carbon fiber / epoxy samples and the glass 
fiber / epoxy samples to measure the thickness of the pure epoxy surface layer and to look 
for percolation paths. The carbon black particles were not visible under either the optical 
microscope or the scanning electron microscope. In preparation for the microscopy the 
samples were potted in Buehler EPO-Kwick two part epoxy, which cures at room 
temperature in 24 hours. Once the epoxy had cured, the samples were polished using 
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Figure 6.21 : Typical transverse versus longitudinal strain curve 



diamond grit polish on a Buehler Ecomet IV polisher with an Euromet I power head. The 
samples were examined using an Olympus model BH-2 optical microscope with a 
maximum magnification of 750 times. Micro-graphs were taken of the samples, from 
which the surface layer thickness was measured several times and averaged. 
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CHAPTER 7 


EXPERIMENTAL RESULTS AND 
DISCUSSION 


This section presents the results from the experimental portion of this research. 
The through-thickness electrical properties of carbon fiber / epoxy are discussed first. 
Then the development of a conductivity-tailorable glass fiber / epoxy composite is 
discussed, and the change in the electrical and mechanical properties due to adding 
conductive carbon black to the composite are presented. Some typical results will be 
presented here, while all the results can be found in Appendix C. 

7.1 CARBON FIBER / EPOXY ELECTRICAL PROPERTIES 

The data is shown in Figures 7.1 through 7.3, and Tables 7.1 through 7.3. The 
same formatting is used on all three figures. The sample electrical resistances are plotted 
versus sample area, and the sample breakdown voltages are plotted versus sample 
thickness. The mean of each group of data, (with each thickness and layup), is shown. 
The error bars represent one standard deviation. The unidirectional samples are displayed 
with open symbols, and the quasi-isotropic samples are displayed with filled-in symbols. 
The thickness of the samples is represented by increasing size of the symbol; a triangle 
for the 4 ply samples, a circle for the 8 ply samples, and a square for the 32 ply samples. 
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Figure 7. 1 : Unsanded carbon fiber / epoxy resistance data 
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Figure 7.2: Sanded carbon fiber / epoxy resistam 
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Figure 7.3: Unsanded carbon fiber / epoxy breakdown voltage data 






Table 7.1: 


Carbon fiber / epoxy resistance and conductivity thickness comparison 


Test Case 

Thickness 

Resistance (Q) 

Conductivity (1/Q-m) 


(mm) 

Average 

S.D. 

Average 

S.D. 

Unsanded 

0.52 

0.634 

0.279 

0.817 

0.273 


0.99 

0.620 

0.226 

1.476 

0.605 


3.75 

0.556 

0.351 

6.791 

2.485 

Sanded 

0.38 

0.305 

0.209 

1.587 

0.705 


0.84 

0.150 

0.073 

7.071 

3.722 


3.63 

0.216 

0.252 

33.727 

17.345 

Surface 

0.52 

0.305 

0.209 

0.015 

0.005 

Layer 

0.99 

0.150 

0.073 

0.014 

0.006 : 


3.75 

0.216 

0.252 

0.017 

0.007 


S.D.: Standard deviation 
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Table 7.2: 


Carbon fiber / epoxy resistance and conductivity area comparison 


Test Case 

Area 

(mm 2 ) 

Resista 

Average 

nee (Q) 
S.D. 

Conductiv 

Average 

ity (1/Q-m) 
S.D. 

Unsanded 

634 

0.875 

0.221 

3.723 

3.514 

1269 

0.557 

0.213 

3.372 

3.795 

2558 

0.378 

0.099 

1.990 

1.897 

Sanded 

634 

0.384 

0.261 

11.449 

12.472 

1269 

0.165 

0.065 

13.079 

17.580 

2558 

0.122 

0.087 

17.848 

23.104 

Surface 

Layer 

634 

0.875 

0.221 

0.020 

0.005 

1269 

0.557 

0.213 

0.016 

0.006 

2558 

0.378 

0.099 

0.011 

0.003 


S.D.: Standard deviation 


Table 7.3: Carbon fiber / epoxy resistance and conductivity layup comparison 


Test Case 

Layup 

Resista 

Average 

nee (Q) 
S.D. 

Conductiv 

Average 

ity (1/Q-m) 
S.D. 

Unsanded 

Unidirectional 

0.737 

0.286 

2.417 

2.392 

Quasi-isotropic 

0.469 

0.195 

3.639 

3.697 

Sanded 

Unidirectional 

0.304 

0.238 

10.611 

16.319 

Quasi-isotropic 

0.143 

0.090 

17.639 

18.614 

Surface 

Layer 

Unidirectional 

0.737 

0.286 

0.012 

0.004 

Quasi-isotropic 

0.469 

0.195 

0.019 

0.005 


S.D.: Standard deviation 
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For clarity, the horizontal placement of the data has been shifted slightly so the x-axis 
variation within each area or thickness group is not meaningful. The data in the tables is 
the average over all samples with the specified parameter (thickness, area, and stacking 
sequence or layup), and S.D. represents the standard deviation of all these samples. 

7.1.1 Through-Thickness Conductivity 

The raw data shows high scatter. As shown in Figure 7.1, there are no observable 
thickness trends. The resistance decreases with increasing area, but it is unclear if there is 
the expected (inverse) area dependency. Also, there seems to be a stacking sequence 
trend; the resistance of the quasi-isotropic samples are lower than that of the 
unidirectional samples. These resistance values are all low, and all close to the noise 
floor of the electrometer of approximately 0.1 Q. 

Figure 7.2 shows the resistance of the samples after the epoxy rich surface layer 
has been removed. Comparison of Figures 7.1 and 7.2 shows that the resistance of the 
samples decreased after sanding, as would be expected since the bulk composite 
conductivity is much greater than the pure epoxy conductivity (refer to Table 2.1). There 
is still a fair amount of scatter in the data. The sanded resistances are for the most part in 
or near the noise floor of the electrometer. The data still seems to have no apparent 
thickness trend. There seems to be the expected trend of decreasing resistance with 
increasing area. Also, the stacking sequence trend observed in the unsanded samples is 
less clear. 
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When the resistance data is interpreted as conductivity, there are non-classical 
thickness and surface area dependencies. If conductivity is a true material property, all 
data from the various samples with different thickness, areas, and stacking sequences 
would reduce (using equation 8 1 ) to the same value of conductivity. Instead, there seems 
to be some decrease in conductivity with area, and a distinct increase in conductivity with 
thickness. The range of values of conductivity fit within the range, 0.1 to 106 1/Q-m, of 
through-thickness conductivities for carbon fiber / epoxy composites previously reported 
(refer to Table 2.1). When the resistance of the sanded samples is interpreted as 
conductivity, there seems to be no area dependency, but there is still a non-classical 
thickness dependency. Most of the sanded conductivity values are close to the upper 
previously reported value of 106 1/Q-m. 

As noted in section 4.6.1, it is possible that the resistance of the sample is 
dominated by low resistance surface layers. The apparent thickness of these were 
measured, and the unsanded resistance data was interpreted as the conductivity of the 
surface layer. The apparent epoxy surface layer thickness was measured using 
microscopy on samples from the same laminate as the conductivity samples. The 
thickness of the surface layer was measured 9 times on each surface and averaged. For 
the unidirectional laminates it was approximately 9.0 pm, and for the quasi-isotropic 
laminates it was approximately 10.0 pm. Results are shown in Tables 7.1 through 7.3. A 
reasonably consistent set of conductivities results. There is a slight non-classical area 
dependency, but no thickness dependency. The values of the conductivity computed in 
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this way are above the range of reasonable values for the conductivity of epoxy by 
approximately one order of magnitude, indicating that the surface layer is not a uniform 
insulator. 

It can be seen that the measurement of the through-thickness conductivity of 
carbon fiber / epoxy composites is difficult to achieve accurately. The measured values 
of the conductivity are within the range of previously reported values. However, there 
are non-classical area and thickness dependencies, and a high degree of scatter. The 
through-thickness conductivity appears to be influenced by both the layup and the 
presence of an epoxy surface layer. The conductivity does not appear to behave as a 
continuum property; it is more likely dominated by percolation effects, in which the 
properties of the lowest resistance percolation path are being measured, which will vary 
from sample to sample. 

7.1.2 Through-Thickness Dielectric Strength 

The through-thickness dielectric strength was measured using samples with no 
applied conductive surface layer. When the raw breakdown voltage was looked at there 
was no observable trend with thickness. The data seemed to be scattered around a value 
of approximately 1500 volts, as shown in Figure 7.3. It was postulated that the epoxy 
surface layer might again be dominating the data, so the surface layers were sanded off as 
done above with the conductivity samples. When the sanded samples were tested, no 
breakdown voltage could be measured as there was a current flow at the lowest possible 
voltage increment. 
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When the breakdown voltage data is interpreted as a material property (dielectric 
strength), there is a non-classical thickness dependency. Instead of being a constant value 
with some scatter, there seems to be a decrease in dielectric strength with thickness. The 
range of data however is within the previously reported range for carbon fiber / epoxy 
composites. When the breakdown voltage data is interpreted as the dielectric strength of 
the epoxy surface layer (using equation 87), more consistent results are obtained. The 
dielectric strength of the surface layer measured in this way shows a fair bit of scatter 
around a constant value of approximately 150 MV/m, which is within the range of 
previously reported values for organic polymers. 

The through-thickness dielectric strength of carbon fiber / epoxy composites 
seems to be epoxy surface layer dominated. The measured breakdown voltages make 
sense if interpreted as the dielectric strength of the epoxy surface layer. As with the 
conductivity, it may not be meaningful to assign a continuum value to the dielectric 
strength of carbon fiber / epoxy composites. 

7.1.3 Summary 

The through-thickness electrical properties, the conductivity and the dielectric 
strength, of carbon fiber / epoxy composites are dominated by the epoxy surface layer. 
This effect is observed in carbon fiber / epoxy composites because of the large difference 
in the conductivity of the bulk composite and the surface layer. Even when this layer is 
sanded off, there is still a high degree of scatter and non-classical thickness dependencies 
in the resistance data, making the determination of a bulk conductivity difficult. 
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7.2 CONDUCTIVITY TAILORABLE GLASS FIBER / EPOXY 


The results of tests on the conductivity-tailorable glass fiber / epoxy composites 
were as desired — the electrical properties were increased significantly while the 
mechanical properties displayed minimal change. Figures 7.4 through 7.10 represent the 
effects of adding conductive carbon black on the material properties of the composite. 
The actual values can be found in Table 7.4, where the listed property is the average of 
the six samples measured, and S.D. is the standard deviation of these six samples. 

Some manufacturing difficulties were encountered during the cure process. The 
main reason for these was that the higher percentages of carbon black made the epoxy 
very viscous. This resulted in the epoxy not flowing well during the cure process, which 
resulted in thicker laminates, as shown in Figure 7.11. This increase in thickness was 
created when the extra epoxy on the top surface of the laminate did not flow out around 
the top cure plate during the curing process. The surface layer was measured using 
microscopy and the results can be seen in Figure 7.12, and Table 7.5. The total surface 
layer thickness of the carbon black filled laminates was greater than that of the pure glass 
fiber / epoxy laminate. In particular, the surface layer thickness for the 15% and 20% 
carbon black laminates increased significantly. Micrographs of the samples are shown in 
Figure 7.13. The black circles represent voids in the potting epoxy. 

A side effect of this surface layer was to cause the laminates to bend when they 
cooled from the cure temperature to room temperature, as shown in Figure 7.14. The 
reason for this effect was that the surface layer caused the laminate to appear to be 
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Young's Modulus (GPa) 



Figure 7.7: Glass fiber / epoxy Young’s modulus 
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Poisson Ratio 



Figure 7.8: Glass fiber / epoxy Poisson ratio 
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Figure 7.9: Glass fiber / epoxy failure stress 
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Figure 7. 1 0: Glass fiber / epoxy failure strain 
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Table 7.4: Conductivity-tailorable glass fiber / epoxy composite results 


Property 


Percentage Carbon Black By Mass 



0% 

5% 

10% 

15% 

20% 

Conductivity (1/Q-m) 

8.82 x 10' 13 

4.37 x 10' 12 

2.86 x 1 0' 12 

4.40 x 1 O'" 

1.23 x 10‘ 9 

S.D. 

2.45 x 10’ 13 

1.45 x 1 O' 12 

2.54 x 10' 12 

2.19 x 10'" 

1.49 x 10' 9 

Dielectric Constant 

5.30 

6.33 

9.20 

29.81 

82.23 

S.D. 

0.12 

0.16 

0.36 

2.08 

54.53 

Density (kg/m 3 ) 

1607 

1618 

1618 

1565 

1523 

S.D. 

18 

30 

12 

16 

10 

Young's Modulus (GPa) 

1 1.92 

10.47 

10.40 

10.43 

10.47 

S.D. 

0.49 

0.37 

0.25 

0.23 

0.48 

Poisson Ratio 

0.572 

0.552 

0.543 

0.589 

0.569 

S.D. 

0.045 

0.017 

0.028 

0.029 

0.030 

Failure Stress (MPa) 

105.8 

114.9 


110.4 

103.1 

S.D. 

9.1 

6.5 


5.4 

2.7 

Failure Strain (ps/e) 

59590 

63940 

53200 

49100 

49250 

S.D. 

7890 

10270 

6190 

9120 

3660 


S.D.: Standard deviation 
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Figure 7.11: Glass fiber / epoxy thickness 
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Figure 7. 12: Glass fiber / epoxy surface layer thickness 
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Table 7.5: 


Glass fiber / epoxy surface layer thickness 


Surface Layer 


Percentage Carbon Black By Mass 


Thickness (mm) 

0% 

5% 

10% 

15% 

20% 

Upper Surface 

0.075 

0.131 

0.129 

0.224 

0.486 

S.D. 

0.029 

0.053 

0.082 

0.027 

0.046 

Lower Surface 

0.065 

0.097 

0.102 

0.079 

0.087 

S.D. 

0.041 

0.053 

0.057 

0.046 

0.030 


S.D.: Standard deviation 


247 












0 % Carbon Black 


5% Carbon Black 





1 0 % Carbon Black 


15% Carbon Black 



20% Carbon Black 


Figure 7.13: Micro-graphs of the surface layer of the conductivity-tailorable glass fiber 

/ epoxy composite 


248 










unsymmetric; therefore thermal bending strains were induced in the laminate as it cooled. 
This curvature had no effect on the mechanical properties since only a small load was 
required to straighten out the laminate. 

7.2.1 Conductivity 

The conductivity, as shown in Figure 7.4, increased with the addition of carbon 
black to the epoxy resin. It increased by up to 3 orders of magnitude with the 20% 
carbon black samples. The conductivity seems to be constant until approximately 10% 
carbon black, at which point there is a rapid increase. This is consistent with percolation 
theory, where there is no significant increase in properties until there are sufficient 
particles to create percolation paths from one edge of the sample to the other, as shown in 
Figure 4.14. Therefore it would appear that the percolation limit for this carbon black 
system is around 1 0% carbon black by mass of epoxy resin. 

7.2.2 Dielectric Constant 

The dielectric constant, as shown in Figure 7.5, also displayed the same trend as 
the conductivity, remaining fairly constant up to 10% carbon black then increasing 
rapidly. Exponential increases were observed by Bent 64 , who examined up to 5% carbon 
black samples, and by Yacubowicz 65 , who observed dielectric constants in the hundreds 
with a lower percolation limit. Both the large variation in the dielectric constant at 20% 
carbon black and the very high values may be due to unusual percolation paths, as shown 
in Figure 7.15. These paths would have the same effect as reducing the apparent 
thickness of the sample used in the calculation of the dielectric constant from the 
measured capacitance (refer to equation 85). If this geometry existed, the dielectric 
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Figure 7.15: Schematic of unusual percolation paths which can lead to high dielectric 

constant 
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constant of the material (K actua ,) should be calculated using the effective thickness (h eff ). 
However, the dielectric constant (K measured ) is actually calculated using the sample 
thickness (h); this results in the measured dielectric constant being magnified by the ratio 
of the sample thickness over the effective thickness (h/h cfT ). Therefore, very big effective 
values for the dielectric constant of the sample are possible. This variation also shows the 
randomness of carbon aggregation that can radically affect the electrical properties, even 
with extremely careful dispersion during manufacturing of the samples, and careful 
measurement techniques. 

7.2.3 Dielectric Strength 

When the dielectric strength was measured using samples with no added 
conductive surface layer, no data was obtainable as none of the samples had broken down 
when the maximum voltage, 1 0,000 volts, was applied across the samples. Therefore, the 
dielectric strength of all samples was greater than 10.5 MV/m. This result makes sense 
since the published values for the dielectric strength of glass fiber / epoxy composites is 
between 17.7 and 21.7 MV/m. 

7.2.4 Density 

Theoretically the density should increase linearly with increasing percentage of 
carbon black. This trend is not seen in Figure 7.6, probably due to the high viscousness 
of the epoxy, which limited its ability to soak completely through the glass fiber bundles, 
resulting in voids. Also as mentioned earlier, the increase in surface layer thickness is 
also affecting the density measurements. 
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7.2.5 Tensile Tests 


Tensile tests were performed on the laminates to determine the Young’s modulus, 
Poisson ratio, failure stress, and failure strain. To verify that epoxy was still usable, an 
initial set of tests were performed on unidirectional glass fiber / epoxy laminates. The 
results are presented in Table 7.6, and they agreed with data measured by Ed Wolf when 
he used the same system. It should be noted that all the failures initiated near one of ends 
close to the loading tabs, as shown in Figures 7.16 through 7.21. This indicates that the 
data is not a very good absolute strength measurement. It is used here for comparison 
with the carbon black filled data. The data was also verified using MCLAM, an in-house 
CLPT code using a Tsai-Wu failure criterion, with typical glass fiber / epoxy material 
properties. The results are shown in Table 7.7. There was excellent agreement for the 
Young’s modulus, Poisson ratio, and failure stress of the [±45] s laminates. There was 
also good agreement for the Young’s modulus and Poisson ratio of the [0] 4 laminates. 
The failure stress measured was 50% of the calculated failure stress, due to the splitting 
failure mode initiated near the loading tabs (refer to Figure 7.16). 

The Young’s modulus, as shown in Figure 7.7, decreases slightly with increasing 
percentage of carbon black. The Poisson ratio, as shown in Figure 7.8, remained fairly 
constant around 0.56 with a little scatter above and below this value. The failure stress, 
as shown in Figure 7.9, remains fairly constant, slightly increasing with increasing 
percentage of carbon black then slightly decreasing. The failure strain, as shown in 
Figure 7.10, decreases linearly with increasing percentage of carbon black. Even though 
there is a significant decrease in the failure strain from just over 6% strain to 5% strain, 
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Table 7.6: [0] 4 glass fiber / epoxy tensile test results 


Property 

Average 

S.D. 

Young's Modulus (GPa) 

49.62 

1.23 

Poisson Ratio 

0.216 

0.018 

Failure Stress (MPa) 

554.4 

23.9 

Failure Strain (ps/s) 

22134 

1091 


S.D.: Standard deviation 
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Figure 7.19: Broken 10% carbon black tensile specimens 
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Figure 7.21 : Broken 20% carbon black tensile specimens 
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Table 7.7: 


MCLAM results for [0] 4 and [±45] s glass fiber / epoxy laminates 


Property 

[0] 4 

[±45] s 

Young's Modulus (GPa) 

38.59 

12.55 

Poisson Ratio 

0.260 

0.516 

Failure Stress (MPa) 

1061.79 

99.91 
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these values are well above any acceptable amount of strain that should be seen in an 
actual aerospace structure. It should also be noted that the high values of failure strain 
was calculated from the stroke data since at failure the strain gages had long since 
stopped working. This value can typically be inaccurate, as using stroke as a 
measurement of strain is a gross approximation, as with the [0] s laminates (refer to Figure 
C. 1 ). However, in some cases there can be quite good agreement between the strain gage 
strain and the stroke strain, as with the [±45] s laminates which had acceptable to excellent 
agreement (refer to Figure C.29). 

7.2.6 Summary 

With the addition of conductive carbon black into the insulative matrix of a 
composite, the conductivity can be increased significantly, as was desired. The dielectric 
constant was also increased. Due to testing limitations no trends could be found for the 
dielectric strength, however the value is greater than 10.5 MV/m in all cases. The 
increase in electrical properties had little effect on the mechanical properties of the 
composite. They remained fairly constant across all values of carbon black added. The 
only problem with the addition of carbon black to the matrix of the composite was the 
manufacturing difficulties caused by the increase in the epoxy viscosity. This problem 
could be overcome by using a less viscous epoxy, or by developing a manufacturing 
technique which would allow the epoxy to flow out of the composite during the curing 
process. Possible manufacturing technique variations include heating the epoxy before 
mixing so that it is less viscous, using pressure during the cure to press the top cure plate 
down into the laminate and thus force the extra epoxy out, or modifying the cure cycle to 
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include a hold time at the flow temperature of the epoxy (as is done with AS4/3501-6 
cure cycle) instead of ramping straight up to the cure temperature. 
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CHAPTER 8 


CONCLUSIONS 


A tool now exists which can determine the charge density, electric field, and 
voltage as a function of depth into a composite material, in a specified orbital 
environment, and thus predict if an electrostatic discharge will occur. The program, 
CoDDCA, can be run for any orbit and material system, by specifying the relevant 
parameters. The program can also include geomagnetic substorms and solar particle 
events of various intensities specified by the user. The space environment is determined 
using the Environmental Workbench software by specifying the orbit parameters (the 
apogee height, perigee height, inclination and solar cycle condition). 

For conducting materials, such as carbon fiber / epoxy, there is no apparent 
problem with deep dielectric charging, however, peculiarities of the composite such as 
surface layers can be a concern. For insulating materials, such as glass fiber / epoxy, 
polymers, and Kevlar, there is no electrostatic discharge concern with average day-to-day 
flux levels. Problems can arise during geomagnetic substorms and solar particle events, 
where flux levels increase by several orders of magnitude. There is also a concern for 
electrostatic discharges if materials which have extremely low conductivities are used. 
The conductivity enhancements, such as high field conductivity and radiation induced 
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conductivity, have a minimal effect on the deep dielectric charging analysis and in any 
case only help to reduce the likelihood of an electrostatic discharge. The 
micromechanical details between the fiber and matrix on the electric field are in the most 
part negligible, with a maximum electric field magnification of 3 times. Unlikely fiber 
arrangements in percolation paths can lead to electric field increases on the order of 50 
times. Another related problem is that of the electric field in the insulative surface layer 
of a conductive composite. The results show that there can be a large electric field in the 
surface layer even if the bulk composite is not charged significantly, due to the large 
difference in conductivities between these two areas. 

Parametric studies showed several relationships exist between the deep dielectric 
charging parameters and the material properties. The time to steady state is proportional 
to the dielectric constant and inversely proportional to the conductivity. The electric field 
is inversely proportional to the conductivity and proportional to the amount of particles 
present in the space environment. The charge density, electric field, and voltage 
distributions are affected by the location of particles in the material and the relative 
amounts of electrons and protons. 

From the case study orbits, the worst orbit was the base case orbit, a 7000 km 
altitude circular orbit. A similar environment could be expected in a geosynchronous 
transfer orbit. Next came the geosynchronous orbit, the Global Positioning System orbit, 
and the Molniya orbit. These all produced an order of magnitude smaller electric fields in 
the base case material, as compared to the base case orbit. The space shuttle parking orbit 
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and the low altitude retrograde or polar orbit produced negligible electric fields due to the 
small numbers of high energy particles at low altitudes. 

A complete set of tests were done on the through-thickness electrical properties of 
a carbon fiber / epoxy composite. The results were incompatible with the idea that the 
electrical properties were continuum material properties on the scale of the specimen. 
The actual numerical values were within the range of previously measured values (0.1 to 
106 1/Q-m for the through-thickness conductivity), however the thickness, area, and 
stacking sequence dependencies were either non-classical or indiscernible. The through- 
thickness electrical properties appear to be affected by the surface layer, since the 
conductivity increased and the dielectric strength goes away when the surface layer is 
sanded off. The other factor affecting the through-thickness electrical properties is 
percolation effects. Therefore, when dealing with electrical properties of composites one 
has to be wary of previously published numbers, and attention must be paid to surface 
effects. 


A conductivity tailorable composite can be manufactured using a glass fiber / 
epoxy composite with conductive carbon black added to the epoxy resin. The 
conductivity can be increased by 3 orders of magnitude (from 10 12 1/Q-m to 10' 9 1/Q-m) 
with minimal change to the mechanical properties of the composite. The only drawback 
to using high percentages of carbon black in the epoxy resin is the increased 
manufacturing difficulties due to the high viscosity of the epoxy. This tailorable system 
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may be very useful if a material used in the space environment needs to be an insulator, 
but deep dielectric charging also needs to be minimized. 
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APPENDIX A 


FORTRAN SOURCE CODE AND USAGE 

INTRODUCTION 

CoDDCA is an easy-to-use computer code which will calculate the charge 
density, electric field, and voltage distributions through the thickness of the material, as a 
function of time. All the equations used in this code are given in Chapter 4. The code 
can solve the problem in closed form or using a central difference finite difference time 
stepping routine. 


INPUTS 

The inputs required are the material properties requested by the code and entered 
by the user (use D to represent a power of 10, hence 1.0D5 = l.Ox 10 5 ), the electric field 
or voltage boundary conditions, and the environment data. The environmental data is 
calculated by the Environmental Workbench software. There are two files, one for the 
electrons and one for the protons. The first column is the energy of the particles in MeV, 
and the second column is the particle fluence in #/cm 2 for 1 year. There should be 100 
data entries. 
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ENVIRONMENTAL WORKBENCH DATA 

This section explains how to get the required environment data from EWB. Once 
EWB is running, to define the desired orbit select ‘Orbit’ from the ‘System Definition’ 
menu; this will open the orbit definition window. Enter the Apogee height, Perigee 
height, and inclination of the desired orbit, and set the mission duration to 365 days. 
When done press the ‘Update Form’ button. To close the window press the ‘Dismiss’ 
button. Next open the trapped electron and proton fluences windows. This is done by 
selecting ‘Trapped Electron Fluence’ and ‘Trapped Proton Fluence’ from the 
‘Environments’ menu. The only thing that need to be done in these windows is to set the 
solar cycle condition to either ‘1’ for solar maximum or ‘2’ for solar minimum, followed 
by pressing the ‘Update Form button. To generate the actual data files, the XY Plot 
generator windows needs to be opened. This is accomplished by selecting XY Plots 
from the ‘Param-Studies’ menu. For simplicity the instructions will be given only for 
generating the electron fluence data file; the same procedure is followed for the protons. 
To generate the data, select one of the ‘Max Energy Range’ cells in the ‘Trapped Electron 
Fluence’ window, then press the ‘Select’ button in the ‘Primary Independent Variable’ 
section of the ‘XY Plot Specification’ window. This sets the first column of the output to 
be the electron energy. In the same section of the XY Plot Specification window, set 
the minimum value to 0.01, the maximum value to 100, the number of points to 100, and 
select a logarithmic spacing. The dependent variable needs to be set; this is done by 
selecting the ‘Fluence’ cell to the right of the ‘Max Energy Range’ cell selected earlier, in 
the ‘Trapped Electron Fluence’ window. Now press the ‘Select’ button in the ‘Dependent 
Variable’ section of the ‘XY Plot Specification’ window to set the second column of the 
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output as the electron fluence. Press the ‘Calculate’ button in the ‘XY Plot Specification’ 
window to calculate the data. The data can be viewed by pressing the ‘Display’ button in 
the ‘XY Plots’ window. To save the data press the ‘File Table’ button in the ‘XY Plots’ 
window, and enter a name for the file when asked. 


SAMPLE INPUT 

The following is an example of an EWB output file. The output file is in Arial 
font and comments on the file are to the right of the file and in “Times New Roman” font 
and contained within quotation marks. For convenience the data has been truncated by 
removing the middle rows. 


Energy Column 

Fluence Column 


[MeV] 

[#/cm2] 


En_Max_(MEV) 

Ele_Fluence_ 

Row 1 , Column heading 

9.9999998E-03 

O.OOOOOOOE+OO 

Row 2 

1.0974987E-02 

0.0000000E+00 

Row 3 

1.2045036E-02 

0.0000000E+00 

Row 4 

1.321 941 2E-02 

0.0000000E+00 

Row 5 


7.5646576E+01 

8.8328767E+15 

Row 98 

8.3022034E+01 

8.8328767E+15 

Row 99 

9.1 1 16585E+01 

8.8328767E+15 

Row 100 

1.0000034E+02 

8.8328767E+15 

Row 101 
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SAMPLE SESSION 


The following section is a sample session using the base case data inputs. Everything 
displayed on the screen are in Arial font, the user inputs are bold, and comments added 
for explanation are in “Time New Roman” font and contain within quotation marks. 


Select Model Type (enter no.): 

1: For Closed Form Solution 
2: For Finite Difference Solution 

1 

Select type of boundary conditions (enter no.) 


Type Boundary Condition 

1 E(0) V(0) 

2 E(0) V(L) 

3 E(L) V(0) 

4 E(L) V(L) 

5 V(0) V(L) 

4 

Enter value of boundary condition 1 (V/m) 

0D0 

Enter value of boundary condition 2 (V) 

0D0 

Enter EWB electron data filename 

e-07k-max 


E stands for electric field 
V stands for voltage 
0 stands for at the front surface 
L stands for at the back surface 


D sets number to double precision 


Filenames up to 25 characters are accepted 


Enter EWB proton data filename 

p-07k-max 

Enter the material dark conductivity (1/ohm-m) 

ID-10 

Enter the dielectric constant 

3D0 

Enter the dielectric strength (V/m) 

15D6 

Enter the material density (kg/m3) 

1600D0 
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Enter the material thickness (m) 

2.5D-3 

Include radiation induced conductivity (y/n)? Upper or lower case is acceptable 

n 

Include substorms (y/n)? 

n 

Include solar particle events (y/n)? 

n 

Calculating 

OUTPUTS 

The code outputs 2 or 3 ASCII files depending on the user options selected. Two 
files are always created: ‘Incoming Environment.txt’ which contains space environment 
data, and ‘CoDDCA Output.txt’ which contains the charge density, electric field, and 
voltage distributions. The third file, ‘Conductivity.txt’, is only created if high field 
conductivity is selected. It contains the conductivity and high field conductivity scaling 
factor distributions. The distributions are output as a set of tables, the first column is the 
depth and the following columns are the distributions at various times (one time per 
column). The times are in the first row of the table. The data the columns are separated 
by tabs. 

SAMPLE OUTPUT 

The section contain examples of the output files. As above the output file is in 
Arial font and comments are to the right and in Times New Roman font, enclosed by 
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quotation marks. The tabs between columns is represented using ». Also for clarity only 
a portion of the output files are present, this is accomplished by omitting the middle rows. 


The first file is ‘Incoming Environment.txt’. The first 3 rows are used for 
pertinent data used in the run. The first column of the first table is the energy bins 
selected by EWB, the second column in the depth an electron penetrates the material for 
the given energy in column 1, the third column is the electron flux, which is simply the 
EWB fluence divided by 365 days (3.153 x 10 7 seconds). The fourth and fifth columns 
are the same as the second and third columns but for protons. Following this table of data 
a second table is displayed containing the incoming charge density rate data. The first 
column is the depth, the second column in the incoming charge density rate due to the 
electrons, the third column is the incoming charge density rate due to the protons, and the 
last column is the total incoming charge density rate. 


‘Incoming Environment.txt’ 


Electron data file e 

Proton data file P 

Density (kg/m3) 1600.000000000000000 

Energy (MeV) E Depth (m) E Flux (#/m2-s) P Depth (m) P Flux (#/m2-s) 

99999998E-02 .1725361 IE-04 .OOOOOOOOE+OO 20449374E-06 .00000000E+00 

.1 0974987 E-01 20653364E-04 OOOOOOOOE+OO 21585172E-06 .00000000E+00 

12045036E-01 24469699E-04 OOOOOOOOE+OO 22839133E-06 .OOOOOOOOE+OO 

.13219412E-01 .28734697E-04 OOOOOOOOE+OO .24222959E-06 .00000000E+00 


Row 1 
Row 2 
Row 3 
Row 4 
Row 5 
Row 6 
Row 7 
Row 8 
Row 9 


75646576E+02 10000000E+04 .OOOOOOOOE+OO .11244411E+00 .83095465E+06 Row 102 
'83022034E+02 10000000E+04 OOOOOOOOE+OO 12894112E+00 .68189026E+06 Row 103 
.91116585E+02 10000000E+04 .00000000E+00 .14808055E+00 .48842372E+06 Row 104 
.10000034E+03 .10000000E+04 .OOOOOOOOE+OO .17040497E+00 .55502696E+06 Row 105 
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Depth (m) ChargelnE (#/m3-s) 
.OOOOOOOOE+OO .OOOOOOOOE+OO 
.50000000E-05 .OOOOOOOOE+OO 
.1 OOOOOOO E-04 .OOOOOOOOE+OO 
.15000000E-04 .OOOOOOOOE+OO 


ChargelnP (#/m3-s) Chargeln (#/m3-s) 
59039949E+1 1 .1 1807990E+17 
.78446859E+16 .78446859E+16 
56096617E+16 .56096617E+16 
.45846771 E+16 .45846771 E+16 


Row 106 
Row 107 
Row 108 
Row 109 
Row 1 1 0 
Row 1 1 1 


24850000E-02 - 62344377E+13 
24900000E-02 -.61302806E+13 
24950000E-02 -.60261234E+13 
.25000000E-02 -.59219663E+13 
End of file 


.22244680E+13 -.40099697E+13 
22225795E+13 -.39077011E+13 
.2220691 0E+1 3 -.38054324E+13 
.22188025E+13 -.37031638E+13 


Row 605 
Row 606 
Row 607 
Row 608 
Row 609 


The second file is ‘CoDDCA Output. txf. The first 25 rows are used for pertinent 
data used in the run, the actual number of rows used depends on the options selected. 
The following three tables are of the same format, the first one contains the charge 
density distribution, the second table contains the electric field distribution, and the last 
table contains the voltage distribution. The first column in the depth, and the number of 
columns used for the distribution is always 1 0 for the closed form solution, but it can vary 
for the finite difference solution. 


‘CoDDCA Output.txt’ 


Closed Form Solution 


Row 1 

Conductivity (1/ohm-m) 

1.000000000000000E-010 

Row 2 

Dielectric Constant 

3.000000000000000 

Row 3 

Dielectric Strength (V/m) 

1 500000000000000E+007 

Row 4 

Thickness (m) 

2.500000000000000E-003 

Row 5 

Density (kg/m3) 

1600.000000000000000 

Row 6 

Delta x (m) 

5.000000000000000E-006 

Row 7 

Surface Thickness (m) 

5.000000000000000E-006 

Row 8 

Electron data file 

e-07k-max 

Row 9 

Proton data file 

p-07k-max 

Row 10 

Time Constant (sec) 

2.656256339999999E-001 

Row 1 1 

Steady-state time (sec) 

1.223251250330985 

Row 12 
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Rad. Conductivity 


F 


Other information goes here depending on options selected 


Steady-state reached 

Column 1 Column 2 Column 3 .... Column 1 1 

Charge Density (#/m3) at times (sec) indicated at the top of the column 
Depth (m) 12232513E+00 24465025E+00 • .12232513E+01 

.00000000E+00 .1 1575041E+16 .18878397E+16 ■ ■ ■ ■ .31051397E+16 

50000000E-05 76899251E+15 12541940E+16 • • • ■ .20629122E+16 

10000000E-04 .54989937E+15 89686242E+15 • ■ • ■ .14751693E+16 


Row 13 
Row 14 
Row 15 
Row 16 
Row 17 
Row 18 
Row 19 
Row 20 
Row 21 
Row 22 
Row 23 
Row 24 
Row 25 
Row 26 
Row 27 
Row 28 
Row 29 
Row 30 
Row 31 
Row 32 
Row 33 


.24850000E-02 -39308606E+12 - 64110659E+12 ■ -.10544992E+13 

24900000E-02 - 38306095E+12 -.62475607E+12 • ■ • -10276057E+13 

24950000E-02 -37303585E+12 -60840555E+12 • • -.10007 122E+ 13 

25000000E-02 -.36301074E+12 -59205503E+12 • • • -97381869E+12 

Electric Field (V/m) at times (sec) indicated at the top of the column 
Depth (m) .12232513E+00 .24465025E+00 - .12232513E+01 

.00000000E+00 .14909883E+04 24317382E+04 39997500E+04 

50000000E-05 .15200387E+04 .24791 182E+04 .4077681 2E+04 

10000000E-04 .15399268E+04 .251 15549E+04 • .41310333E+04 

.15000000E-04 .15549959E+04 .25361321E+04 ■ ■ ■ • .41714581E+04 


Row 528 
Row 529 
Row 530 
Row 531 
Row 532 
Row 533 
Row 534 
Row 535 
Row 536 
Row 537 
Row 538 


24850000E-02 .34204457E-01 .5578601 0E-01 • ■ • ■ .91757447E-01 

.24900000E-02 22500626E-01 .36697561 E-01 • • •• .60360555E-01 

.24950000E-02 .1 1099140E-01 18102225E-01 • ■ ■ • .29774740E-01 

.25000000E-02 OOOOOOOOE+OO 00000000E+00 - ■ • • .OOOOOOOOE+OO 

Voltage (V) at times (sec) indicated at the top of the column 
Depth (m) 12232513E+00 24465025E+00 • ■ ■ • .12232513E+01 

.00000000E+00 .1 1 141258E+01 . 181 7091 7E+01 • • • • .29887725E+01 

.50000000E-05 .11065982E+01 .1 80481 45E+01 ■ • • ■ .29685790E+01 

.10000000E-04 .10989483E+01 .17923379E+01 • • • • .29480572E+01 

15000000E-04 .109121 10E+01 .1 77971 86E+01 ■ • • ■ .29273009E+01 


Row 1032 
Row 1033 
Row 1034 
Row 1035 
Row 1036 
Row 1037 
Row 1038 
Row 1039 
Row 1040 
Row 1041 
Row 1042 


281 



.24850000E-02 
24900000E-02 
.24950000E-02 
.25000000E-02 
End of file 


.25350997E-06 
.1 1 174727E-06 
.27747850E-07 
.00000000E+00 


41346395E-06 .68007009E-06 

18225503E-06 ■ • • ■ 29977509E-06 

.45255561 E-07 • • • 74436848E-07 

.00000000E+00 • ■ • • .00000000E+00 


Row 1536 
Row 1537 
Row 1538 
Row 1539 
Row 1540 


The third file is ‘Conductivity.txt’. The first 4 rows are used for pertinent data 
used in the run, the actual number of rows used depends on the options selected. The 
following two tables are of the same format, the first one contains the high field 


conductivity scaling factor distribution, and the second table contains the conductivity 
distribution. The first column is the depth, and the number of column used for the 


distributions can vary since this file can only be created when the finite difference 


solution is used. 


‘Conductivity.txt’ 


Conductivity (1/ohm-m) 
Dielectric Constant 
Temperature (degree C) 


1 .000000000000000E-01 0 

3.000000000000000 

120.000000000000000 

Column 3 • • • 

24400000E+00 • • • 
10001694E+01 • ■ • 
10001727E+01 • • • 
10001750E+01 • • 

.10001767E+01 ■ ■ 


Column 8 

Row 1 
Row 2 
Row 3 
Row 4 
Row 5 

.79300000E+00 

Row 6 
Row 7 

.10002675E+01 

Row 8 

.10002727E+01 

Row 9 

.10002762E+01 

Row 10 

.10002789E+01 

Row 1 1 


Column 1 Column 2 

High Field Conductivity Scaling Factor 


Depth (m) 

.00000000E+00 

50000000E-05 

.10000000E-04 

.15000000E-04 


.12200000E+00 

.10001039E+01 

.10001059E+01 

.10001073E+01 

.10001083E+01 


.24850000E-02 

.24900000E-02 

24950000E-02 

.25000000E-02 


.10000000E+01 

.10000000E+01 

.10000000E+01 

.10000000E+01 


.10000000E+01 • 
.10000000E+01 • 
.10000000E+01 ■ 
10000000E+01 ■ 


10000000E+01 Row 505 
.10000000E+01 Row 506 
.10000000E+01 Row 507 
.10000000E+01 Row 508 
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Conductivity (1/ohm-m) 


79300000E+00 

Row 509 
Row 5 1 0 

Depth (m) .12200000E+00 

24400000E+00 

Row 51 1 

00000000E+00 10001039E-09 

.10001694E-09 • ■ • 

10002675E-09 

Row 5 1 2 

.50000000E-05 10001059E-09 

.10001727E-09 • • • 

10002727E-09 

Row 5 1 3 

10000000E-04 .10001073E-09 

.10001750E-09 • • 

10002762E-09 

Row 5 1 4 

.15000000E-04 .10001083E-09 

10001767E-09 • • 

• 10002789E-09 

Row 5 1 5 

24850000E-02 .10000000E-09 

10000000E-09 • • 

. 1 0OOOOOOE-09 

Row 1009 

24900000E-02 10000000E-09 

.10000000E-09 ■ • 

• 10000000E-09 

Row 1010 

.24950000E-02 .10000000E-09 

.10000000E-09 ■ • 

• 10000000E-09 

Row 1011 

25000000E-02 .10000000E-09 

End of file 

.10000000E-09 

10000000E-09 

Row 1012 
Row 1013 


SOURCE CODE 

! Program CoDDCA, Composite Deep Dielectric Charging Analyzer 
! © 1997 Steven A. Czepiela 
! Massachusetts Institute of Technology 
! 77 Massachusetts Avenue 
! Cambridge MA 02139 

i Permission to use, copy, and modify this software and its documentation for internal purposes 
i only without fee is hereby granted provided that the above copyright notice and this permission 
I appear on all copies of the code and supporting documentation. For any use of this software, in 
! original or modified form, including but not limited to, adaptation as the basis of a commercial 
! software or hardware product, or distribution in whole or part, specific prior permission and/or the 
! appropriate license must be obtained from MIT. This software is provided “as is” without any 
! warranties whatsoever, either expressed or implied, including but not limited to the implied 
! warranties of merchantibility and fitness for a particular purpose. This software is a research 
! program and MIT does not represent that it is free of error or bugs or suitable for any particular 

! task. 

I This program calculates the electric field in a composite material, caused by radiation 
! using either a closed form solution for the time or a finite difference solution. 

! The initial data comes from the Environmental Workbench software (EWB), and 
! is read by the program, the first column is energy [MeV], and the second column is 
! fluence [#/cm2], this data is for 1 year. Radiation induced conductivity and high 
I field conductivity may be included. As well as substorms and solar particle events, 

! which are modeled by using a scaling factor to increase the amount of incoming 
! fluence. The program outputs the electric field, voltage, and charge density every x 
! seconds, for an approximate output of 10 different times, in the file called 'CoDDCA 
! Output.txt'. The program also outputs the incoming particle data in the file called 
! 'incoming Environment.txt'. If high field conductivity is selected, a third output 
! file is created containing the conductivity distribution for the same time intervals 
! as mentioned above, in the file called 'Conductivity.txt'. The code terminates when 
i the electric field exceeds the dielectric strength of the material, or when the 
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! electric field reaches steady state. The surface (Surface) is considered to be 1/500 
! of the thickness. 


! Variable Definitions 


! Again 

Logical variable if iterating again on finite difference routine 

! Bad 

Logical variable if invalid user input 

! BC1 

Boundary condition 1 

! BC2 

Boundary condition 2 

! BCType 

Boundary Condition option selection 

! Beta 

Frenkel parameter [C-N-m2] A 0.5 

I Break 

Logical variable if breakdown has occurred 

I BreakTime 

Time breakdown occurred [s] 

! Cl 

First integration constant 

!C2 

Second integration constant 

! Charge1(501) 

Charge density 1 through material [#/m3] 

! Charge2(501) 

Charge density 2 through material [#/m3] 

! ChargeE(IOO) 

Incoming electron charge density rate [#e/m3-s] 

! Chargeln(501) 

Net incoming particle charge density rate [#/m3-s] 

! ChargelnE(501) 

Incoming electron charge density rate [#e/m3-s] 

! ChargelnP(501) 

Incoming proton charge density rate [#p/m3-s] 

! ChargeP(IOO) 

Incoming proton charge density rate [#p/m3-s] 

I Charges(501,20) 

All charge densities through material [#/m3] 

! Conduc 

Conductivity of the material [1/ohm-m] 

! Count 

Number of iterations in finite difference solution 

! Count2 

Number of iterations in FD solution since the last data output 

! C e 

Charge of one electron [C] 

! C_K 

Conversion from Celsius to Kelvin 

! DarkCon 

Dark conductivity of material [1/ohm-m] 

1 Delta 

Radiation induced conductivity exponent 

I DeltaX 

Space step [m] 

! DeltaT 

Time step [s] 

! Density 

Density of material [kg/m3] 

! Depth(501) 

Depth into material [m] 

! DepthE(IOO) 

Depth electrons penetrate into material [m] 

! DepthP(IOO) 

Depth protons penetrate into material [m] 

! Dielec 

Permittivity of the material [C2/N-m2] 

! DieStr 

Dielectric Strength of the material [V/m] 

! Dose 

Net dose rate [rad/s] 

! DoseE 

Electron dose rate [rad/s] 

! DoseP 

Proton dose rate [rad/s] 

! Energy(IOO) 

Energy of particles from EWB [MeV] 

! Error 

Value below which the residual is small enough, ie steady-state 

! FDCon(501) 

Conductivity [1/ohm-m] 

! FDCons(501,20) 

All conductivities [1/omh-m] 

! Field(501) 

Electric field through material [V/m] 

! Fields(501,20) 

All electric fields through material [V/m] 

! FileE 

EWB electron fluence filename 

! FileP 

EWB proton fluence filename 

! FluenceE(IOO) 

Electron fluence from EWB [#e/cm2] 

! FluenceP(IOO) 

Proton fluence from EWB [#p/cm2] 

! FluxE(IOO) 

Electron flux [#e/m2-s] 

! FluxP(IOO) 

Proton flux [#p/m2-s] 
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! HFC 

! HFCon(501) 

! HFCons(501,20) 

! I, J 

! Integral1(501) 

! Integral2(501) 

! J_MeV 
! K 

! KRIC 
! LETE(IOO) 

! LETP(IOO) 

! MDie 
! N1 
! N2 

! Ocount 
! Option 
! Permittivity 
! PCount 
! Pie 
! QE 
! QP 

! RadCon 
! Rad_Jkg 
! Residual 
! RIC 
! sec yr 
! SPE 

! SPEFactor 
! SS 

! SSTime 
! SSFactor 
! Surface 
! Temp 
! Thickness 
! Time 
! Times(20) 

! TimeConst 
! TimeFactor 
! Tmpl 
! Tmp2 
! TmpC 
! Voltage(501) 

! Voltages(501,20) 


Logical variable if high field conductivity is included 
High field conductivity scaling factor 
All high field conductivity scaling factors 
Do loop counters 

First integration of charge density [V/m] 

Second integration of charge density [V] 

Conversion factor from joules to mega electron volt [J/MeV] 
Boltzmann's constant [J/K] 

Coefficient of radiation induced conductivity [s/ohm-m-rad] 
Electron LET [MeV-m2/kg] 

Proton LET [MeV-m2/kg] 

Dielectric constant of material 
Number of EWB data points 

Number of points used in the space direction (thickness direction) 

Counter for the number of outputs 

Closed form or finite difference option selection 

Permittivity of free space [C2/N-m2] 

Number of iterations to be performed before data send to file 
Pie, 3.1415.... 

Electron quality factor 
Proton quality factor 

Radiation induced conductivity [1/ohm-m] 

Conversion factor from Rad to J/kg [Rad-kg/J] 

Finite difference residual 

Logical variable if radiation induced conductivity is included 
Conversion factor from seconds to a year [s/yr] 

Logical variable if solar particle events are included 

Scaling factor for solar particle events 

Logical variable if substorms are included 

Time to steady-state 

Scaling factor for substorms 

Thickness of surface layer [m] 

Temperature [K] 

Thickness of material [m] 

Time [s] 

All output times [s] 

Charging time constant [1/s] 

Scaling factor based on time in closed form solution 

Temporary Variable 1 

Temporary Variable 2 

Character temporary variable 

Voltage distribution through material [V] 

All voltage distributions through material [V] 


I The Main Program 


Implicit Logical (A-Z) 

Integer N1, N2, 1, J, BCType, Option, Count, Count2, PCount, OCount 
Double Precision DeltaX, SSTime, Conduc, Dielec, C_e, MDie, Density, DieStr, Tmpl, & 
sec_yr, Tmp2, Energy(IOO), FluenceE(IOO), FluenceP(IOO), FluxE(IOO), & 
FluxP(IOO), DepthE(IOO), DepthP(IOO), ChargeE(IOO), ChargeP(IOO), & 
Chargeln(501), Thickness, Charge1(501), Depth(501), Field(501), Permittivity, & 
TimeConst, Time, TimeFactor, J_MeV, Rad_Jkg, QE, QP, DoseE, DoseP, Dose, & 
LETE(IOO), LETP(IOO), RadCon, DarkCon, Delta, KRIC, Surface, Voltage(501), & 
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BC1, BC2, Integral-1 (501), Integral2(501), Cl, C2, BreakTime, ChargelnE(501), & 
ChargelnP(501), Temp, k, Beta, C_K, Pie, DeltaT, Residual, Error, HFCon(501) & 
FDCon(501), Charge2(501), SSFactor, SPEFactor, Charges(501,20), & 

Fields(501 ,20), Voltages(501,20), Times(20), HFCons(501,20), FDCons(501 20) 
Character FileE *25, FileP *25, TmpC *1 


! Set Constants 
N1 = 100 
N2 = 500 

Permittivity = 8.8541878D-12 
QP= 10 
QE = 1 

Pie = 3.14159266D0 
Error = 0.1 DO 
k = 1.380658D-23 
C_e = 1.6021 9D-1 9 
sec_yr = 3.153D7 
J_MeV = 1.602D-13 
Rad_Jkg = 100 
cm2_m2 = 1 0000 
C_K = 273.15D0 


1 Size of input files for electron & proton fluences 
1 Number of depth steps, size of output files 
Permittivity of free space (C2/N-m2] 

Quality factor for protons, used in Dose calculation 
Quality factor for electrons, used in Dose calculation 
Pie 

Value at which the residual is small enough 

Boltzmann's constant (J/KJ 

Number of coulombs per electron [C/#] 

Number of seconds per year [sec/yr] 

Number of joules per MeV [J/MeV] 

Number of rads per J/kg [rad-kg/J] 

Number of cm2 per m2 [cm2/m2] 

Conversion Factor from Celcius to Kelvin 


! Set I/O formats 

15 Format (E15.8, A1, E15.8, A1, E15.8, A1, E15.8, A1, E15 8) 
25 Format (E15.8, A1, E15.8, A1, E15.8, A1, E15.8) 

35 Format ('Depth (m)', 20(A1, El 5.8)) 

45 Format (E15.8, 20(A1, E15.8)) 


! Get user input 
Bad = .True. 

Do While (Bad) 

Print *, 'Select Model Type (enter no.):' 

Print *, ' 1: For Closed Form Solution' 

Print *, ’ 2: For Finite Difference Solution' 

Read *, Option 
If (Option <= 2) Then 
Bad = .False. 

Else 
Print * 

Print *, Invalid type entered, enter number between 1 and 2' 
Print * 

End If 
End Do 


Print * 

Bad = .True. 

Do While (Bad) 

Print *, 'Select type of boundary conditions (enter no.)' 
Print *, 'Type Boundary Condition' 

Print*,' 1 E(0) V(0)' 

Print *, ' 2 E(0) V(L)' 

Print *, ' 3 E(L) V(0)' 

Print *, ' 4 E(L) V(L)' 

Print *, ' 5 V(0) V(L)' 
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Read *, BCType 
If (BCType <= 5) Then 
Bad = False. 

Else 

Print * 

Print *, 'Invalid type entered, enter number between 1 and 5’ 
Print * 

End If 
End Do 

Print * 

If (BCType == 5) Then 

Print *, 'Enter value of boundary condition 1 (V)' 

Else 

Print *, 'Enter value of boundary condition 1 (V/m)' 

End If 

Read *, BC1 
Print * 

Print *, 'Enter value of boundary condition 2 (V)’ 

Read *, BC2 

Print * 

Print *, ’Enter EWB electron data filename’ 

Read *, FileE 
Print * 

Print*, ’Enter EWB proton data filename’ 

Read *, FileP 
Print * 

Print *, 'Enter the material dark conductivity (1/ohm-m)' 

Read *, DarkCon 
Print * 

Print *, 'Enter the dielectric constant' 

Read *, MDie 
Print * 

Print *, 'Enter the dielectric strength (V/m)’ 

Read *, DieStr 
Print * 

Print *, 'Enter the material density (kg/m3)' 

Read *, Density 
Print * 

Print *, 'Enter the material thickness (m)' 

Read *, Thickness 

Print * 

Print *, 'Include radiation induced conductivity (y/n)?' 

Read *, TmpC 

If ((TmpC == 'y') .OR. (TmpC == 'Y')) Then 
RIC = .True. 

Print * 

Print *, 'Enter coef. of radiation induced conductivity (s/ohm-m-rad)' 
Read *, KRIC 
Print * 

Print *, 'Enter the value for delta' 

Read *, Delta 
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Else 

RIC = False. 

End If 

If (Option == 2) Then 
Print * 

Print *, 'Include high field conductivity (y/n)?' 
Read *, TmpC 

If ((TmpC == 'y') .OR. (TmpC == 'Y')) Then 
HFC = True. 

Print * 

Print *, 'Enter temperature (degree C)' 

Read *, Temp 
Else 

HFC = False. 

End If 
Else 

HFC = .False. 

End If 

Print * 

Print *, 'Include substorms (y/n)?' 

Read *, TmpC 

If ((TmpC == 'y') OR. (TmpC == 'Y')) Then 
SS = .True. 

Print *, 'Enter substorm scaling factor:' 

Read *, SSFactor 
Else 

SS = .False. 

End If 
Print * 

Print *, 'Include Solar particle events (y/n)?' 

Read *, TmpC 

If ((TmpC == 'y') OR. (TmpC == 'Y')) Then 
SPE = True. 

Print *, 'Enter Solar particle event scaling factor:' 
Read *, SPEFactor 
Else 

SPE = .False. 

End If 

! Initialize variables 
DeltaX = Thickness / N2 
Surface = Thickness / 500 
Dielec = MDie * Permittivity 
Beta = (C_e**3 / (Pie * Dielec))**0.5 
Do 1=1, N2+1 

Depth(l) = (1-1)* DeltaX 
End Do 

! Open files 

Open (Unit = 1, File = FileE) 

Open (Unit = 2, File = FileP) 

Open (Unit = 3, File = 'Incoming Environment.txt') 


![m] 

! Surface layer thickness [m] 

! Dielectric of Material [C2/N-m2] 

! Frenkel Parameter [C-N-m2] A 0. 5 

![m] 
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Open (Unit = 4, File = 'CoDDCA Output.txt') 

If (HFC) Then 

Open (Unit = 5, File = 'Conductivity.txt') 

End If 

! Read data 
Read (1,*) 

Read (2,*) 

Do 1=1, N1 

Read (1,*) Energy(l), FluenceE(l) 

Read (2,*) Energy(l), FluenceP(l) 

End Do 

! Close files 
Close (Unit = 1) 

Close (Unit = 2) 

! Include substorms and solar particle events 
If (SS) Then 
Do 1=1, N1 

FluenceE(l) = FluenceE(l) * SSFactor 
End Do 
End If 

If(SPE) Then 
Do 1=1, N1 

FluenceP(l) = FluenceP(l) * SPEFactor 
End Do 
End If 

! Calculate fluxes 
FluxE(l) = 0 
FluxP(l) = 0 
Do 1=1, N1-1 

FluxE(l+1) = (FluenceE(l+1) - FluenceE(l)) * cm2_m2 / sec_yr ! [#/m2-s] 

FluxP(l+1 ) = (FluenceP(l+1) - FluenceP(l)) * cm2_m2 / sec_yr ! [#/m2-s] 

End Do 

! Calculate depth particles penetrate material 
Do 1=1, N1 

Tmpl = DLog10(Energy(l)) 

If (Energy(l) <= 7) Then ! The electron LET equation is only valid below 7 MeV 
Tmp2 = 0.0439D0 * Tmp1**6 + 0.1513D0 * Tmp1**5 + 0.0861D0 * Tmp1‘*4 - & 
0.1666D0 * Tmpl **3 + 0.0801 DO * Tmpl **2 + 0.4055D0 * Tmpl -0.6288D0 

LETE(I) = 10 ** Tmp2 ! [MeV-m2/kg] 

Depth E(l) = Energy(l) / (LETE(I) * Density) ! [m] 

Else 

LETE(I) = 0 
DepthE(l) = 1000 
End IF 

Tmp2 = - 0.004D0 * Tmp1**6 - 0.01 04D0 * Tmp1**5 + 0.0344D0 *Tmp1**4 + & 
0.1376D0 * Tmp1**3 - 0.2025D0 * Tmp1**2 - 0.8642D0 * Tmpl + 1.0404D0 
LETP(I) = 10** Tmp2 ! [MeV-m2/kg] 

DepthP(l) = Energy(l) / (LETP(I) * Density) ! [m] 

End Do 
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! Calculate dose rate, radiation induced conductivity, net material conductivity 
If (RIC) Then 
DoseE = 0 
DoseP = 0 
Do 1=1, N1 

DoseE = DoseE + FluxE(l) * LETE(I) * QE * J_MeV * Rad_Jkg ! [rad/s] 

DoseP = DoseP + FluxP(l) * LETP(I) * QP * J_MeV ‘ Rad_Jkg ! [rad/s] 

End Do 

Dose = DoseE + DoseP 
RadCon = KRIC * Dose**Delta 
Conduc = DarkCon + RadCon 
Else 

Conduc = DarkCon 
End If 

! Initialize more variables 
TimeConst = Dielec / Conduc 
SSTime = -1 * TimeConst * DLog(O.OIDO) 

Temp = Temp + C_K 

! Get more user input 
If (Option == 2) Then 
Print * 

Print *, The estimated time to steady-state is in seconds: ', SSTime 
Print *, 'Enter the time increment (seconds): ' 

Read *, DeltaT 

PCount = Nint(SSTime / (10 * DeltaT)) 

End If 

! Set-up output files 

Write (3,*) 'Electron data file', Char(9), Char(9), FileE 
Write (3,*) 'Proton data file', Char(9), Char(9), FileP 
Write (3,*) 'Density (kg/m3)', Char(9), Char(9), Density 
Write (3,*) 

Write (3,*) 'Energy (MeV)', Char(9), 'E Depth (m)', Char(9), 'E Flux (#/m2-s)', Char(9), & 

■P Depth (m)', Char(9), 'P Flux (#/m2-s)’ 

Do 1=1, N1 

Write (3,15) Energy(l), Char(9), DepthE(l), Char(9), FluxE(l), Char(9), DepthP(l), Char(9), & 
FluxP(l) 

End Do 

If (Option == 1 ) Then 

Write (4,*) 'Closed Form Solution’ 

Else 

Write (4,*) 'Finite Difference Solution' 

End If 

Write (4,*) 'Conductivity (1/ohm-m)', Char(9), Char(9), Char(9), Char(9), DarkCon 
Write (4,*) 'Dielectric Constant', Char(9), Char(9), Char(9), Char(9), MDie 
Write (4,*) 'Dielectric Strength (V/m)’, Char(9), Char(9), Char(9), Char(9), DieStr 
Write (4,*) 'Thickness (m)’, Char(9), Char(9), Char(9), Char(9), Thickness 
Write (4,*) 'Density (kg/m3)', Char(9), Char(9), Char(9), Char(9), Density 
Write (4,*) 'Delta x (m)\ Char(9), Char(9), Char(9), Char(9), DeltaX 
Write (4,*) 'Surface Thickness (m)', Char(9), Char(9), Char(9), Char(9), Surface 
Write (4,*) 'Electron data file', Char(9), Char(9), Char(9), Char(9), FileE 


! [rad/s] 

! [1/ohm-m] 
I [1/ohm-mj 


! [sec] 
! [sec] 
! [K] 
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Write (4,*) 'Proton data file', Char(9), Char(9), Char(9), Char(9), FileP 
Write (4,*) 'Time Constant (sec)’, Char(9), Char(9), Char(9), Char(9), TimeConst 
Write (4,*) 'Steady-state time (sec)', Char(9), Char(9), Char(9), Char(9), SSTime 
Write (4,*) 'Rad. Conductivity', Char(9), Char(9), Char(9), Char(9), RIC 
If (RIC) Then 

Write (4,*) 'Delta', Char(9), Char(9), Char(9), Char(9), Delta 

Write (4,*) 'Coef. Rad. Conductivity (s/ohm-m-rad)', Char(9), Char(9), Char(9), Char(9), KRIC 
Write (4,*) 'Electron dose rate (rad/s)', Char(9), Char(9), Char(9), Char(9), DoseE 
Write (4,*) 'Proton dose rate (rad/s)', Char(9), Char(9), Char(9), Char(9), DoseP 
Write (4,*) 'Total dose rate (rad/s)', Char(9), Char(9), Char(9), Char(9), Dose 
Write (4,*) 'Radiation induced conductivity (1/ohm-m)’, Char(9), Char(9), Char(9), Char(9), & 
RadCon 

End If 

If (Option == 2) Then 

Write (4,*) Time step (sec)', Char(9), Char(9), Char(9), Char(9), DeltaT 
Write (4,*) ’High Field Conductivity 1 , Char(9), Char(9), Char(9), Char(9), HFC 
If(HFC) Then 

Write (4 *) 'Temperature (degree C)', Char(9), Char(9), Char(9), Char(9), Temp - C_K 
Write (4^*) 'Frenkel Parameter (C-N-m2) A 0.5\ Char(9), Char(9), Char(9), Char(9), Beta 
End If 
End If 

If (SS) Then 

Write (4,*) 'Substorm scaling factor', Char(9), Char(9), Char(9), Char(9), SSFactor 
End If 

If (SPE) Then 

Write (4,*) 'Solar Particle Event scaling factor', Char(9), Char(9), Char(9), Char(9), SPEFactor 
End If 

If ( NOT. RIC) Then 
Write (4,*) 

Write (4,*) 

Write (4,*) 

Write (4,*) 

Write (4,*) 

Write (4,*) 

End If 

If (Option /= 2) Then 
Write (4,*) 

Write (4,*) 

End If 

If (.NOT. HFC) Then 
Write (4,*) 

Write (4T) 

End If 

If ( NOT. SS) Then 
Write (4,*) 

End If 

If ( NOT. SPE) Then 
Write (4,*) 

End If 

If (HFC) Then 

Write (5,*) 'Conductivity (1/ohm-m)’, Char(9), Char(9), Char(9), DarkCon 
Write (5T) 'Dielectric Constant', Char(9), Char(9), Char(9), MDie 
Write (5,*) 'Temperature (degree C)', Char(9), Char(9), Char(9), Temp - C_K 
If (RIC) Then 
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Write (5,*) 'Radiation induced conductivity (1/ohm-m)', Char(9), Char(9), Char(9), & 
Char(9), RadCon 

Else 

Write (5,*) 

End If 
End If 

! Calculate charge density 

ChargeE(l) = -1 * FluxE(l) / (DepthE(l) + (DepthE(2) - DepthE(l)) / 2) 

ChargeP(1 ) = FluxP(1 ) / (DepthP(1 ) + (DepthP(2) - DepthP(1 )) / 2) 

Do 1=2, N1-1 

If (DepthE(l) == 1000) Then 
ChargeE(l) = 0 

Else If (DepthE(l+1) == 1000) Then 

ChargeE(l) = -1 * FluxE(l) / (DepthE(l) - Depth(l-I)) 

Else 

ChargeE(l) = -1 * FluxE(l) / ((DepthE(l+1) - DepthE(l-l)) /2) 

End If 

ChargeP(l) = FluxP(l) / ((DepthP(l+1) - DepthP(l-l)) / 2) 

End Do 

If (DepthE(NI) == 1000) Then 
ChargeE(NI) = 0 
Else 

ChargeE(NI) = -1 * FluxE(NI) / ((DepthE(NI) - DepthE(NI-l)) / 2) 

End If 

ChargeP(NI) = FluxP(NI) / ((DepthP(NI) - DepthP(NI-l)) / 2) 

! Calculate surface charge density 
Tmpl = 0 
Tmp2 = 0 
Do 1=1, N1 

If (DepthE(l) < Surface) Then 

Tmpl = Tmpl + FluxE(l) ! [#/m2-s] 

End If 

If (DepthP(l) < Surface) Then 

Tmp2 = Tmp2 + FluxP(l) ! [#/m2-s] 

End If 
End Do 

Tmpl = -1 * Tmpl 
ChargelnE(l) = Tmpl 
ChargelnP(l) = Tmp2 

Chargeln(l) = (Tmpl + Tmp2) / Surface ! [#/m3-s] 

! Bin charge densities 
Do 1=2, N2+1 
Tmpl = 0.0 
Do J=1, N1-1 

If ((DepthE(J) < Depth(l)) And. (DepthE(J+1) > Depth(l))) Then 

Call Interp (Depth(l), DepthE(J), DepthE(J+1), ChargeE(J), ChargeE(J+1), Tmpl) 
Else If (DepthE(J) == Depth(l)) Then 

Tmpl = ChargeE(J) ! [#/m3-s] 

End If 
End Do 
Tmp2 = 0 


! [#/m3-s] 
! [#/m3-s] 

! [#/m3-s] 
! [#/m3-s] 
! [#/m3-s] 

! [#/m3-s] 
! [#/m3-s] 
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Do J=1, N1-1 

If ((DepthP(J) < Depth(l)) And. (DepthP(J+1) > Depth(l))) Then 

Call Interp (Depth(l), DepthP(J), DepthP(J+1), ChargeP(J), ChargeP(J+1), Tmp2) 
Else If (DepthP(J) == Depth(l)) Then 

Tmp2 = ChargeP(J) ! [#/m3-s] 

End If 
End Do 

ChargelnE(l) = Tmpl 
ChargelnP(l) = Tmp2 

Chargeln(l) = Tmpl + Tmp2 ! [#/m3-s] 

End Do 

! Output incoming charge density 
Write (3,*) 

Write (3,*) 'Depth (m)', Char(9), 'ChargelnE (#/m3-s)', Char(9), 'ChargelnP (#/m3-s)\ Char(9), & 
'Chargeln (#/m3-s)' 

Do 1=1, N2+1 

Write (3,25) Depth(l), Char(9), ChargelnE(l), Char(9), ChargelnP(l), Char(9), Chargeln(l) 
End Do 

! Perform Closed Form Solution 
If (Option == 1 ) Then 
OCount = 10 
Print * 

Print *, 'Calculating' 

! Calculate 1st and 2nd integral of incoming charge density 
Integrall(l) = 0 
Integral2(1) = 0 
Do 1=2, N2+1 

Integrall (I) = Integrall (1-1 ) + (Chargeln(l) + Chargeln(l-I)) 12 * DeltaX * C_e / Conduc 
I [V/m] 

Integral2(l) = Integral2(l-1) - (Integrall (I) + Integrall (1-1)) / 2 * DeltaX I [V] 

End Do 

! Calculate Charge, Field, Voltage for times equal to 10% through 100% steady-state 
Break = False. 

Do J=1, 10 

! Calculate charge density 

Time = J / 10.0 * SSTime ! [sec] 

TimeFactor = 1 - DExp(-1 * Time / TimeConst) 

Do 1=1, N2+1 

Chargel(l) = Chargeln(l) * TimeConst * TimeFactor ! [#/m3] 

End Do 

! Calculate constants of integration based on boundary conditions 
If ((BCType == 1 ) OR. (BCType == 2)) Then 
Cl = BC1 

Else If ((BCType == 3) .OR. (BCType == 4)) Then 
Cl = BC1 - TimeFactor* Integrall (N2+1) 

End IF 

If ((BCType == 1) .OR. (BCType == 3)) Then 
C2 = BC2 
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Else If ((BCType == 2) .OR. (BCType == 4)) Then 

C2 = BC2 - TimeFactor * lntegral2(N2+1) + Cl * Depth(N2+1) 

End If 

If (BCType == 5) Then 
C2 = BC1 

Cl = (TimeFactor * lntegral2(N2+1) + C2 - BC2) / Depth(N2+1) 

End If 

! Calculate voltage and electric field 
Do 1=1, N2+1 

Field(l) = TimeFactor * Integrall (I) + Cl ! [V/m] 

Voltage(l) = TimeFactor * Integral2(l) - Cl * Depth(l) + C2 ! [V] 

If (DAbs(Field(l)) >= DieStr) Then 
Break = .True. 

BreakTime = Time 
End If 
End Do 

! Store data 
Times(J) = Time 
Do 1=1, N2+1 

Charges(I.J) = Chargel(l) 

Fields(l,J) = Field(l) 

Voltages(I.J) = Voltage(l) 

End Do 
End Do 

! Perform Finite Difference Solution 
Else If (Option == 2) Then 

! First charge density is the input charge density 
Do 1=1, N2+1 

Chargel(l) = Chargeln(l) * DeltaT ! [#/m3] 

End Do 

! Iterate until steady-state is reached or breakdown occurs 

Count = 1 

Count2 = 1 

OCount = 1 

Break = .False. 

Again = True. 

Do While (Again) 

! Calculate 1st and 2nd integral of charge density 
Integrall (1) = 0 
Integral2(1) = 0 
Do 1=2, N2+1 

Integrall (I) = Integrall (1-1) + (Chargel(l) + Chargel(l-I)) 12 * DeltaX * C_e / Dielec 
! [V/m] 

Integral2(l) = Integral2(l-1) - (Integrall (I) + Integrall (1-1 )) 12 * DeltaX ! [V] 

End Do 

! Calculate constants of integration based on boundary conditions 
If ((BCType == 1 ) .OR. (BCType == 2)) Then 
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Cl = BC1 

Else If ((BCType == 3) .OR. (BCType == 4)) Then 
Cl = BC1 - lntegral1(N2+1) 

End IF 

If ((BCType == 1) OR. (BCType == 3)) Then 
C2 = BC2 

Else If ((BCType == 2) .OR. (BCType == 4)) Then 
C2 = BC2 - lntegral2(N2+1) + Cl * Depth(N2+1) 

End If 

If (BCType == 5) Then 
C2 = BC1 

Cl = (lntegral2(N2+1) + C2 - BC2) / Depth(N2+1) 

End If 

! Calculate voltage and electric field 
Do 1=1 , N2+1 

Field(l) = Integrall (I) + Cl 
Voltage(l) = Integral2(l) - Cl * Depth(l) + C2 
If (DAbs(Field(l)J >= DieStr) Then 
Break = True. 

Again = False. 

End If 
End Do 

! Calculate high field conductivity 
If (HFC) Then 
Do 1=1, N2+1 

HFCon(l) = (2 + DCosh(Beta * DSqrt(DABS(Field(l)))/ (2 * k * Temp))) / 3 
If (RIC) Then 

FDCon(l) = DarkCon * HFCon(l) + RadCon 
Else 

FDCon(l) = DarkCon * HFCon(l) 

End If 
End Do 
End If 

! Calculate charge density 
If (HFC) Then 

Charge2(1) = Chargel(l) + (Chargeln(l) - (FDCon(l) / Dielec * & 

Chargel (1)) - (Field(l) * ((FDCon(2) - FDCon(l)) / (DeltaX)))) * DeltaT 
Charge2(N2+1) = Charge1(N2+1) + (Chargeln(N2+1) - (FDCon(N2+1) / Dielec * & 
Chargel (N2+1)) - (Field(N2+1) * ((FDCon(N2+1) - FDCon(N2)) / (DeltaX)))) * & 
DeltaT 
Do 1=2, N2 

Charge2(l) = Chargel (I) + (Chargeln(l) - (FDCon(l) / Dielec * Chargel (I)) - & 
(Field(l) * ((FDCon(l+1) - FDCon(l-l)) / (2 * DeltaX)))) * DeltaT 

End Do 
Else 

Do 1=1, N2+1 

Charge2(l) = Chargel (I) + (Chargeln(l) - (Conduc / Dielec * Chargel (I))) * DeltaT 
End Do 
End If 

! Calculate Residual 
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Tmpl = 0 
If (HFC) Then 

If (Chargel(l) == 0) Then 

Tmpl = Tmpl + DAbs(Chargeln(1) - (FDCon(l) / Dielec * Chargel(l)) & 

- (Field(l) * ((FDCon(2) - FDCon(l)) / (DeltaX)))) 

Else 

Tmpl =Tmp1 + DAbs((Chargeln(1) - (FDCon(l) / Dielec * Chargel(l)) - & 
(Field(l) * ((FDCon(2) - FDCon(l)) / (DeltaX)))) * DeltaT/ Chargel(l)) 

End If 
Do 1=2, N2 

If (Chargel(l) == 0) Then 

Tmpl = Tmpl + DAbs(Chargeln(l) - (FDCon(l) / Dielec * Chargel(l)) - & 
(Field(l) * ((FDCon(l+1) - FDCon(l-l)) / (2 * DeltaX)))) 

Else 

Tmpl = Tmpl + DAbs((Chargeln(l) - (FDCon(l) / Dielec * Chargel(l)) - & 
(Field(l) * ((FDCon(l+1) - FDCon(l-l)) / (2 * DeltaX)))) * DeltaT / & 
Chargel (I)) 

End If 
End Do 

If (Chargel (N2+1)== 0) Then 

Tmpl = Tmpl + DAbs(Chargeln(N2+1) - (FDCon(N2+1) / Dielec * & 

Chargel (N2+1)) - (Field(N2+1) * ((FDCon(N2+1) - FDCon(N2)) / (DeltaX)))) 

Else 

Tmpl = Tmpl + DAbs((Chargeln(N2+1) - (FDCon(N2+1) / Dielec * & 

Chargel (N2+1)) - (Field(N2+1) * ((FDCon(N2+1) - FDCon(N2)) / & 
(DeltaX)))) * DeltaT/ Chargel (N2+1)) 

End If 
Else 

Do 1=1, N2+1 

If (Chargel (I) == 0) Then 

Tmpl = Tmpl + DAbs(Chargeln(l) - (Conduc / Dielec * Chargel(l))) 

Else 

Tmpl =Tmp1 + DAbs((Chargeln(l) - (Conduc/ Dielec * Chargel (I))) * & 
DeltaT / Chargel (I)) 

End If 
End Do 
End If 

Residual = Tmpl 
If (Residual <= Error) Then 
Again = .False. 

End If 

If (Count2 == PCount) Then 
! Store data 

Times(OCount) = Count * DeltaT 
Do 1=1, N2+1 

Charges(I.OCount) = Charge2(l) 

Fields(I.OCount) = Field(l) 

Voltages(l,OCount) = Voltage(l) 

If (HFC) Then 

HFCons(I.OCount) = HFCon(l) 

FDCons(l,OCount) = FDCon(l) 
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End If 
End Do 

OCount = OCount + 1 
Count2 = 0 
End If 

Print *, Count * DeltaT, Residual 

If (Again) Then 

Count = Count + 1 
Count2 = Count2 + 1 
Do 1=1, N2+1 

Chargel(l) = Charge2(l) 

End Do 
End If 
End Do 

! Store data 

Times(OCount) = Count * DeltaT 
Do 1=1, N2+1 

Charges(l, OCount) = Charge2(l) 

Fields(l, OCount) = Field(l) 

Voltages(l, OCount) = Voltage(l) 

If (HFC) Then 

HFCons(l, OCount) = HFCon(l) 

FDCons(l, OCount) = FDCon(l) 

End If 
End Do 
End If 

! Output data to output file 
Write (4,*) 

If (Break) Then 

Write (4,*) 'Breakdown occurred at (sec): ', BreakTime 
Else 

Write (4,*) 'Steady-state reached' 

End If 
Write (4,*) 

Write (4,*) 'Charge Density (#/m3)' 

Write (4,35) (Char(9), Times(J), J=1, OCount) 

Do 1=1, N2+1 

Write (4,45) Depth(l), (Char(9), Charges(U), J=1, OCount) 
End Do 
Write (4,*) 

Write (4,*) 'Electric Field (V/m)' 

Write (4,35) (Char(9), Times(J), J=1, OCount) 

Do 1=1, N2+1 

Write (4,45) Depth(l), (Char(9), Fields(I.J), J=1, OCount) 
End Do 
Write (4,*) 

Write (4,*) 'Voltage (V)' 

Write (4,35) (Char(9), Times(J), J=1, OCount) 

Do 1=1, N2+1 

Write (4,45) Depth(l), (Char(9), Voltages(I.J), J=1, OCount) 
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End Do 
If (HFC) Then 
Write (5,*) 

Write (5,*) 'High Field Conductivity Scaling Factor' 

Write (5,35) (Char(9), Times(J), J=1, OCount) 

Do 1=1, N2+1 

Write (5,45) Depth(l), (Char(9), HFCons(I.J), J=1, OCount) 
End Do 
Write (5,*) 

Write (5,*) 'Conductivity (1/ohm-m)' 

Write (5,35) (Char(9), Times(J), J=1, OCount) 

Do 1=1, N2+1 

Write (5,45) Depth(l), (Char(9), FDCons(I.J), J=1, OCount) 
End Do 
End If 

! Close output file 
Write (3,*) 'End of file’ 

Write (4,*) 'End of file' 

Close (Unit = 3) 

Close (Unit = 4) 

If (HFC) Then 

Write (5,*) 'End of file' 

Close (Unit = 5) 

End If 

End 


! Interp Subroutine 

! This subroutine performs a linear interpolation 

Subroutine lnterp(x, xl, x2, y 1 , y2, y) 

Double Precision x, xl, x2, y 1 , y2, y 

Y = ((x - xl) / (x2 - xl)) * (y2 - yl) + yl 

Return 

End 
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APPENDIX B 


CODE RUN PLOTS 


This appendix contains the output plots from CoDDCA for high field 
conductivity, radiation induced conductivity, and specific cases where the results differ 
significantly from the base case. The cases that were significantly different are those 
from the sensitivity studies of the thickness, density, geomagnetic substorms, and solar 
particle events. The plots included are the charge density distribution, electric field 
distribution, and voltage distribution as functions of depth into the material and time. 
The data used for inputs can be found by referring to Chapter 5, and Tables 5.1, 5.2, 5.4, 
5.5, 5.6, and 5.7. 
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Figure B. 1 : High field conductivity charge density comparison with base case 
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Figure B.2: High field conductivity electric field comparison with base case 
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Figure B.3: High field conductivity volt', 
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Figure B.4: Radiation induced conductivity charge density comparison with base case 
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Figure B.5. Radiation induced conductivity electric field comparison with base 
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Figure B.6: Radiation induced conductivi 
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Figure B.7: Electric field of base case with thickness of 0.01 mm 
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Figure B.9: 
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Figure B. 1 0: Electric field of base case with thickness of 3 mm 
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Figure B . 1 2 : Electric field of base case with thickness of 1 0 mm 
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Figure B.13: Electric field of base case wi 
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Figure B. 1 5: Electric field of base case with a substorm intensity of 100 
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Figure B. 1 6: Electric field of base case with a substorm intensity of 1 000 
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Figure B. 1 7. Electric field of base case with a solar particle event intensity of 5 
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Figure B . 1 8 : Electric field of base case with a solar particle event intensity of 1 0 
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Figure B.19. Electric field of base case with a solar particle intensity of 100 
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Figure B. 20: Electric field of base case with a solar particle event intensity of 1000 
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APPENDIX C 


EXPERIMENTAL DATA 


This appendix contain the conductivity and dielectric strength measurements, and 
the stress-strain and transverse versus longitudinal strain plots from the tensile tests. The 
slope of the linear portion of stress-strain curve is the Young’s modulus in terapascals, 
‘Stroke’ represents strain calculated from the stroke and ‘Gage’ represents strain obtained 
from a strain gage. The slope of the linear portion of the transverse versus longitudinal 
strain curve is the negative of the Poisson ratio. 

The numbering system used for conductivity measurement samples is ‘R 
the ‘R’ stands for resistance measurements, the first number is the layup number (‘ 1 ’ is 
for [0] 4 , ‘2’ is for [0] 8 , ‘3’ is for [0] 32 , ‘4’ is for [±45] s , ‘5’ is for [0/±45/90] s , and ‘6’ is for 
[0/±45/90] 4S ), the second number is the sample area number (‘1’ is for 645 mm 2 (1 in 2 ), 
‘2’ is for 1290 mm 2 (2 in 2 ), and ‘3’ is for 2580 mm 2 (4 in 2 )), and the last number is the 
sample number in the specified test group. 

The numbering system used for the dielectric strength measurement samples is 
‘B-#-#’; the ‘B’ stands for breakdown voltage measurements, the first number is the 
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layup number defined above, and the second number is the sample number for the 
specified test group. 


The numbering system used for the tensile test specimen is the first 

number is the cure set number (‘3’ is for [0] 4 0% carbon black, the rest of the numbers are 
for [±45] s layups, ‘4’ is for 0% carbon black, ‘5’ is for 5% carbon black, ‘6’ is for 10% 
carbon black, ‘7’ is for 15% carbon black, and ‘8’ is for 20% carbon black), the ‘4’ 
represent that the layup contains 4 plies, the third number represents the laminate number 
(3 laminates of the same layup were cured in one cure set), and the last letter represents 
the sample designator. 
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Table C. 1 : Unidirectional unsanded conductivity resistance data 


Sample 

Thickrn 

Average 

jss (mm) 
| SD 

Area (mm 2 ) 
Average | SD 

Resistance (Q) 
Average | SD 

R-1-1-1 

0.52 

0.007 

630.5 

3.03 

0.71 

0.049 

R-1-1-2 

0.52 

0.005 

624.0 

1.47 

1.43 

0.046 

R-1-1-3 

0.52 

0.005 

621.3 

3.58 

0.86 

0.115 

R-1-2-1 

0.52 

0.007 

1240.9 

3.88 

0.90 

0.199 

R-1-2-2 

0.52 

0.007 

1240.7 

5.47 

1.05 

0.010 

R-1-2-3 

0.52 

0.008 

1230.0 

0.94 

0.45 

0.087 

R-1-3-1 

0.50 

0.011 

2541.3 

23.77 

0.48 

0.010 

R-1-3-2 

0.51 

0.011 

2540.6 

5.26 

0.47 

0.182 

iMimsm 

0.51 

0.008 

2541.3 

16.66 

0.45 

0.035 

R-2-1-1 

0.99 

0.008 

644.4 

4.72 

0.70 

0.093 

R-2-1-2 

1.00 

0.005 

632.9 

13.72 

0.94 

0.176 

R-2-1-3 

0.96 

0.023 

643.6 

1.41 

1.09 

0.059 

R-2-2-1 

0.96 

0.019 

1277.9 

1.10 

0.98 

0.046 

R-2-2-2 

1.01 

6.009 

1270.6 

1.14 

0.81 

0.212 

R-2-2-3 

1.00 

0.004 

1276.1 

0.99 

0.70 

0.156 

R-2-3-1 

0.99 

0.016 

2567.9 

4.33 

0.58 

0.040 

R-2-3-2 

0.97 

0.025 

2585.7 

3.84 

0.58 

0.040 

R-2-3-3 

0.97 

0.036 

2577.9 

4.72 

0.43 

0.078 

R-3-1-1 

3.71 

0.005 

645.7 

1.91 

0.55 

0.051 

R-3-1-2 

3.67 

0.013 

647.0 

0.38 

2.22 

0.012 

R-3-1-3 

3.67 

0.013 

641.7 

2.16 

0.90 

0.055 

R-3-2-1 

3.71 

0.013 

1285.7 

5.55 

0.40 

0.023 

R-3-2-2 

3.70 

0.020 

1281.2 

2.33 

0.47 

0.040 


3.73 

0.005 

1287.7 

4.57 

0.63 

0.010 

R-3-3-1 

3.64 

0.047 

2583.0 

7.06 

0.40 

0.052 

R-3-3-2 

3.66 

0.041 

2574.2 

14.09 

0.34 

0.078 

R-3-3-3 

3.69 

0.029 

2585.0 

1.80 

0.39 

0.046 
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Table C.2: 


Quasi-isotropic unsanded conductivity resistance data 


Sample 

Thickness (mm) 
Average | SD 

Area (mm 2 ) 
Average [ SD 

Resista 

Average 

nee (Q) 
SD 

R - 4 - 1-1 

0.54 

0.007 

624.8 

5.63 

1.30 

0.141 

R - 4 - 1-2 

0.52 

0.009 

627.5 

2.71 

0.55 

0.051 

R - 4 - 1-3 

0.53 

0.005 

609.1 

1.25 

0.61 

0.096 

R - 4 - 2-1 

0.53 

0.008 

1269.9 

0.96 

0.42 

0.021 

R - 4 - 2-2 

0.52 

0.016 

1271.6 

2.80 

0.46 

0.010 

R - 4 - 2-3 

0.54 

0.005 

1267.9 

2.93 

0.42 

0.105 

R - 4 - 3-1 

0.52 

0.009 

2558.5 

3.25 

0.34 

0.139 

R - 4 - 3-2 

0.51 

0.011 

2560.3 

6.14 

0.26 

0.017 

R - 4 - 3-3 

0.52 

0.007 

2549.6 

2.95 

0.27 

0.006 

R - 5 - 1-1 

0.98 

0.011 

625.5 

0.72 

0.67 

0.065 

R - 5 - 1-2 

0.94 

0.015 

626.3 

1.77 

0.75 

0.159 

R - 5 - 1-3 

1.00 

0.013 

627.0 

1.74 

0.58 

0.051 

R - 5 - 2-1 

1.01 

0.025 

1269.6 

1.18 

0.39 

0.049 

R - 5 - 2-2 

1.03 

0.008 

1272.3 

0.66 

0.54 

0.020 

R - 5 - 2-3 

0.99 

0.028 

1267.9 

8.68 

0.50 

0.021 

R - 5 - 3-1 

0.99 

0.033 

2539.3 

25.02 

0.33 

0.040 

R - 5 - 3-2 

0.97 

0.027 

2565.1 

3.05 

0.32 

0.031 

R - 5 - 3-3 

1.01 

0.011 

2556.2 

2.74 

0.29 

0.006 

R - 6 - 1-1 

3.83 

0.021 

647.8 

3.46 

0.66 

0.087 

R - 6 - 1-2 

3.84 

0.011 

648.3 

4.05 

0.70 

0.015 

R - 6 - 1-3 

3.88 

0.007 

648.3 

5.28 

0.53 

0.026 

R - 6 - 2-1 

3.88 

0.008 

1272.0 

2.39 

0.28 

0.010 

R - 6 - 2-2 

3.78 

0.075 

1276.4 

3.39 

0.31 

0.006 

R - 6 - 2-3 

3.79 

0.068 

1282.1 

4.59 

0.31 

0.000 

R - 6 - 3-1 

3.74 

0.072 

2513.0 

4.60 

0.31 

0.021 

R - 6 - 3-2 

3.82 

0.033 

2542.3 

3.92 

0.32 

0.015 

R - 6 - 3-3 

3.73 

0.069 

2567.1 

5.85 

0.28 

0.000 
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Table C.3: Unidirectional sanded conductivity resistance data 


Sample 

Thickne 

Average 

ss (mm) 
SD 

Area 

Average 

mm 2 ) 

SD 

Resist: 

Average 

ince (Q) 
SD 

R - 1 - 1-1 

0.39 

0.031 

630.5 

3.03 

1.49 

0.179 

R - 1 - 1-2 

0.37 

0.050 

624.0 

1.47 

0.33 

0.067 

R - 1 - 1-3 

0.34 

0.048 

621.3 

3.58 

0.30 

0.110 

R - 1 - 2-1 

0.39 

0.029 

1240.9 

3.88 

0.16 

0.062 

R - 1 - 2-2 

0.41 

0.046 

1240.7 

5.47 

0.15 

0.029 

R - 1 - 2-3 

0.34 

0.045 

1230.0 

0.94 

0.17 

0.068 

R - 1 - 3-1 

0.29 

0.054 

2541.3 

23.77 

0.24 

0.120 

R - 1 - 3-2 

0.36 

0.051 

2540.6 

5.26 

0.30 

0.108 

R - 1 - 3-3 

0.34 

0.055 

2541.3 

16.66 

0.30 

0.035 

R - 2 - 1-1 

0.78 

0.048 

644.4 

4.72 

0.20 

0.078 

R - 2 - 1-2 

0.80 

0.017 

632.9 

13.72 

0.36 

0.086 

llrMBC 

0.77 

0.020 

643.6 

1.41 

0.15 

0.020 

R - 2 - 2-1 

0.82 

0.028 

1277.9 

1.10 

0.32 

0.156 

R - 2 - 2-2 

0.82 

0.045 

1270.6 

1.14 

0.28 

0.080 

R - 2 - 2-3 

0.78 

0.045 

1276.1 

0.99 

0.14 

0.099 

R - 2 - 3-1 

0.85 

0.041 

2567.9 

4.33 

0.17 

0.107 

R - 2 - 3-2 

0.84 

0.053 

2585.7 

3.84 

0.16 

0.119 

R - 2 - 3-3 

0.82 

0.054 

2577.9 

4.72 

0.07 

0.055 

R - 3 - 1-1 

3.41 

0.179 

645.7 

1.91 

0.92 

0.284 

R - 3 - 1-2 

3.56 

0.042 

647.0 

0.38 

0.81 

0.425 

R - 3 - 1-3 

3.57 

0.031 

641.7 

2.16 

0.41 

0.099 

R - 3 - 2-1 

3.58 

0.064 

1285.7 

5.55 

0.25 

0.155 

R - 3 - 2-2 

3.56 

0.052 

1281.2 

2.33 

0.20 

0.085 

R - 3 - 2-3 

3.60 

0.047 

1287.7 

4.57 

0.15 

0.031 

R - 3 - 3-1 

3.55 

0.062 

2583.0 

7.06 

0.14 

0.070 

R - 3 - 3-2 

3.56 

0.040 

2574.2 

14.09 

0.01 

0.006 

R - 3 - 3-3 

3.56 

0.036 

2585.0 

1.80 

0.06 

0.031 
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Table C.4: Quasi-isotropic sanded conductivity resistance data 


Sample 


SI 




R - 4 - 1-1 

0.45 

0.023 

624.8 

5.63 

0.38 

0.085 

R - 4 - 1-2 

0.41 

0.022 

627.5 

2.71 

0.41 

0.275 

R - 4 - 1-3 

0.44 

0.016 

609.1 

1.25 

0.21 

0.035 

R - 4 - 2-1 

0.38 

0.029 

1269.9 

0.96 

0.25 

0.095 

R - 4 - 2-2 

0.38 

0.047 

1271.6 

2.80 

0.23 

0.057 

R - 4 - 2-3 

0.39 

0.015 

1267.9 

2.93 

0.13 

0.031 

R - 4 - 3-1 

0.41 

0.048 

2558.5 

3.25 

0.11 

0.015 

R - 4 - 3-2 ] 

0.39 

0.029 

2560.3 

6.14 

0.18 

0.099 

R - 4 - 3-3 j 

0.42 | 

0.024 "1 

2549.6 

2.95 

0.15 

0.035 

R - 5 - 1-1 

0.89 

0.012 

625.5 

0.72 

0.14 

0.087 

R - 5 - 1-2 

0.81 

0.042 

626.3 

1.77 

0.17 

0.006 

R - 5 - 1-3 

0.90 

0.021 

627.0 

1.74 

0.07 

0.025 

R - 5 - 2-1 

0.87 

0.029 

1269.6 

1.18 

0.07 

0.010 

R - 5 - 2-2 

0.88 

0.040 

1272.3 

0.66 

0.09 

0.026 

R - 5 - 2-3 

0.91 

0.022 

1267.9 

8.68 

0.13 

0.010 

R - 5 - 3-1 

0.87 

0.044 

2539.3 

25.02 

0.09 

0.059 

R - 5 - 3-2 

0.80 

0.070 

2565.1 

3.05 

0.03 

0.010 

R - 5 - 3-3 

0.87 

0.048 

2556.2 

2.74 

0.08 

0.006 

R - 6 - 1-1 

3.72 

0.023 

647.8 

3.46 

022 

0.021 

R - 6 - 1-2 

3.76 

0.018 

648.3 

4.05 

0.20 

0.021 

R - 6 - 1-3 

3.80 

0.025 

648.3 

5.28 

0.15 

0.075 

R - 6 - 2-1 

3.79 

0.042 

1272.0 

2.39 

0.04 

0.012 

R - 6 - 2-2 

3.68 

0.061 

1276.4 

3.39 

0.10 

0.012 

R - 6 - 2-3 

3.69 

0.054 

1282.1 

4.59 

0.11 

0.021 

R - 6 - 3-1 

3.61 

0.073 

2513.0 

4.60 

0.07 

0.031 

R - 6 - 3-2 

3.73 

0.037 

2542.3 

3.92 

0.02 

0.015 

R - 6 - 3-3 

3.64 

0.072 

2567.1 

5.85 

0.04 

0.021 
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Table C.5: Dielectric strength breakdown voltage data 


Sample 

Thickne 

Average 

ss (mm) 
SD 

Breakdown 

Average 

Voltage (V) 
SD 

B-1-1 

0.51 

0.013 

1100 

446.8 

B-1-2 

0.51 

0.013 

810 

424.6 

B-1-3 

0.50 

0.013 

817 

127.0 

B-2-1 

1.01 

0.000 

1533 

57.7 

B-2-2 

1.01 

0.004 

1400 

264.6 

B-2-3 

0.94 

0.025 

1367 

472.6 

B-3-1 

3.70 

0.011 

1000 

264.6 

B-3-2 

3.71 

0.005 

1100 

100.0 

B-3-3 

3.68 

0.015 

1020 

170.9 

B-4-1 

0.54 

0.004 

830 

270.0 

B-4-2 

0.53 

0.004 

1200 

173.2 

B-4-3 

0.55 

0.005 

1700 

556.8 

B-5-1 

0.96 

0.020 

1767 

378.6 

B-5-2 

1.01 

0.007 

1767 

503.3 

B-5-3 

1.02 

0.005 

1533 

404.1 

B-6-1 

3.73 

0.033 

1157 

567.0 

B-6-2 

3.71 

0.033 

953 

331.3 

B-6-3 

3.80 

0.023 

897 

179.5 
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Transverse Strain (|is/e) 




Figure C.l: Tensile test 3-4-1 -a 
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Figure C.4: Tensile test 3-4-2-b 
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Figure C.5: Tensile test 3-4-3-a 
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Figure C. 7: Tensile test 4-4-1 -b 



Transverse Strain (jae/s) Stress (MPa) 




Figure C.8: Tensile test 4-4- 1-c 
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Figure C.9: Tensile test 4-4-2-a 
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Figure C. 1 0: Tensile test 4-4-2-b 
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Figure C. 1 4: Tensile test 5-4- 1 -b 
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Figure C. 16: Tensile test 5-4-2-b 
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Figure C. 18: Tensile test 5-4-3-b 
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Transverse Strain (\xelt) Stress (MPa) 




Figure C. 19: Tensile test 6-4-1 -a 
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Figure C.20: Tensile test 6-4-1 -b 
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Figure C.2 1 : Tensile test 6-4-2-a 
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Transverse Strain (^s/e) Stress (MPa) 




Figure C.22: Tensile test 6-4-2-b 
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Figure C.23: Tensile test 6-4-3-a 
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Transverse Strain (j lele) Stress (MPa) 




Figure C.24: Tensile test 6-4-3-b 
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Longitudinal Strain (j iele ) 


Figure C. 26: Tensile test 7-4-1 -b 
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Figure C.28: Tensile test 7-4-2-b 
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Longitudinal Strain (ne/e) 


Figure C.29: Tensile test 7-4-3-a 
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Figure C. 31: Tensile test 8-4-1 -a 
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Longitudinal Strain (\izle) 


Figure C.33: Tensile test 8-4-2-a 
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Figure C.34: Tensile test 8-4-2-c 
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Figure C.35: Tensile test 8-4-3-a 
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